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THE GEOMAGNETIC LATITUDE EFFECT THE NUCLEON 
AND MESON COMPONENT COSMIC RAYS SEA 


ABSTRACT 


Measurements have been taken the changes intensity the nucleon 
and meson components cosmic rays during cruise the Canadian Naval 
Icebreaker Labrador into the Arctic, through the North West Passage, and 
circumnavigating the North American Continent. The geomagnetic latitudes 
covered extend from 18°N. 89°N. The latitude knee clearly shown 
geomagnetic latitude about 52° the case the nucleon component and less 
definitely between 40° and 50° the case the meson component. The rigidity 
particles arriving vertical direction 52° 2.1 Bv. and 45° 3.7 Bv. 
Meyer and Simpson have shown that changes the primary spectrum between 
1948 and 1954 probably extend these rigidities and such changes should, 
therefore, sea level. The longitude effect low latitudes 
clearly shown differences intensity between the measurements the east 
and west sides North America. the case the meson component, the 
magnitude the longitude effect these longitudes was found greater 
than that shown Millikan and Neher 1936. The interpretation the 
meson component results above the knee complicated difficulties 
temperature correction. the case the nucleon component, apparent 
longitude effect exists above the knee that there was small difference the 
intensity high latitudes the eastern and western parts the North 
American Arctic. satisfactory explanation offered for this. The diurnal 
variation the nucleon component high latitudes shown but unusual 
features were found. Appreciation expressed the Royal Canadian Navy 
for making these measurements possible. 


INTRODUCTION 


The interaction between primary cosmic rays and the earth’s magnetic 
field explained the theory allowed and forbidden 
orbits. Detailed agreement between theory and measurements complicated 
the fact that the earth’s field not quite equivalent centered dipole 
and because the effects the atmosphere make extensive direct measurements 
primary cosmic rays virtually impossible. Extensive measurements covering 
appreciable part the earth’s surface involve the organization large 
expeditions and therefore have not very often been undertaken. 

The new measurements described here were taken sea level board 
the Canadian Navy Icebreaker Labrador during its 1954 cruise through the 
North West Passage and circumnavigating the North American Continent. 
Both the nucleon and the meson component were measured. The range 
latitudes covered (from 18° 89° geomagnetic latitude) includes the so- 
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called knee the intensity vs. geomagnetic latitude curves, that the point 
above which the intensity does not increase increases very slowly. 

The position the knee both sea level and high altitude the atmos- 
phere very important analyzing cosmic ray data. Neher, Peterson, and 
Stern (1953) using measurements taken near the top the atmosphere 
1951 show that the knee close geomagnetic latitude 58°N. with 
reference the centered dipole magnetic field. Van Allen (1953) and Ellis, 
Gottlieb, and Van Allen (1954) taking measurements outside the atmosphere 
the same epoch confirm cutoff the primary spectrum. More recent 
measurements Neher and Stern (1955) show that the position the knee 
near the top the atmosphere not fixed, being much further north 1954 
than 1951. That this shift the knee extends well down into the atmosphere 
was shown Simpson (1951), Simpson, Fonger, and Treiman (1953), and 
Meyer and Simpson (1955), taking measurements airplane altitudes. The 
change the position the knee between 1948 and 1954 was quite appre- 
ciable and represents estimated change total intensity much 13%. 

The position the knee low altitudes gives some information about the 
minimum energy primary particles that can contribute through the atmos- 
phere the observed effects. Since much easier take extensive measure- 
ments ground sea level than higher the atmosphere, accurate 
knowledge the latitude effect and its variations still sufficiently important 
justify extensive expeditions. Further, there the fact that, though 
assumed that the interaction the earth’s magnetic field primary cosmic 
rays effective far from the earth’s surface that the field may treated 
dipole, the effective position the dipole the extent which this assump- 
tion true not all certain. taking early measurements the latitude 
effect, Compton (1933) pointed out that the north and south knees constant 
intensity correlate better with geomagnetic latitude than with geographic 
and has been the custom since then relate variations intensity cosmic 
rays with the geomagnetic axis the axis the nearest equivalent uniformly 
magnetized sphere. The longitude effect due the eccentricity the nearest 
equivalent dipole well known but the analysis Millikan and Neher (1936) 
shows that the cosmic ray data such were then available are not agree- 
ment with the position the dipole indicated magnetic analyses (Chapman 
and Bartels 1940), the minimum intensity the equator being approxi- 
mately 90° longitude rather than 161.8° estimated Chapman and 
Bartels (1940). 

Probably the most comprehensive previous study the latitude effect 
sea level that Compton and Turner (1937) who measured the hard com- 
ponent during several trips between Vancouver and Sydney. The latitude 
range covered between 40° and 52° geomagnetic and indicates knee 
between 35° and 40° though some earlier work placed little higher 
latitude (Clay 1933). The effect the atmosphere makes interpretation 
measurements the hard meson component high latitude difficult since 
this component subject the well-known temperature effect, the corrections 
depending the temperature throughout the atmosphere and for which there 
simple practical method correcting large body data. 
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The nucleon component measured neutron production 
believed not subject changes due temperature variation 
the atmosphere. The geomagnetic effect for this component has been ob- 
served airplane altitudes Simpson (1951), and Simpson, Fonger, and 
Treiman (1953) have also shown that the temperature effect this component 
zero very small. the present experiment both the nucleon component 
and the meson component were recorded and the advantages working with 
the nucleon component because its lack temperature effect are very 
apparent. 

During the period the measurements described here were taken (July 
November 21, 1954), neutron monitor and large counter telescope 
were also kept continuous operation Ottawa. The neutron monitor 
particular showed evidence appreciable changes intensity the type 
usually associated with solar activity, that variations intensity with 
27-day recurrence interval. The sun during that period apparently was un- 
usually quiet. The results, therefore, should represent background without 
disturbances the type studied Simpson (1954). 


APPARATUS AND EXPERIMENTAL TECHNIQUE 


(1) The Nucleon Component 

The nucleon component was measured using the technique described 
Simpson, Fonger, and Treiman (1953) wherein the nucleon component, 
believed protons and neutrons, detected the production neutrons 
lead. The neutrons are slowed down paraffin and detected boron tri- 
fluoride proportional counters. Fig. shows schematic diagram the 


PARAFFIN 


C-BF, PROPORTIONAL COUNTERS 
2-1/2" DIAM. 
SENSITIVE LENGTH 


Fic. Schematic arrangement the neutron monitor. 


monitor” used this experiment. Two BF; (enriched with 
counters were used, their diameters being in. and the effective length being 
in. The counters were mounted tubular opening the center alu- 
minum tanks in. long and in. square. These two tanks were surrounded 
lead for length ft. the center along the axis the pile, the lead being 
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sheets in. thick. There were in. lead between the central paraffin tanks 
carrying the counters, in. top, in. underneath, and in. the sides. 
Surrounding this outer layer paraffin in. thick. The outer paraffin was 
also carried aluminum tanks in. long built that there were two tanks 
top, two underneath, and one each side. Since the lead occupied only 
in. the 45-in. length the pile, the spaces the end the lead were 
filled with paraffin again suitably shaped aluminum tanks. The aluminum 
tanks containing the paraffin were welded from 1/8-in. aluminum sheet 
(the smaller tanks from 1/16-in. sheet). The whole pile when fitted together 
had over-all dimensions in. long, in. wide, and in. high. This was 
housed steel frame, the sides which could drawn together clamp 
the whole pile solidly. The steel frame was then bolted suitable supports 
welded the ship’s deck. 

The use aluminum tanks hold the paraffin was adopted make the 
structure reasonably free from danger due fire and keep the paraffin from 
flowing when operating summer temperatures. The tank structure was readily 
demountable, the largest tank being in., and could handled 
easily two men. aluminum was chosen was the purest commercial 
grade and hence would represent the best construction for minimum absorp- 
tion slow neutrons. 

The amplifier and recording equipment was follows. Each neutron counter 
had preamplifier located just outside the pile the head the counter. The 
preamplifier was conventional three-tube ring with feedback and needs 
further description here. The preamplifiers had voltage amplification factor 
about 65. These were followed the main amplifier which consisted 
second ring three tubes with the same voltage amplification. Operating the 
counters 2500 volts gave output pulses due the boron plus neutron 
alpha particles about fifteen volts. The amplifiers were designed accept 
pulse with time rise about 0.5 usec. and had clipping time about 
usec. The output the amplifier was fed into biased diode discriminator 
and pulse shaping circuit. From this the pulses were fed into binary scaling 
circuits six stages making scaling factor 64. The output from the scalers 
was fed 20-pen Esterline Angus operations recorder. The two outer pens, 
Nos. and 20, were fed minute pulses from precision chronometer. The 
Esterline Angus tape was operated one foot per hour. Message registers were 
connected parallel with the 20-pen recorder. 

The counting rate sea level northern latitudes averaged about pips 
the tape per hour, that counts per hour from each neutron counter. 
The hourly values were usually read counting the pips for each hour the 
tape. The interval between pips the hour was estimated the nearest 10th 
the distance between pips making measurement error not more than 
+0.2% for hourly values. Since the hourly counting rate gave standard 
deviation 1.4% for each counter when the two were added, the mea- 
surement error was, therefore, trivial. For 12-hr. daily averages, this error 
would even more negligible. The message registers were read every hr. 
could have been read oftener anything had gone wrong with the recorder. 
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This method recording has been found very useful the cosmic ray 
observatory Ottawa. Counting the pips the Esterline Angus tape some- 
what tedious but operating the recorder continuously 0.2 in. per min. and 
adjusting scaling factors give between and 100 pips per hour gives very 
good record what happening. power troubles cause the equipment 
off for few minutes, even seconds, intense line transients produce 
extra counts thing that almost impossible eliminate completely 
shipboard though they were never more than few day), always apparent 
the Esterline Angus record. The bad minutes can dropped and the statis- 
tical error for that hour not increased appreciably. Further, intensity 
fluctuations such solar flare phenomena occur, the recorder tape gives 
great deal detail. Averaging times down minute can used, course, 
with corresponding loss accuracy. 

the Ottawa observatory where similar neutron monitor kept 
continuous operation, additional scaler added each channel, recording 
additional set pips four times fewer number but matching every fourth 
pip. This eases the counting the hourly values and have found this 
method simpler than photographic recording message registers. The neutron 
counters had been thoroughly tested and were tested and calibrated periodi- 
cally throughout the voyage. cadmium sheath could slipped the tube 
carrying the counters surround the counters. The sheath consisted two 
layers mils cadmium sheet which would reduce the count neutrons 
negligible numbers the neutrons were all the energy band below the 
absorption edge for cadmium. Actually tests cosmic-ray-produced neutrons 
different fluxes moderated paraffin this type experiment show that 
the reduction factor more like 40. The fact that boron has 
appreciable cross-section for neutrons energies above the cadmium absorp- 
tion energy makes this quite understandable. The actual cadmium background 
high latitudes was 140+13 counts per hour but the majority these were 
those transmitted the cadmium and absorbed the boron. The true 
background believed not more than counts per hour and 
consistent. This represents only 0.25% the counting rate. Neglecting 
would not affect any the results discussed here and, therefore, 
correction was made for background. 

Tests were also made periodically for the flatness the plateau counting 
rate with varying the high voltage the neutron counters. The plateau was 
not flat had been hoped, variation the high voltage causing 
average variation 0.9% the counting rate. The regulation the high volt- 
age was accomplished using Sorensen electronic regulator the a-c. 
side the high voltage power pack. This method had been found very satis- 
factory Ottawa and involved complicated high voltage electronic 
regulators. accurate Weston voltmeter was kept the high voltage circuit 
continuously and the variation observed many times per day. The variations 
noted were practically always less than +0.25%. 

Tests were also made for variation with biasing the pulse height dis- 
criminator. The mean pulse height was about fifteen volts. The noise level was 
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less than one volt. Varying the discriminator 0.1 volt made 0.16% change 
the counting rate and regular checking the bias voltage indicated that 
did not change more than +0.1 volt. 

The flatness these plateaus for neutron counters this sort can shown 
function the amplifier design. relatively slow amplifier with long 
clipping time gives much better plateau than fast amplifier. rise time 
long one microsecond and clipping time twenty more microseconds 
would have given better plateaus. Some experiments carried out with faster 
amplifiers having rise time 0.1 better and short clipping times 
(0.5 demonstrate this effect very clearly. For such fast amplifiers 
and short clipping times the pulse shape depends considerably where the 
particle occurs the counter, and wide range shapes can seen 
fast oscilloscope. The present amplifiers were too fast give the best plateaus 
but were quite satisfactory for the purpose and very simple design. 

For plateau tests and calibration neutron source was available. This 
introduced some difficulties the shipboard monitor since strong source 
would have been troublesome cramped quarters and very weak source 
would not give high counting rate. The spontaneous fission uranium metals 
has been found useful steady source neutrons and for shipboard 
calibration box uranium strips* bound aluminum case was occasion- 
ally put hole the monitor pile cut the lead and paraffin for that pur- 
pose. The calibration rate was little too slow for good results, being neces- 
sary leave the source for six hours get 0.3% statistics. The actual 
calibrations were not successful expected since the method was planned 
only just before the equipment was installed, and was not properly tested. 
The calibration was very sensitive the position the uranium source which 
could not sufficiently accurately maintained during different calibrations. 
Calibrations every two three weeks showed standard deviation +0.9% 
rather than the 0.3% statistics and believed the stability the counters 
and circuitry was much better than this. attempt was made correct for 
individual calibrations. 


(2) The Meson Component 


The meson component was observed triple coincidence telescope. Each 
the three trays Geiger counters contained eight counters in. diameter 
and in. effective length. The telescope was made very compact, the top two 
trays were separated only in. with shielding between them except the 
aluminum case carrying the counters. The second and third trays were sepa- 
rated in., there being in. lead shielding between them. was considered 
first that cubical design telescope should adopted such has been 
recommended for cosmic ray observations for the international geophysical 
year. Such telescope has been operation Ottawa for some time but 
sea and remote stations has been decided keep the equipment simple 
possible that will operate reliably with minimum testing. Using 


*The uranium were loaned Atomic Energy Canada Limited and wish 
ewis for arranging this loan. 
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the eight counters per tray compact telescope, the triple coincidence 
counting rate was close 14,000 per hour northern latitudes. 

Each counter tray contained preamplifier and cathode follower stage 
the end the tray. The high voltage circuit was split that either four the 
eight counters could have their voltage taken off without affecting the capacity 
the coupling circuit the preamplifier. This gave simple method testing 
because when scaler operated directly from counter tray, the difference 
between three counters (if one had failed) and four could detected minute 
so, whereas the difference between seven and eight would have required 
much more careful test. The operation the whole telescope could checked 
this way few minutes. 

The pulses from the trays were fed pulse shaping circuit and coinci- 
dence circuit which has been used successfully for several years the Ottawa 
observatory. slightly unconventional that coincidences are formed 
adding the biased grid the mixing tube rather than the more usual 
Rossi circuit. The resolving time was about give thorough check 
the operation the circuits, the coincidence channels and scalers were 
duplicated. have found this duplication very useful detecting 
trouble. Almost any trouble the pulse shaping coincidence part the 
circuit will put the two recording channels out synchronism and glance 
the recorder chart shows this up. This value because normally this 
equipment the observatory Ottawa runs continuously for months with 
minimum routine checking and when large quantity such equipment 
kept continuous operation, necessary keep the time required for 
testing minimum. 

Triple coincidences were recorded using scaling factor 256, the output 
being put the Esterline Angus 20-pen operations recorder. Double coinci- 
dences were recorded between trays and using scaling factor 1024. 
The counting rate for doubles was about 40,000 per hour. This included the 
total component except for the shielding the roof the hut housing the 
equipment. The double telescope also almost omnidirectional since the two 
trays were very close together. Correction for the total component very 
uncertain, and though great deal data have been collected, little use 
has been made them and they will not reported here. 


(3) Location the Equipment 

The neutron monitor and the counter telescope were housed hut about 
ft. square built into the base the mast H.M.C.S. Labrador. This was the 
best position available, being the highest position the ship and free from 
changes surrounding material that might produce secondary neutrons, such 
moving cargo boats. The hut was heated and insulated and though 
there was automatic temperature control the temperature was kept 
80+10°F. The temperature the neutron monitor would change only very 
slowly since weighed all over 3000 Ib. There was evidence any effect 
due temperature changes. 

The Geiger counter telescope was housed special box which was insulated 
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and temperature controlled. The temperature was kept 90°F. Coincidence 
counting this sort has been found quite insensitive the telescope 
temperature provided does not low enough for the alcohol the counters 
start condense. Our experience that these counters not change their 
rate sufficiently affect coincidences until the temperature below about 
55°F. The failure thermostat one time brought the temperature the 
box 110°F. with apparent effect the coincidence rate. 

The shielding the roof the hut was in. cork and in. steel. The 
mast was open tower structure and carried several radar antennae. This 
made the exact shielding rather uncertain but since all the equipment was 
either omnidirectional accepted very large solid angle, the effect the 
mast structure and the radar antennae overhead was probably very small. 

For comparison with the results during the voyage around the North 
American continent identical neutron monitor was kept operation 
Ottawa, the shielding over being the same order and the amplifiers and 
recording equipment identical. 

large cubical counter telescope was also operated the cosmic ray 
observatory Ottawa during the time the voyage. 

Resolute (lat. 74° 41’ N., long. 94° 54’ W., geomag. lat. 83° N.) small 
triple coincidence telescope had been operation since 1949. This was replaced 
August larger telescope identical with that aboard the 
described above. The results the Resolute observations during the Arctic 
part the voyage will discussed later this paper. 


(4) Meteorological Data 


The Labrador carried full surface meteorological station but did not have 
sonde equipment for upper atmosphere measurements. Pressure and tempera- 
ture data were supplied the meteorological office. 

accurate measure pressure aboard ship difficult. There were 
several barographs and mercury barometer. barograph was located the 
cosmic ray hut and checked regularly with the mercury barometer the 
meteorological office. The meteorological office was below the main deck 
the ship but sufficiently near open hatches that ventilation pressures should 
not affect the readings. When the cosmic ray hut was closed cold weather 
whistling winds the ventilators often caused barometer movements 
fraction millibar. Rolling the instrument relying partially 
gravity balance inaccurate even well damped. Carefully comparing 
hourly readings the mercury barometer with the barograph the cosmic 
ray hut indicated that pressure effects due wind could not eliminated 
and variable differences between the two positions fraction millibar 
and occasionally high millibar which must due eddy pressures 
the wind were apparent. For hourly values the mercury barometer was taken 
standard except when there appeared obvious occasional errors 
when hourly readings were missed. The daily means were taken from the 
average the hourly readings. Although the barometric data were kept 
the nearest 10th millibar, because the above difficulties felt that 
accuracy greater than +0.3 millibars cannot claimed. This problem 
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that will always occur shipboard experiments this type. This makes 
small contribution the statistical error. 

Since upper atmosphere temperatures were not available and generally 
accepted that surface temperature not satisfactory variable correcting 
for meteorological effects the meson component, detailed record the 
temperature was not kept with the cosmic ray records. was, course, kept 
the meteorological office and hourly values and thermograph records are 
available. The temperature effect the meson component will discussed 
after the results are presented. temperature effect was expected the case 
the nucleon component. 


(5) Latitude and Longitude 

The position the ship was reported every three hours the meteorological 
office and this record was used for the cosmic ray results. Navigation problems 
the Arctic are difficult because charts are imperfect, and often the sun was 
not visible for days. Navigation close waters was visual radar 
shore prominences which may may not correctly charted. Further, while 
icebreaking the dead-reckoning method navigation useless. The posi- 
tion the ship was kept well possibie and for the most part was known 
within two three miles but sometimes this uncertainty would high 
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Fic. Map North America showing the course where cosmic ray measurements were 
taken. 
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ten twenty miles. This represents about minutes latitude but may 
degree more longitude high latitudes. 

When travelling north south open water cruising speeds, the ship 
covered about latitude per day. analyzing the results for the latitude 
effect, therefore, 12-hr. averages were used. The position latitudes and longi- 
tudes the three-hourly records were averaged give mean latitude and 
longitude for the 12-hr. period. Fig. shows map North America and the 
course the ship. The voyage started from Halifax July 23, 1954. The 
highest latitude was reached August and the ship stayed that area for 
four days. Considerable time was spent the North West Passage 
latitude about 75° and between longitudes 85° and 100° Some time 
was also spent the Chukchee Sea area north Alaska. The cruise reached 
Victoria September 28, Panama November Grenada November 10, 
and returned Halifax November 21, 1954. 


III. RESULTS 


(1) The Nucleon Component 

Fig. shows the results the measurements the nucleon component 
plotted day day throughout the voyage. Corrections have been made for 
pressure using barometer coefficient —9.9% per cm. mercury —0.74% 
per millibar. This corresponds absorption lengths 136 This was 


140 
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90 30 10 20 30 Te) 20 30 10 20 30 0 20 

AUG. SEPT. OcT. NOV. 


Fic. The nucleon component. Twelve-hour averages throughout the voyage. 


the actual coefficient found during the northern part the expedition, that 
from July lat. 67° N., long. 55° September 23, lat. 57° N., long. 
167° linear regression was calculated between pressure and intensity 
giving the value —0.744+0.007% per mb. For the Ottawa neutron data the 
barometer coefficient evaluated for each month. The mean during the months 
August, September, and October, 1954, was —0.76% per mb. previous 
set measurements six-month period 1952-53 time low solar 
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activity gave value —0.72% per mb. Simpson, Fonger, and Treiman 
(1953) also give value —0.72% per mb. but point out the difficulty 
getting good barometer when other intensity variations are present. 
The monthly values Ottawa show much greater variation than the 
statistical variance derived from each month’s analyses. 

The method used making the barometer corrections merits some discus- 
sion. Using linear expression the correction should 


being the mean intensity over the whole period observations. The actual 
absorption is, however, closely exponential which case for small variations 
the correction 


(2) 


can used, being the intensity observed during the immediate period for 
which the correction made. For larger variation pressure, the integral 
value should used and 


(3) 


where the subscript again refers means over the whole period the 
experiment. Equation (2) first approximation equation (3) and this 
method correction was used, being the coefficient found linear regression 
during the part the cruise north the latitude knee. Using equation (1) 
would lead difficulties below the knee where the intensity change due 
latitude was rapid. The difference between (2) and the more accurate ex- 
pression (3) was considered too small justify the computation necessary 
use (3). 

attempt was made derive barometer coefficient the tropics during 
the present voyage but was complicated the rapid changes due latitude. 
Picking few days nearly equal mean intensity near Panama and using 
three-hour averages produced barometer coefficient per 
mb. This not significantly different from that the north but very low 
accuracy partly because the extremes barometer variation that time were 
only few millibars. 

should noted Fig. that the intensities Halifax July and 
November are significantly different. This may seasonal annual 
change intensity and consistent with changes the Ottawa value the 
same time. The changes the Ottawa value were taken represent world- 
wide annual variation and the shipboard results were corrected proportion- 
ately. This brought the Halifax values the beginning and end the voyage 
together and the corrected results are shown Fig. the intensity being 
plotted against geomagnetic latitude. 

The geomagnetic latitude for each 12-hr. average was estimated using 
McNish (1936) curves, but correcting for change the position the 
geomagnetic axis given the more recent calculations Vestine (1947). 
McNish’s curves were based pole lat. 78.5° N., long. 69.0° W., while 
the later analysis gives lat. 78.6° N., long. 70.1° The latitude correction 
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Fic. nucleon component. Intensity plotted function geomagnetic latitude. 
The solid circles represent measurements the eastern side the North American Con- 
tinent and crosses, the western side. 


was not used since was small and would involve considerable computation, 
but the longitude correction was readily made curves. few 
actual calculations showed that his method was well within the accuracy 
known positions the ship. Fig. the solid circles represent the east side 
the continent and the crosses the west. The east side gives truer picture 
the geomagnetic latitude effect since the curve from Grenada right 
within degree the geomagnetic pole was almost one longitude circle, 
deviating only few degrees. 

Each point the curve represents minimum hr., that standard 
deviation due counting variations Allowing for errors 0.2% due 
voltage variation, and 0.1% due difficulties getting good pressure 
readings (averaged over hr., this probably fair estimate the error) 
over-all standard deviation 0.4% quite generous. The value 
shown Fig. the counting rate error only. For points where the ship was 
anchor for some days when the geomagnetic latitude was changing only 
very slowly, longer averages were used, otherwise some parts the curve 
would too cluttered with points intelligible. The counting rate error 
for such points is, course, smaller than that shown. low latitudes where 
the longitude was changing rapidly, the curve not completed. 

The important features this result are: (1) the position the latitude 
sea level, (2) the demonstration the longitude effect, and (3) the 
difference the curves high latitude representing the voyage north and the 
return south. These points will discussed separately. The position the 
latitude knee appears further north than the case the meson com- 
ponent. Its position outside the earth’s atmosphere during 1951 was shown 
Neher, Peterson, and Stern (1953) about 58°. 30,000 ft., Meyer and 
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Simpson (1955) show the knee about 55° 1954 though was low 
50° 1948. The present result places the knee sea level around 52° but 
very difficult define sharply. The rigidity particles for these 
latitudes according the Stormer theory such that the minimum energies 
for protons arriving vertically are 0.6 Bev. for 58°, 1.0 Bev. for 55°, and 1.5 
Bev. for 52° little less than half that value per nucleon for particles mass 
2Z. one can assume that the yield neutrons sea level closely related 
the number vertical primaries, these figures indicate that appreciable 
portion primary particles with energies low 1.5 Bev. can contribute 
neutron production sea level. points where the geomagnetic cutoff 
and the cutoff the primary spectrum are near each other, the effect low 
down the atmosphere may depend some extent the angular distribution 
secondaries produced primary particles. The Stormer and shadow cone 
still opening for particles the east and north large zenith angles with 
energies well above the supposed spectral cutoff. 

That primaries energies low 1.5 Bev. may contribute neutron 
production sea level not surprising. The energy loss singly charged 
ionizing particle penetrating the atmosphere sea level about Bev. 
The absorption length about 140 gm./cm.? but the geometric cross-section 
for energetic nucleon interactions air about gm./cm.? The sea level 
nucleon probably generations from the primary particle and since 
the interactions involved are probably charge independent, half the ancestors 
are likely neutral. The limiting energy loss ionization, therefore, will 
much smaller than for the mesonic component cosmic rays. the other 
hand, each collision the secondary nucleon will usually have only fraction 
the energy the colliding particle. 

Meyer and Simpson (1955) show that the change primary spectrum 
between 1948 and 1954 extends quite high rigidities, probably high 
Bv. For such rigidities the energy vertically incident nucleons would 
quite sufficient contribute appreciably neutron production sea level. 
Measurements the latitude effect latitudes around the knee should 
repeated often possible sea level since they will give information 
changes the primary spectrum. 

The exact position the knee also depends longitude the curves 
Fig. show. The calculations for minimum vertical energy given above are 
based centered dipole without correction for longitude. The present mea- 
surements the east side the continent and those Neher, Peterson, and 
Stern (1953) and Neher and Stern (1955) were taken approximately the 
same course are comparable. If, however, the position the 
magnetic axis which controls the motion cosmic rays differs appreciably 
from that derived from surface magnetic measurements, these energies are 
wrong. more quantitative analysis the behavior around the knee best 
left until more complete knowledge becomes available. The same equipment 
used for these measurements has recently been taken the Antarctic aboard 
the U.S.N. Icebreaker Atka and much information will available from that 
expedition. 
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The longitude effect shown Fig. referring particularly the differ- 
ences between the two curves low latitude. For comparison with other 


measurements the value 40° has been calculated, where 


and are geomagnetic longitude and latitude. airplane altitudes Simpson 
(1951) found value 0.33 for geographic longitudes between 95° and 110° 
The present data gave value considerably lower, 0.1 (geographic longitude 
was used for rather than geomagnetic since the difference would not affect 
the present discussion). Close agreement with Simpson’s value should not 
expected sea level where the total latitude effect much smaller than 
airplane altitudes. 

The third prominent feature these results, namely the difference 
intensity above the knee the east and west sides the Canadian Arctic, 
very interesting. The effect shown Fig. the difference the east 
and west curves between geomagnetic latitudes 55° 80° and large enough 
significant. careful analysis all possible errors indicates that this 
difference real, but should pointed out that the lower curve from 
geomagnetic latitudes 56° 80° represents only the first five days the 
equipment was operating, and two these days there were some minor 
circuit troubles. These, however, resulted the loss data for only few hours. 
Making all possible corrections, the mean difference between the eastern and 
western curves, including all data between geomagnetic latitudes 55° and 80°, 
1.8%. Omitting the correction for world-wide changes based the monitor 
Ottawa, this difference 0.7% which still significant amount con- 
sideration the data used. 

This difference may due (1) the seasonal annual variation being 
different northern regions from that Ottawa (the eastern and western 
measurements this area were separated time over month), (2) some 
residual meteorological instrumental temperature effect, (3) distortion 
the earth’s magnetic field. Unfortunately none these can eliminated 
completely though the second, that meteorological instrumental 
temperature effect, seems very unlikely. regard (3), the true geomagnetic 
axis is, course, not the center the earth and the field not quite 


equivalent simple dipole. 


(2) The Meson Component 

Fig. shows the meson component measurements corrected only for baro- 
metric pressure. The correction used was that actually found the Arctic 
part the voyage. The same period was used the case the nucleon 
component. The value per cm. mercury 0.19+.01% per mb. 
The difference the two curves above the latitude knee assumed 
largely due differences temperature. The need for temperature correction 
is, therefore, very obvious. The surface temperatures between latitudes 60° 
and 70° were uniformly 11°C. higher the eastern part the voyage than 
the western. Since the intensity was also fairly uniform over these latitudes, 
correction coefficient could worked out from the data this interval. 
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COUNTS PER HOUR 256 


_Fic. meson component. Intensity plotted function geomagnetic latitude. 
The solid circles represent measurements the eastern side the North American Continent 
and crosses, the western side. 


Such correction applied the results throughout the voyage based sur- 
face temperatures gave impossible result and not reproduced. Another 
method was tried, namely using the Resolute data. Assuming that during 
the period the ship was Arctic waters, that above geographic latitudes 
about 55° N., the mean vertical distribution meterological conditions would 
roughly the same over the whole area the same over the ship and over 
Resolute, the meson intensity Resolute could then taken index 
temperature and annual variations. smooth curve was therefore drawn 
through the Resolute data (corrected, course, for pressure variations) and 
the ship results were corrected (after pressure correction) proportionately 
the variation the Resolute data from its mean. 

The difficulty with this correction was that purely arbitrary tapering off 
the correction had adopted the ship moved out the Arctic area 
into the Pacific, here the air mass was bound different. The results 
shown Fig. have the Arctic parts the curve corrected the above 
manner but corresponding correction for the parts the curve below 
geomagnetic latitudes about 50°. The only justification for this that 
gives result that looks reasonable. makes the measurements the eastern 
and western parts the voyage comparable, but the magnitude the inten- 
sity change from north south will unreal. This magnitude also unreal 
without any temperature correction Fig. Fig. and Fig. the 
case the nucleon data, 12-hr. averages were used but where the latitude was 
changing slowly that too many points had placed the curve, con- 
secutive groups were averaged. The smooth line Fig. was drawn through 
the points above the knee without reference the easterly and westerly 
differences since the spread the points was considerable. 


e 
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COUNTS PER HOUR 256 


GEOMAGNETIC LATITUDE 


The meson component. Fig. but corrected for intensity variations Resolute. 


Two features noted these results are: (1) They show the advantages 
working with the nucleon component where the complication temperature 
correction throughout the atmosphere very small exists all. (2) The 
position the latitude knee (best shown Fig. not very definite although 
considerably farther north than the most comprehensive work published 
previously this type measurement: Compton and Turner (1937) place 
the knee between 35° and 40° but show considerable increase intensity 
above it, the most northerly measurements shown being about 52° geo- 
magnetic latitude. The present data would place between 40° and 50°. The 
rigidity particles arriving vertically 45° 3.7 Bv. the minimum energy 
for protons 3.1 Bev. still possible that the variations the primary 
spectrum could affect measurements these latitudes. 

The longitude effect shown clearly lower latitudes and agreement 
direction with that found Millikan and Neher (1936). The value 


0.14. 


compare these with Millikan and curves, more satisfactory 
examine the quantity 5A, where and refer geographic latitude and 
longitude and the intensity above the knee. The present results give 
value 0.078% per degree longitude for this quantity 20°, while 
Millikan and curves give value 0.040%. The present results, 
therefore, indicate greater slope the variation intensity with longitude 
than shown Millikan and Neher. This would suggest that the maximum 
intensity the geomagnetic equator farther east than shown 
Millikan and Neher more agreement with the position that would 
expected from the eccentricity the dipole calculated from geomagnetic 


measurements. 
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IV. THE DIURNAL VARIATION THE NUCLEON COMPONENT 


Since changes the diurnal variation show some relation solar activity 
(Firor, Fonger, and Simpson 1954), differences the diurnal variation 
different latitudes might expected due the sun possible source 
least part the cosmic ray flux (Firor 1954). The results the Arctic part 
the cruise were examined for diurnal variation. The hourly values were used 
for this analysis and only the data from July September 22, that 
during the part the cruise well above the latitude knee. The diurnal values 
lower latitudes where the intensity was changing rapidly due latitude 
would impossible sort out. The barometer variation for the northern 
part was eliminated taking the differences from the hourly values and the 
value indicated the regression line for the pressure during that hour. 

The diurnal variation shown Fig. the interval being divided into two 
parts, the period July August and September September 22. 
The first part covers the period the ship was higher geomagnetic latitude 


OTTAWA JULY 1954 


OTTAWA AVERAGE 1952-1953 


UNIVERSAL TIME 


Fic. The diurnal variation the nucleon component. The first curve (at the top) 
represents the measurements above geomagnetic latitude 78°N. from July August 31, 
1954. The second curve represents measurements the western Arctic area from September 
September 22, 1954. The third curve shows the diurnal variation Ottawa from July 
September 22, 1954, and the fourth curve shows the diurnal variation Ottawa for nine 
months during 1952-53. 


than 78° During this period the cruise covered large range geomagnetic 
longitudes and for few days was within one degree the geomagnetic pole. 
Since the effective position the pole not certain, try convert geo- 
magnetic time seemed useless and, therefore, the results are plotted uni- 
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versal time. The usual curve for local time variation (see, for instance, Firor, 
Fonger, and Simpson 1954, Fonger 1953) shows variation the order 
amplitude with maximum just after local noon though the amplitude and 
shape the diurnal curve varies with solar activity. The data for the top curve 
Fig. were also converted local solar time but the smooth curve looks 
practically the same except for phase shift about six hours and that 
would roughly the difference between universal time and the local time for 
the mean position the ship. This curve does show diurnal variation but 
with almost the opposite phase that usually found during the past few years 
and with low amplitude little over and not very high statistical 
significance. The second curve shows the results from September 22. Here 
the local geomagnetic time variation due changing position the ship 
would spread over only about three hours and, therefore, the universal 
time plot would not appreciably different from that for local time except 
for appropriate phase shift. Actually this curve shows too wide scatter 
interpret any diurnal variation. 

The third curve represents the data for Ottawa during the whole period 
July September 22. The variation small and irregular. Since solar 
activity was low during that period, this consistent with other results. The 
fourth curve shows some previous data taken Ottawa with older neutron 
monitor much lower counting rate. The statistics were low that four-hour 
averages were used and the period covered nine months during 1952-53. 
The diurnal variation with peak little after local noon shown. Noon 
Ottawa hr. U.T. 


DISCUSSION AND SUMMARY 


The data presented here show the variation intensity cosmic rays from 
geomagnetic latitudes about 18° 89° N., both the nucleon compo- 
nent (measured local neutron production) and the meson component 
(triple coincidence counter telescope) sea level. The measurements were 
taken both sides the North American Continent and show the longitude 
effect very clearly. The series the eastern side represent narrow range 
longitudes and give the true latitude effect longitude approximately 
60° The knee the latitude effect about 52° the case the nucleon 
component and between 40° and 50° the case the meson component. 
The minimum energy protons allowed orbits reaching the earth vertically 
these latitudes 1.5 Bev. and about Bev. Since Meyer and Simpson 
(1955) have shown that the low energy primary spectrum this range has 
changed considerably since 1948, the exact position the knee and the shape 
the curve probably changes from time time, and such sea level measure- 
ments should repeated several times during solar cycle. The knee the 
meson component variation not sharp and the interpretation the results 
complicated difficulties temperature correction. The knee is, however, 
considerably farther north than was placed previous comprehensive 
set measurements Compton and Turner (1937). 
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The nucleon component shows ant difference intensity between 
measurements going north roughly along the geographic meridian containing 
the geomagnetic pole and going westward from there through the North West 
Passage and the Beaufort Sea. This differenc. interesting and the present 
data indicate high probability that real though definite explanation 
for presented. Residual meteorological effects, annual variation differ- 
ing different latitudes, distortions the earth’s magnetic field are sug- 
gested possible causes. 

The diurnal variation the nucleon component Arctic regions 
and irregular would expected during low solar activity. The diurnal 
variation the meson component was not analyzed since would compli- 
cated temperature effects. 

These data represent the northern set data taken with the same 
equipment going from the Arctic the Antarctic. The southern part was 
carried out collaboration with Simpson and with the assistance the 
U.S. Navy and will published separately. The complete measurements will 
give more accurate presentation the effect the earth’s field cosmic 
rays than has been available previously. 
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SOME NOTES THE THEORY THERMAL-NEUTRON 


ABSTRACT 


Equations for the asymptotic steady-state distribution neutrons homo- 
geneous and lattice-type reactors are derived without making any assumptions 
about the mechanism diffusion, except the obviously necessary one that the 
probability for neutron which born one given point captured 
second given point function only the distance between these two points. 
The equations are seen form that admits exponential solutions, these 
are written down, and equations for the Laplacians are derived. clear-cut 
definition the migration area lattice reactor given, and pointed out 
that reactor this type there unique value the Laplacian but rather 
range values. 


INTRODUCTION 


assume that the reader has some knowledge reactor theory (Glasstone 
and Edlund 1952; Goodman 1947; Soodak and Campbell 1950; Weinberg 
1952). The reactor theory discussed this paper restricted steady states 
thermal-neutron reactors either the homogeneous (part the lattice 
(part II) type, but should noted that some the ideas may applicable 
other types reactor. The necessity applying boundary conditions 
avoided assuming that the reactors extend infinity; the solutions the 
equations for neutron density etc. are then said asymptotic, and they can 
negative some parts space outside the actual reactor. any real re- 
actor the finite boundary introduces perturbation that appreciable over 
fraction the reactor, and this fraction increases the excess the multi- 
plication factor over increases. Very roughly speaking 1.3 the 
asymptotic solution will not approximate the true density appreciable 
fraction the reactor. 

The reader will recall that one can make several different assump- 
tions about the diffusion process reactor (elementary one-group diffusion 
theory, two-group theory, elementary diffusion theory thermal neutrons 
plus age theory slowing-down, etc.) and still arrive the same equation 
for the asymptotic thermal-neutron density 


(1) 


where the three-dimensional Laplacian operator, the position vector, 
and called the Laplacian the reactor. The reason all these assumptions 
lead the same equation that they have important component 
common: the probability that neutron which born one point will 
captured second point function only the distance between the two 
points. this paper take this necessary common component our only 
assumption about the diffusion process and derive equation equivalent 
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(1); substituting exponential solution our derived equation obtain 
expression for the Laplacian terms the moments the above prob- 
ability. Reference made the appropriate places the work Placzek and 
Volkoff (1947) and Plass (1952). 


EQUATIONS FOR NEUTRON DISTRIBUTION AND LAPLACIAN 
HOMOGENEOUS REACTOR 


Let p(r) number thermal neutrons being captured per cc. per sec. 
(the neutron density proportional p), 


multiplication factor the reactor 


number thermal neutrons produced per thermal neutron 
captured, 


probability that neutron which produced fission 
caused the capture thermal neutron and which 
escapes resonance capture suffers thermal capture the element 
volume dr’ r’. 


The desired equation for p(r) obtained equating p(r) the sum all 
the neutrons being created per second throughout the reactor and being cap- 
tured per cc. thermal neutrons 


can once write down exponential solution (2): 
where must satisfy 
(Note that (3) satisfies (1).) 

order get explicit expression for proceed follows. Choose 


the x-axis parallel Then K-r Ax. Expand the exponential powers 
and change the order integration and summation; the result 


j=0 


where 


times the (27)th moment P(r) the x-direction. called the 
migration area the neutrons and usually denoted Using the expan- 
sion the reciprocal the series (5) get 


(6a) 


@w 
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Finally invert this series, i.e. get series powers instead 


The above treatment would seem the most direct way obtaining 
equation (7) and showing that one gets solutions the form (3) indepen- 
dently the diffusion mechanism that assumed. One could, however, 
obtain the same results simple extension the work Placzek and Vol- 
koff (1947); they showed the Fourier-transform method that infinite 
uniform multiplying medium with the rate creation fission neutrons 
falls off exp great distances |z| from plane source neutrons 
independently the diffusion mechanism, and they also obtained the equiva- 
lent the first two terms (7). 


EQUATIONS FOR NEUTRON DISTRIBUTION AND LAPLACIAN 
LATTICE-TYPE REACTOR 
exact theory the distribution neutron density more complicated 

for this type than for homogeneous reactor. The method that usually used 
get approximate equation for this distribution that the equivalent 
homogeneous reactor. Plass (1952) has treated the problem lattice-type 
reactor using the elementary equation for microscopic diffusion and solving 
approximately the cellular method. this paper treat the problem 
essentially the same method have used for the homogeneous reactor. 
addition spatial periodicity make some other assumptions order 
get simple derivations the equations for neutron distribution and La- 
placian; one could, however, drop these assumptions without changing the 
essential form the derived equations, but the theory would then con- 
siderably more complicated. accordingly assume that 

(i) the reactor consists lattice parallel and exactly similar rods 
fissile material; 

(ii) the medium between the rods contains fissile material and has the 
same periodicity the lattice rods; 

(iii) the cross-section the reactor, the neutron density, etc., are indepen- 
dent distance parallel the rods (this makes the problem two-dimensional) 
(iv) the thermal-neutron absorption rod can considered taking 
place its axis. 

The notation essentially the same for the homogeneous reactor but 

adapted follows our lattice reactor: 

the number thermal neutrons one generation captured 
fissile material per thermal neutron the preceding genera- 
tion captured fissile material, 

the position vector plane normal the rod axes, 

the number thermal neutrons captured per second per unit 
length the rod whose axis position 
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the probability that neutron which has escaped resonance 
capture and also thermal capture the medium between rods 
suffers thermal capture the rod whose position with 
respect the rod which the parent thermal neutron was 
captured (note that P(s) 1), 

1 23 

P(s). 

get the counterpart equation (2) simply equate the total 

number thermal neutrons that originate per second throughout the reactor 

and are captured 


(8) M2; = 


which admits solutions the form (3), provided that 


(10) the counterpart (4), and (7) follows from from (4) when the 
moments are defined (8). should noted that the values the 
will depend the direction chosen for the x-axis. This means that the lattice- 
type reactor has unique Laplacian but rather range Laplacians. Usually, 
however, this range small that for practical purposes can considered 
zero. 

get expression for the rate capture thermal neutrons any 
point p(r), let the value p(r) due unit value p,(r’). Then 


where Q(r) has the periodicity the lattice. Plass (1952) quotes formula 
(his number (10)) which equivalent (11). His justification would appear 
based the assumption the validity elementary diffusion theory 
both moderator and fissile material, whereas our derivation independent 
this assumption. 
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RADIATION RESISTANCE DIPOLES INTERFACE 
BETWEEN TWO DIELECTRICS! 


ABSTRACT 


Exact expressions are derived for the radiation resistance electric and 
magnetic dipoles located the plane interface between two semi-infinite 
lossless dielectrics. 


The radiation resistances both electric and magnetic dipoles situated 
lossless media infinite extent are well-known quantities (Schelkunoff 1943). 
For example, the length the current element electric dipole the 


where the wavelength the radiation. The corresponding value for 
small loop area magnetic dipole given 


(2) 320 ohms. 


When the dipoles are situated inhomogeneous medium the values 
and can somewhat different. Sommerfeld and Renner (1942) have 
studied the case where electric dipoles are situated medium dielectric 
constant height over flat interface, below which there medium 
dielectric constant Their solution for the radiation resistance for 
electric dipole oriented parallel the interface given 


and 


Sommerfeld and Renner were mainly interested the case where was com- 
plex with large imaginary part corresponding highly dissipative medium. 
They have evaluated the integral approximate manner with the restric- 
tion that was much greater than the real quantity and the height 
not too small. 

possible derive exact expression for the ratio for the interesting 
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case when the dipole situated within the interface, with the 
restriction that and are real. That is, both media are perfect dielectrics. 
The ratio now given by: 
(5) imaginary part 


The integral can broken into two symmetrical parts that 
with 
k,* 2 2 2,2 


where identical form but with interchanged with These 
integrals can evaluated standard means, give 


The ratio plotted Fig. function the relative refractive 
index defined 
N = ko/ky = (€2/€1)*/?. 


500 


RELATIVE 
RADIATION 
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Relative Refractive 


Fic. The ratio the radiation resistance electric and magnetic dipoles interface 
the corresponding value free space. 


100 
& 
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values greater than about the following expression valid within 
one per cent: 


The corresponding result for magnetic dipole small wire loop can also 
derived similar manner. For magnetic dipole oriented normally 
the interface was shown previously (Wait 1952) that the radiation resistance 
given 


where 


When the height the dipole allowed approach zero easily follows 
that 


( ) Ro ki (k2—k}) (ki p) pdp (ke p) pdp 
which evaluated give 
This ratio for the vertically oriented magnetic dipole small loop situated 
the interface also plotted Fig. function the relative refractive 
index For large values greater than about the corresponding approximate 
relation 
probable that these results can generalized radiating systems 


finite size, lying dielectric medium (i.e., frozen ground ice), long 
their dimensions are small compared the wavelength. 
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TRANSIENT FIELDS VERTICAL OVER 
HOMOGENEOUS CURVED 


ABSTRACT 


Expressions are derived for the transient fields short vertical antenna, 
situated smooth spherical conducting earth, and energized current 
which discontinuous time. When the antenna current linear function 
time, the radiation field flat perfectly conducting earth step-function 
form. The departure from the step shape the field shown due the 
finite conductivity and dielectric constant the ground, the induction and 
static fields the antenna, and the curvature the earth. 


INTRODUCTION 

Some interest has been shown recently the fields electric dipole 
excited transient current. Van der Pol (1955) has recently presented 
solution for the dipole, whose moment proportional step function 
time, located the plane interface two different dielectrics. similar prob- 
lem was presented concurrently Poritsky (1955) who employed the artifice 
representing the radiation from the dipole superposition plane waves. 
The case where the transient dipole located over imperfect dielectric 
dissipative medium even greater interest. approach this problem 
was published recently Wait and Froese (1955) who considered the source 
plane wave incident obliquely the plane interface semi-infinite 
homogeneous dissipative medium. The electric vector which was taken the 
plane incidence was proportional step function time. Sometime earlier 
the author (1953) derived exact solution for magnetic dipole small loop 
situated infinite conducting medium and excited step-function cur- 
rent. The case (Wait 1951) the transient magnetic dipole over stratified 
ground has also been treated using approximate method which valid 
when the observer near the source. The latter problem has recently been 
extended Belluigi (1954). 

the purpose this paper present derivation for dipole short 
vertical antenna located just above spherical dissipative smooth earth. The 
solution facilitated employing new type series expansion due 
Bremmer (1949) for the fields dipole whose current varies har- 
monically The first term this expansion corresponds the field 
dipole over flat earth given Sommerfeld (1926) and Norton (1937). 
The succeeding terms are proportional ascending powers (1/a) where 
the radius the earth. employing techniques equivalent Laplace 
Fourier analysis, the transient fields the dipole are found using 
suitable frequency spectrum for the antenna current 


FLAT EARTH CASE 


starting point the earth considered flat and homogeneous with 
conductivity and relative dielectric constant For short vertical antenna 
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height excited current the vertical electric field 
distance the ground, given (Bremmer 1949; Norton 1937; 
Sommerfeld 1926) 


where 
tp 
with 


not exact (Norton 1937), being approximation valid for 
2 


The antenna current now expressed function time f(t) 


other words I(s) the Laplace transform where replaces tw. Written 
real part insure convergence for various forms f(t). The corresponding 
transient electric field then related E(iw) 


E(s) Since f(t) and are specified, E(s) can obtained from 
equation (1). The problem then find from E(s), being equivalent 
the evaluation the Bromich contour integral 


1 c+ too 
which written symbolically 
(7) e(t) 


For convenience, the function replaced that current the 
antenna linearly increasing function time being equal amp. 
to. Then 

(8) I(s) = Lele t/to = Ip/s? to. 


then follows without difficulty that 
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where t—D/c and 


For perfectly conducting ground the factor unity and consequently 
F(s)/s the above equation unity. It, therefore, concluded that for 
large distances small times, such that D/c, the field characterized 
step-function form. the antenna current was step function 
instead linearly rising function, the response would 


where 


with being the unit impulse function t—D/c With similar 
reasoning follows that the current impulsive the response 


Todt; 


with 


The core the problem then evaluate the inverse transform F(s)/s 
obtain the response function approach expand F(s)/s 
asymptotic expansion p(s) follows: 


Furthermore, displacement currents the ground are neglected such that 
p(s) 


where then follows that 


The inversion this transform readily carried out give 


m=1 


which course identical 


This equation can checked directly taking the Laplace transform the 
right-hand side leading 
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which identical the function F(s)/s given equation (2) with the sub- 


and 
These integrations can carried out yield 
(22) (D/c) 1 
and 
with 


now proposed obtain correction for the displacement currents 
the ground. This effected expanding 1/p(s) series ascending powers 


(24) 


doubtful the coefficients powers higher than one are meaningful 
since the stated error the defining equation, (3), for p(s) the order 
the term Inserting this expression into the series formula for F(s)/s 
and retaining only the first correction term leads 


where 
m—1 2 
(26) 
can verified that 
The correction term proportional can neglected when 
BK, 


is, therefore, reasonably justified neglect displacement currents the 


WAIT: TRANSIENT FIELDS 


CORRECTION FOR CURVATURE 
The expansion formula due Bremmer (1949) for the field vertical 
electric dipole spherical earth given equation (1) F(tw) replaced 
terms A‘, etc., 


where 


where the radius the earth. Replacing and expanding F(iw) 
inverse powers follows that 


F(s F(s 


readily carried out. Using this result, the final expression for the response 


2/2 


terms X3, etc. 


anc 


The above formula neglects displacement currents throughout. first order 

correction for these was given equation (25). That is, the quantity 

added the right-hand side equation (31). 
Formally the response functions and for step function and im- 


*Equation (28) its present form was communicated the author privately Dr. 
Bremmer. believed more accurate than equation Bremmer (1949). 
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mho/m 


Response B(t') 


(microseconds) 


The transient response distances 20, 50, and 100 km. for longer times where the 
static and induction fields are significant, for electric dipole energized linearly rising 
current. 


pulse function sources, respectively, can obtained from carrying 
out the operations indicated equations (12) and (14). 

Some numerical calculations were carried out for the function and 
they are illustrated Figs. Fig. the response calculated from 
and km. and and mhos per meter. The solid curves 
correspond when displacement currents are considered negligible. 
The dotted curves show the first order correction for displacement currents 
assuming 10. The effects the curvature corrections the flat earth 
formula are illustrated Fig. The first curvature correction refers equa- 
tion (31) with the term proportional being omitted. The second curva- 
ture correction includes terms proportional both and 

interesting note that the flat earth formula differs only appreciably 
from the curvature corrected formula small times such that 0.1 usec. 
distances greater than km., however, the curvature correction becomes 
significant larger times. Unfortunately the series type solution here 
becomes very poorly convergent and terms proportional etc. would 
needed for greater distances smaller times. 

times greater than the induction field becomes significant. This 
indicated Figs. and the slight rise B(t’) above unity for the 
range For greater times, both the induction and static field 
become quite significant illustrated Fig. for 20, 50, and 100 km. 
this region the time domain, the field almost independent ground 
conductivity and the curvature corrections are negligible these distances. 
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EXPONENTIAL EXCITATION 

rather striking that the transient responses for dipole over flat 

highly conducting ground can expressed quite simple form. Another case 


that yields answer closed form when the dipole energized 
exponential current (Jo Note that 


The response distance along the surface flat earth, neglecting displace- 
ment currents, then given 


where 

fod F(s)s —sD/e 


and where F(s) given (2). 
This can conveniently rewritten 


Therefore 


where 


F(s) 


Using the convolution theorem follows that 


(37) 


and therefore 


, , t’ /2a—al’ 
ot 4a 
When these values for and are inserted into equation (36), the 
response current u(t) completely specified. The above 
exponential source can easily extended suddenly applied sinusoidal 
source. For example, the current sin where real, then the re- 
sponse the imaginary part Similarly the current were 
cos u(t) the response would the real part e(—7Q, interesting 
note that when the response the current, simply the 
steady state value characteristic vertical dipole flat earth. Since 


(40) 
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reached the steady state condition. quite possible, using tables the 
probability integral with complex argument, evaluate the way which the 
response builds its steady state value. For smaller 
times, however, the form the curves would affected curvature and the 
displacement currents, both which have not been considered this section. 


CONCLUDING REMARKS 


The behavior transient dipole over smooth spherical conductor has 
been shown mainly determined the finite conductivity the ground. 
For distances greater than km. and times less than 0.1 following the 
first arrival the signal the curvature the earth begins have marked 
effect. distances between and 100 km. the induction field the source 
dipole becomes appreciable for times greater than usec. 

The extension the results this paper distances out 2000 km. would 
require the use different approach. this case would probably desir- 
able start with the complete residue series (Bremmer 1949) for the field 
and then apply suitable Fourier transformation. 
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SMALL-ANGLE SCATTERING FAST POLARIZED NEUTRONS 
HEAVY 


SAMPLE 


ABSTRACT 


Detailed calculations have been carried out which indicate that the small- 
angle scattering fast neutrons lead depends the polarization, spin 
orientation, the neutrons. When the scattering neutrons whose spin vectors 
point upward observed the horizontal plane, more neutrons should found 
scattered the right than the left. For completely polarized 3.1 Mev. neutrons, 
the theory predicts maximum “right intensity ratio 14.5:1 
scattering angle 0.5°, the ratio decreasing 5°, and approaching 
unity rapidly the scattering angle increases. 


INTRODUCTION 


date, most successful methods measuring the percentage polarization 
fast neutrons (Adair, Darden, and Fields 1954; Huber and Baumgartner 
1953; Ricamo 1953 have involved resonant elastic scattering, with the 
attendant uncertainties due the variation with energy the polarization 
dependence the scattering cross-section. Adair al. (1954) have discovered 
non-resonant polarization dependence the 90° scattering neutrons 
nuclei the region 110. This effect, most 20% difference intensity 
between left- and right-scattered neutrons for 50% polarized 0.4 Mev. neu- 
trons, rather weak for use neutron polarimeter. Schwinger (1948) 
showed, the Born approximation method, that small-angle scattering 
heavy nuclei should strongly polarization dependent because the inter- 
action the Coulomb field the nucleus with the magnetic moment the 
moving neutron. Calculation the perturbation sphere’’ wave 
function this magnetic interaction yields expression for the differential 
cross-section differing from Schwinger’s only that the polarization depen- 
dent term about 25% greater. 


CALCULATION THE SCATTERING CROSS-SECTION 


The Schrédinger equation describing the motion the neutron may 
written 


(1) 
where 2mE/h?, 


S=spin 
neutron mass, 
total neutron energy, 


—1.9135 nuclear magnetons, the neutron magnetic moment, 
charge target nucleus. 
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SAMPLE: SMALL-ANGLE SCATTERING 


The expression representing the magnetic interaction, may derived 
considering the force exerted magnetic dipole the apparent mag- 
netic field arising from its motion the electrostatic field nucleus. (The 
magnetostatic fields nuclei and electrons may safely neglected.) 

Experiment (Whitehead and Snowdon 1953) shows that the differential 
elastic scattering cross-section lead for neutrons energy about Mev. 
resembles closely that perfectly elastic sphere scattering angles less than 
60°. The wave function should then good approximation 
from which begin perturbation calculations. 

the boundary conditions 


where the expression describes the spin orientation 
the neutron, and being the usual spin eigenfunctions, and a(}) and 


(2) 


where kr, 


are the spherical harmonics defined Condon and Shortley 

(1951). 
Equation (1) becomes 


where has been neglected with respect the right. 

obtain the radial equation corresponding (4), convenient ex- 


where 


Jil 
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and 3,0, are the Clebsch-Gordan coefficients defined 
Condon and Shortley (1951); 


i+} J 
substitution equations (5) and (6) (4) gives 
with the boundary conditions 
forp 
The particular integral equation (7) may found terms and 


constructing Green’s function satisfying the boundary conditions. 
The asymptotic expression for 


where S,(d) = G? K,-H;?, T G, 


2 

t 2 


Evaluation the necessary integrals terms tabulated functions carried 
out Appendix 

When (8) substituted equation (6), the being expressed terms 
spherical harmonics, and the result added the difference between equation 
(3) and the incident wave, the resulting asymptotic 
expression for the scattered wave 

- tkr 


For convenience numerical calculation, the spherical harmonics have been 
expressed terms associated Legendre polynomials, and tabulated 
Jahnke and Emde (1945). 


| 


SAMPLE: SMALL-ANGLE SCATTERING 


The polarization vector, the incident beam defined 


where the Pauli spin vector, and *,. indicate respectively complex con- 
jugate and Hermitean adjoint. Introducing the unit vector normal the 
plane the reaction, defined 


the differential cross-section obtained from equation (9) may written 
1 2 
+;(P-n)[ReA 


where and indicate respectively the real and imaginary parts 
the quantity. Evaluation the terms and outlined Appendix 
DISCUSSION 

can seen from equation (10) that polarization has maximum effect 
the cross-section when the direction observation chosen that 
normal the plane the reaction. this case the interference term, the 
third term (10), has its maximum absolute value. This term changes sign 


POLAR ANGLE DEGREES 


Fic. Differential scattering cross-section lead for polarized 3.1 Mev. neutrons. 
C—Schwinger, P.n=1 F—Hard-sphere scattering 
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accordingly the scattering the right left sense, that is, changed 
The differential scattering cross-section for 100% polarized (|P| 
3.1 Mev. neutrons scattered lead plotted function scattering angle 
Fig. The reaction appreciably polarization sensitive only for less 
than 10°. For comparison, Schwinger’s expression for the cross-section 
plotted Fig. with inserted for the nuclear scattering amplitude 
(Schwinger 1948, equation 10). The only difference between the two expressions 
that the polarization dependent term equation (10) 25% larger. 

attempt Longley (1952) measure the polarization neutrons 
from D(d, scattering led inconclusive results. attempt 
somewhat different method this laboratory (Sample, Neilson, and 
Warren 1955) failed because the difficulty collimating fast neutrons 
sufficiently well investigate scattering very small angles. 
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APPENDIX 
Evaluation 
can shown (Watson 1945) that 


where and are any cylinder functions. Expressing the integrands 
and terms Bessel and Neumann functions (equation 8), and putting 


readily expressed terms Bessel and Neumann functions, and 


where being the lowest value for which these approxi- 
mations have the desired degree accuracy. 


SAMPLE: SMALL-ANGLE SCATTERING 
APPENDIX 
Calculation the Terms Equation (10) 
Since converges rapidly, calculation numerical values straight- 
forward. the other hand, converges slowly, that necessary 


obtain approximate closed form for the sum high order terms. From 
equation (13) 


may written one the forms the Laplace integral: 


The second sum equation (14) becomes 


0 lo+1 l= l 
where cos sin cos The sums are 


@ 


l 
Substituting these values equation (15), and rewriting the finite sums 
terms associated Legendre polynomials, gives 


(cos 


(16) 


(The only term with vanishes.) 
The integral equation (16) has the value cot Substituting this (16) 
and the result equation (14), 


calculated applying equation 2.3 the review article Blatt and Bieden- 
harn (1952). 
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MAGNETIC SUSCEPTIBILITY LOW RESISTIVITY 
n-TYPE GERMANIUM! 


ABSTRACT 


Some results the magnetic susceptibility low resistivity poly- 
crystalline germanium observed between room and liquid nitrogen tempera- 
tures are presented. The contribution the susceptibility from both the free 
and bound charge carriers has been derived from the experimental results and 
appears vary inversely with temperature. effective mass for the charge 
carriers m-type germanium found 0.16 the free electron mass. 


INTRODUCTION 


The work the present time the magnetic properties germanium 
has been chiefly concerned with the determination the susceptibility pure 
germanium and only the work Stevens and Crawford (1953) and van 
Itterbeek (1955) are any results given for impure germanium with known 
amounts impurity. For this reason the following measurements were made 
more complete set polycrystalline specimens which include those 
having impurity concentrations that would class the semiconductor essen- 
tially The magnetic susceptibility presumably impure ger- 
manium between 290°K. and 80°K. has also been measured Owen (1912) 
and Squire (1940). and Klemm (1940) report measurements ger- 
alloys but the concentrations antimony are too great 
for the results interpreted terms semiconductor theory. Busch and 
Helfer (1954) have recently reported the magnetic susceptibility pure 
germanium between 1200°K. and 40°K. while van Itterbeek (1955) give 
measurements some samples between room and liquid hydrogen tempera- 
tures. 

explain (1912) results the temperature dependence the 
susceptibility germanium Aigrain, Dugas, and Etzel (1951) calculated the 
amount impurity necessary, the assumption that all the charge 
carriers are bound the temperature liquid nitrogen and are free room 
temperature. These calculations indicate impurity concentrations solid 
solution approximately 1%. This greatly exceeds the solid solubility limits 
commonly found for impurities germanium: for instance Thurmond and 
Struthers (1953) have found the maximum solubility antimony and copper 
X-ray measurements, made the author, the lattice parameters some 
samples, are agreement with Thurmond and Stru- 
thers’ results and indicate that solic solubility antimony germanium 
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high extremely unlikely. fact only after prolonged annealing 
900°C. and subsequent quenching was there any lattice parameter change found 
the alloys nominally containing antimony. The results indicate small 
increase the lattice parameter germanium which would expected 
antimony forms substitutional solid solution germanium low solid 
solubility. 


SPECIMEN PREPARATION AND METHOD MEASUREMENT 


The specimen rods were prepared for Johnson Matthey and Mallory 
Ltd. melting germanium together with various amounts antimony 
weight. The alloys were melted them graphite reducing 
atmosphere, cast ingot, annealed and recast rods in. diameter and 
rods supplied showed the presence excess antimony and other precipitated 
impurity, probably germanium oxide. order put the greatest possible 
amount antimony into solid solution, the rods were subjected long periods 
annealing vacuo (up two weeks 860°C.) and then quenched. Micro- 
scopic examination then showed that some second phase existed the material 
but measurements electrical resistivity indicated increase the amount 
dissolved impurity. The amount antimony actually solid solution was 
deduced from the resistivity values using the results Pearson, Struthers, 
and Theuerer (1950). 

The susceptibility the specimens was measured the Gouy method, 
using specially developed servo-balance apparatus.* Current for the 
was supplied from d-c. generator and stability the magnetic field 
ensured automatic regulation the current the field coils the genera- 
tor.* 

estimated that the measurements the susceptibility have absolute 
accuracy and relative accuracy 0.5%. Any specimens that showed 
ferrous impurity indicated the dependence magnetic susceptibility 
magnetic field strength were rejected. 


DISCUSSION RESULTS 


The susceptibility various specimens function temperature 
shown Fig. The lowest curve refers pure zone-refined 
and the upper curves samples with successively greater amounts impurity 
for which the material becomes increasingly diamagnetic. None these 
specimens approaches Owen’s value e.m.u./gm. 100°K. 
that his material must have been extremely impure. The curve for pure ger- 
manium agrees well with that found previous workers (Stevens and Craw- 
ford 1953; van Itterbeek al. 1955; Busch and Helfer 1954) and the most 
impure samples Stevens and Crawford agree roughly with specimen 
though with rather smaller negative slope. 

The analysis the susceptibility results will follow closely that originally 


*To described separate publications. 
indebted Dr. Rutter the University Toronto, Department Metallurgy, 


for supplying this sample. 


= 
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perature. 


given Busch and Mooser (1953) for their measurements gray tin. xz, 
the total susceptibility semiconductor, may regarded made 
four parts: 


due electrons the valence band and inner shells. 

due free electrons and mobile holes.* 

due electrons and holes bound impurity sites. 

precipitated impurity, lattice dislocations, vacancies, inter- 
stitial atoms, etc. 


may remarked that this method dividing the total susceptibility 
semiconductor into the above components somewhat arbitrary and the 
division and based the simple model electron (or hole) 
being bound the coulomb field impurity site until sufficient thermal 
energy acquired ionize the trapped electron when then assumed free 
move through the crystal lattice. The term due the (neutral) ger- 
manium atoms themselves and can determined experimentally measuring 
the susceptibility high resistivity germanium below the temperature the 
intrinsic region (i.e. approximately 650°K.). One might expect that would 


*In the present n-type specimens the concentration holes negligible that there will 
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temperature independent; however, previous measurements high purity 
samples (Stevens and Crawford 1953; van Itterbeek 1955; Busch and 
Helfer 1954) and the present results for the zone-refined specimen and speci- 
men indicate that varies significantly with temperature. present 
explanation for this anomalous behavior has been given. 

Since the amount precipitated impurity small the temperature depen- 
dence due precipitated impurity should negligible and the tem- 
perature range under consideration unlikely that there will any change 
the number lattice imperfections. shall therefore assume that the 


Now 


and have said, determined the data pure germanium; 
subtracting the curve for the pure substance have: 


which the first two terms are the susceptibility the free and bound charge 
carriers and temperature independent. Now assume quasi free 
electron model for the electrons the conduction band, the susceptibility 
the charge carriers can calculated theoretically the analysis Stoner 
(1935) extended Wilson (1953) take account electron-lattice inter- 
action inclusion the effective mass the electron. non-degenerate 
electron gas the contribution the mass susceptibility given 


(1) 
where the number free charge-carriers, 
the Bohr magneton, 
and are the free and effective electron masses, 
the density, 
Boltzmann’s constant, 
the absolute temperature. 
For degenerate electron gas occurring samples with high impurity concen- 


trations normal temperatures moderately impure samples low tem- 
peratures, statistics must used and Eq. (1) modified 


4m*y? 
where Planck’s constant and the other symbols have the same meaning 
Eq. (1). 


The expression for the susceptibility contribution the bound electrons, 
based hydrogen-like model, has the form (Busch and Mooser 1953) 


(2) 
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the first term being the paramagnetic spin contribution and the second the 
diamagnetic orbital component, 
where the number bound electrons, 
the electron charge, 
the dielectric constant the material, 
the first Bohr radius for the hydrogen atom, 
the velocity light, 


and the other symbols have the same meaning Eq. (1). 

The effective mass electrons the conduction band unique only the 
Fermi surface sphere. Cyclotron resonance experiments (Lax al. 1954) 
have established the need for ellipsoidal energy surfaces represent the con- 
duction band germanium. Enz (1955) has shown that using the appropri- 
ate averaging process the value used (1) for ellipsoidal 
energy surfaces given 


where and are the longitudinal and transverse mass ratios. Using the 
values 0.083 and 1.4 for and determined Lax al. yields 
crystal such germanium can give information the individual 
mass ratios, and 

The experimental results were originally interpreted the assumption that 
and varied with temperature (cf. Putley 1949). With this approach 
had assumed that decreased with temperature. similar assumption 
was found necessary Debye and Conwell (1954) fit theoretical curves 
the data the electrical properties However the work Hung 
and Gleissman (1950) the Hall effect germanium suggests that 
high impurity concentrations bound charge carriers exist down liquid 
nitrogen temperatures that and are independent temperature. 

then assume that and are constant for the high impurity 
samples and and that the free charge carriers behave non- 
degenerate* electron gas, see from equation (1) that vary inversely 
with temperature. can seen Fig. xreq does indeed vary linearly with 
1/T expected the present interpretation the results. Using the slopes 
the lines Fig. for specimens and effective mass may then 
calculated using equation (1). The value found which agrees 
well with the value m/6 deduced Stevens and Crawford (1953) and also 
reasonably well with the value 0.14 calculated from the previously men- 
tioned data Lax al. (1954). 

For sample with electron density determined 
from resistivity measurement, the assumption constant leads 
value 0.14+.01 for m*. This apparent change the effective mass can 
interpreted one two ways; namely, that the effective mass function 


calculation the degeneracy temperature for these impurity concentrations indicates 
between 60° and 100°K. depending the effective mass chosen. 
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Fic. susceptibility the free and bound charge carriers and the precipitated impur- 
ity the alloys function 1/T. 


impurity concentration indicated the experimental evidence Crawford 
al. (1955) that lower temperatures some the electrons are bound 
impurity sites specimen. determination the Hall coefficient 
function temperature for sample indicates that the number charge 
carriers indeed not independent temperature. appears therefore that, 
these moderately high impurity concentrations, some electrons are bound 
that now varies with temperature. For sample the Hall coefficient 
showed change approximately 20% between liquid nitrogen and room 
temperature. Assuming this change the Hall coefficient measure the 
change which probably true for sample with the impurity concentra- 
tion calculation the temperature-dependent part can made 
from equations (1) and (2) using the found for the less pure samples. The 
dashed line Fig. represents the calculated susceptibility the free electron 
gas were constant, while the dotted line includes the calculated suscep- 
tibility the bound electrons, the numbers being based the Hall 
coefficient measurements. can seen Fig. extremely good agreement 
with the experimental curve found when included the expression for 
should perhaps mentioned that the effective mass used both 
expressions for the susceptibility the free and bound charge carriers the 
same, conformity with the simple model which equations (1) and (2) 
are based (cf. Debye and Conwell 1954). 

intended extend these measurements down the temperatures 
liquid helium and should possible, from the results, obtain further 
useful information about the behavior free and bound charge carriers 
impure semiconductors. 
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METEOR ECHO DURATIONS AND VISUAL 


ABSTRACT 


The statistical relation between the radio echo duration and the visual mag- 
nitude has been investigated for approximately 3300 meteors observed the 
combined Dominion National Research Council programat Ottawa. 
Both echo durations and visual magnitudes were reduced absolute values, 
defined those for meteor the zenith height 100 km. For meteors 
the absolute magnitude range straight line relation exists between 
log absolute duration and absolute magnitude, longer durations corresponding 
brighter meteors. Applying current meteor theory these data indicates that 
meteor absolute magnitude produces 2X10" electrons per cm. path 
length. For Perseid meteors the ratio between ionization produced and visual 
luminosity almost constant the range absolute magnitude. For 
other meteors this ratio seems decrease somewhat for the brighter objects. 
These results lead higher electron densities for the bright meteors than had 
previously been estimated. 


Since the summer 1947 the Dominion Observatory and the National 
Research Council Ottawa have cooperated combined program meteor 
observation. the course this work large number meteor echoes has 
been recorded the Ottawa station frequency 32.7 Mc./sec. (McKinley 
and Millman 1949). Many the more prominent echoes, corresponding the 
overdense ionization trails with greater than 10" electrons per cm. path length 
(Greenhow 1952; Kaiser and Closs 1952), were produced meteors for which 
good visual observations were available. 

was early recognized that there was general statistical relation between 
the observed durations these echoes and the brightness the meteors, 
the longer durations corresponding the brighter meteors. Preliminary 
relations, based 700 echoes, were published for the Perseid and Geminid 
meteors Millman (1950). These values were found sorting the meteors 
into classes apparent visual magnitude and then averaging the logarithm 
the duration seconds for each class. These preliminary curves have been 
used number investigators discussing the theory meteor echo 
formation (Greenhow and Hawkins 1952; Kaiser 1953; Lindblad 1952). 

The 1950 curves were published illustrate qualitatively general relation 
and were not meant applied quantitatively other problems. They are 
open criticism the following counts: only apparent visual magnitudes 
the meteors were used and these were not reduced standard absolute 
system; (b) the durations the echoes were not reduced standard zenithal 
distance; (c) taking means log durations artificially decreases the slope 
the curve owing the large statistical spread about the mean and the relatively 
small numbers very long and very short duration echoes; (d) the number 
echoes used was rather small considering the statistical spread observed 
values. 
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MILLMAN AND METEOR ECHO DURATIONS 


The present paper attempt correct these faults. The observational 
material used consists 3283 meteors observed Ottawa during the 
1948, 1949, and 1950. For each meteor recorded radio echo and visual 
magnitude estimate were available. The echoes were recorded Ottawa, lat. 
45° 20’ 32” N., long. 75° W., elevation 320 ft. above sea level, fre- 
quency 32.7 Mc./sec., transmitter peak power 400 kw., pulse width usec., 
receiver sensitivity w., antennae—half-wave horizontal dipoles ori- 
ented and mounted quarter-wave above ground, display—range-time 
sweep. The visual magnitudes were recorded team visual observers, 
six eight number, estimating the nearest half magnitude comparing 
the meteor with standard stars. Where two more independent magnitude 
estimates were available these were averaged. The meteors were 
with the various better-known showers consideration the date 
observation, and the position, orientation, and general appearance the visual 
trail. many cases confirming evidence for the identification appeared the 
radio echo. The distribution numbers among the various showers was 
follows, the generally accepted mean observed velocity km./sec. being 
listed brackets for each shower: 


Number 
Shower meteors 

Lyrids (49) 
Aquarids (64) 
Aquarids (40) 158 
Perseids 1404 
Orionids (66) 
Leonids (72) 
Geminids (35) 256 
Non-shower 1420 

Total 3283 


evident from the above table that Perseids provide the best statistical 
material from the standpoint numbers and homogeneity. 

The absolute magnitude meteor defined the magnitude would 
have were located the observer’s zenith height 100 km. 
convert absolute magnitude the apparent visual magnitude must cor- 
rected for the difference between the actual distance the meteor and 100 km., 
and for the increased atmospheric absorption away from the zenith. con- 
venient table for the first correction can computed assuming that the 
brightness point source falls off inversely the square the distance. 
This reproduced Table applying this table the mean ranges the 
enduring parts the echoes and characteristics) were used. 

The absorption correction was taken from the values given the Handbuch 
der Astrophysik (Lundmark 1932), assuming the mean color index meteors 
—1.8. This correction used listed Table II. may noted here 
that, although Table has been extended elevation for the sake 
completeness, practically meteors recorded the radio equipment had 
elevations under 10° and very few were low 15° elevation since the 
maximum range recorded varied from 270 390 km. Mean elevation the 
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TABLE 
DISTANCE CORRECTION FOR METEOR MAGNITUDE 


Distance Correction Distance Correction Distance Correction 

(km.) (mag.) (km.) (mag.) (km.) (mag.) 

123 195 
+0.5 

128 204 
+0.4 

134 214 
+0.3 

141 224 
+0.2 

148 235 
+0.1 —0.9 

155 245 
0.0 

102 162 256 

107 170 269 

112 177 282 

117 186 296 

123 195 310 

TABLE 


EXTINCTION CORRECTION FOR METEOR MAGNITUDES 


Mean Mean Mean 
elevation Correction elevation Correction elevation Correction 
(h) (mag.) (h) (mag.) (h) (mag.) 


meteor trail was taken from the visual plot. All visual apparent magnitudes 
were reduced visual absolute magnitudes through the use Tables and 

The absolute duration meteor echo was defined the duration that 
would observed with the equipment employed the meteor were the 
zenith height 100 km. make the reduction from the observed dura- 
tion the absolute duration defined above use was made the curves 
published McKinley Fig. 5). These curves were extended cover 
ranges less than 100 km. 

All meteors were now sorted into uniformly spaced classes for every 0.5 


MILLMAN AND McKINLEY: METEOR ECHO DURATIONS 


absolute magnitude and every 0.2 log absolute duration. The resultant 
scatter diagram reproduced Fig. will noted that any one class 
absolute magnitude log absolute duration there wide scatter about 


i 


2.8 


2.0 


0.8 


LOG ABSOLUTE DURATION 


3283 METEORS 
UTE UDF 


Fic. Scatter diagram the total observational material illustrating the statistical 
relation between the absolute visual magnitude meteor and the log absolute duration 
seconds its radio echo. 


the mean. This results from combination real differences among the indi- 
vidual meteors with various forms accidental and observational errors. 
can best handled single error function. For the absolute magnitude 
classes the scatter log duration gives evidence double peak maximum. 
The origin this due primarily the non-homogeneity the observational 
material, the bulk which made two categories, Perseid and Non- 
shower meteors. These have somewhat different duration characteristics and 
most the tendency double maximum for log durations disappears when 
the material separated into smaller, more homogeneous groups. view 
this, attempt was made reduce further the data whole. The four main 
groups, Perseids, Geminids, Aquarids, and Non-shower, were handled in- 
dependently. Tables list the basic observational data for each 
these groups, giving the numbers meteors sorted into the classes abso- 
lute magnitude and log absolute duration. 
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TABLE IV. (256 METEORS) 


Absolute magnitude 
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has been noted above that the relation previously published (Millman 
1950) was found averaging all log durations for each magnitude class. This 
method averaging decreases the slope the curve. The effect can best 
seen giving illustration. Fig. the Perseid data have been plotted 
two ways, (A) averages log duration for each magnitude class—filled circles 

28 


2.4 


20 


0.8 


04 


LOG ABSOLUTE DURATION 


0.0 


1404 PERSEID METEORS 


-6.0 -3.0 -40 -3.0 -2.0 1.0 2.0 3.0 4.0 5.0 
ABSOLUTE MAGNITUDE 


Fic. Mean values for the Perseid meteors: 
A—mean log absolute durations for absolute magnitude classes. 
B—mean absolute magnitudes for log absolute duration classes. 


and solid line, and (B) averages magnitude for each log duration class—open 
circles and broken line. will seen that there marked difference be- 
tween the two curves. Obviously the best statistical relation use any 
theoretical application studies the upper atmosphere something 
between these two extremes. 

preliminary finding the best mean curve between absolute magnitude 
and log absolute duration, the error functions for scatter the two coordinates 
were studied for the Perseids and for the Non-shower meteors. was assumed 
that the form the error function was constant over the entire range for each 
coordinate and that the observational material was essentially complete over 
approximately the central third these ranges. the basis these assump- 
tions mean error functions were determined and are listed Table VII. These 
are given terms percentage the observed material each class listed 
relation the mean class. The figures are rounded off the nearest half 
one per cent. Owing the greater homogeneity the data the mean errors 
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TABLE VII 
Error functions absolute magnitude Error functions log absolute duration 
mag. Perseids Non-shower log dur. Perseids Non-shower 
0.5% 1.4 0.5% 
1.5% 2.5 1.0 1.5 1.5 
4.0 4.5 0.8 4.5 4.5 
8.5 7.0 0.6 8.5 9.0 
13.0 10.5 0.4 11.0 15.0 
16.0 14.0 0.2 13.5 17.0 
0.0 17.5 17.0 0.0 16.5 14.5 
+0.5 14.5 15.0 15.5 11.5 
+1.0 10.0 12.0 13.0 9.0 
+1.5 7.5 8.0 7.5 6.5 
+2.0 4.5 4.5 4.0 4.5 
+2.5 2.0 2.5 2.0 3.0 
+3.0 1.0 0.5 1.0 2.0 
1.0 
0.5 
100.0% 100.0% 100.0% 100.0% 


Mean error 0.91 mag. 0.98 mag. 0.38 log dur. 0.41 log dur. 


for the Perseids are slightly smaller than for the Non-shower meteors but the 
difference not great one would have expected from other considerations. 
The error functions determined from the 3283 meteors taken together give 
mean errors somewhat higher than those listed Table VII. 


Mean errors 


Perseids +0.91 mag. +0.38 log dur. 
Non-shower +0.98 +0.41 
All meteors +1.06 +0.43 


little difficult say how much these mean errors intrinsic the 
meteors themselves and how much observational error. The mean error 
the original magnitude estimates may taken +0.5 magnitude. 
consider the various reduction and reading errors can place the observa- 
tional mean error single absolute magnitude determination between 
+0.6 and +0.7 magnitude and log duration between +0.2 and +0.3. This 
leaves +0.7 magnitude and +0.3 log duration the scatter originally 
present the meteors. These values should regarded only rough approxi- 
mations. 

Assuming the Perseid error functions Table VII and using the methods 
numerical integration, correction was computed for each the points 
Fig. The resulting corrected positions have been plotted Fig. significance 
the data being indicated three sizes dots based numbers meteors 
contributing each point. Although there very slight tendency for the 
points little low for the very bright and very faint meteors seems 
best represent the relation straight line was done previously. The 
points the extremes this line are represented relatively few meteors and 
the magnitude estimates these meteors are difficult make and are more 
subject systematic error. 
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Fic. Mean values for Perseid meteors after correction for the effect the error function. 


The observational material available for the Geminid and Aquarid meteors 
was not enough make possible the determination reliable error functions. 
Hence, the error functions for the Perseids were applied the case the other 
two showers and the mean points corrected the same way for the Perseids. 
The Non-shower meteors were reduced, using their error functions listed 
Table VII. The best straight lines were now drawn through the corrected 
points for the Geminids, Aquarids, and Non-shower meteors. These, together 
with the Perseid line from Fig. are plotted Fig. The equations for these 
lines are given below where the absolute duration seconds and the 
absolute magnitude. 


(2) log 1.15—0.349 
(3) log 1.21—0.322 
(4) log 1.10—0.378 


the case the Geminid and Aquarid curves there were not enough obser- 
vational data make good determination the relation and the results for 
these showers should considered tentative. 

For comparison purposes the Perseid curve published Millman (1950) 
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LOG ABSOLUTE DURATION 


P PERSEIDS (1404) 
G GEMINIDS (256) 
O 8 AQUARIDS (158) 
N NON-SHOWER (1420) 
PERSEID CURVE 


-5.0 -4.0 3.0 -20 -1.0 1.0 2.0 3.0 4.0 5.0 
ABSOLUTE MAGNITUDE 


Fic. Corrected mean relations between absolute magnitude and log absolute duration 
for various groups meteors. 


has been shown Fig. dotted line. The smaller slope this earlier curve 
due the fact that only magnitude classes were averaged without any cor- 
rection for the error function. Its position above the final curve because 
was drawn for apparent magnitude. Conversion absolute magnitude would 
shift approximately one magnitude interval the left. 

will noted that the curves Fig. are fairly close together for the bright 
meteors. absolute magnitude the over-all average duration for all types 
meteors 850 sec. move the fainter meteors the echo duration for 
corresponding visual brightness longer for the slower objects. particular, 
the Geminid meteors have smaller slope than the Perseid meteors, that 
faint Geminids have echoes which last considerably longer than faint Perseids. 
This apparently anomalous effect very likely due the fact that the 
average the enduring ionization occurs lower height for the slower meteors 
(Kaiser 1953; Millman 1950, 1954), and hence subject smaller diffusion 
coefficient. 

The material dealt with this paper consists almost entirely overdense 
ionization trails where the number electrons per cm. path length, 
greater than (Greenhow, 1952). this case, according Greenhow, the 
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where the diffusion coefficient and the critical electron 
density for the wavelength used. McKinley has shown (1953, Fig. that there 
second order correction the value for given equation (5) depending 
the sensitivity the transmitter—receiver system used the observations. 
will assume that equation (5), without correction, applies observations 
made with the Ottawa equipment. 

Greenhow and Neufeld (1955) have recently determined the variation 
diffusion coefficient with height study underdense meteor trails. The 
mean height Perseid echoes measured Ottawa three-station triangula- 
tion 99.3 km. and Geminid echoes 91.6 km. (Millman 1954). The diffusion 
coefficients (Greenhow and Neufeld 1955) corresponding these heights have 
ratio whose log 0.53. For observations made single wavelength, equa- 
tion (5) may written 
(6) log log a—log 


allow for the differences the diffusion coefficients the Geminid curve 
Fig. should lowered 0.53 for comparison with the Perseid curve. This 
brings the Perseid and Geminid curves into coincidence near absolute magni- 
tude and can conclude that the ratio between ionization and visual 
luminosity roughly equal for the two showers this magnitude. Making 
reasonable assumption concerning the mean velocity the Non-shower 
meteors, this curve also, when corrected for change the diffusion coefficient 
with height, agreement with the other two magnitude +5. The 
curve does not agree well but based few meteors that 
not very well determined, particularly the lower end. 

Assuming equation (6) correct, and taking average for any given 
shower, log varies directly log Since each drop magnitude one 
slope —0.400 between log and indicates constant ratio between ioniza- 
tion, and visual luminosity for all magnitudes. seen from equation (1) 
this very nearly satisfied for the Perseid meteors. For our observations 
hence from equation (5) 


take 100 km. the best general average for the height Perseid echoes, 
zero magnitude Perseid meteor from equation (7), round 
numbers 

Since the Perseid curve shows nearly constant ratio between ionization 
and visual luminosity can adopt for Perseids absolute 
magnitude +5. This value can also adopted good average for the 
Geminid and Non-shower meteors. brighter meteors the corrected 
Geminid curve falls below the Perseid curve and magnitude the relative 
ionization for the Geminids has dropped factor between and when 
compared with the Perseids. These results alter somewhat the conclusions 
Greenhow and Hawkins (1952), which were based our preliminary experi- 
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mental curve. The slopes the final curves given here are fair agreement with 
the theoretical curve plotted Greenhow and Hawkin’s paper (1952), and 
the indication that there are more electrons present for the brighter meteors 
than was suggested their paper. upward revision the number 
electrons present bright meteor trails will simplify the explanation certain 
features meteor spectra. Further consideration this point will appear 
elsewhere. 


Among those who assisted various times the reduction basic data, 
the authors would like express their thanks particularly Miss Joan 
Stewart, who was chiefly responsible for compiling Tables III VI. 
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MODIFICATION THE ONE-PARTICLE FORMULA FOR 
NUCLEAR 


WALLACE 


ABSTRACT 


The one-particle theory isomeric y-transitions nuclei consistently over- 
estimates the probability per unit time decay. shown that, modifying 
the one-particle model take account the readjustment the wave functions 
all the nuclear particles during the predominantly one-particle transition, 
satisfactory agreement with experiment can obtained. The fact that odd- 
proton transitions are depressed more than odd-neutron ones attributed the 
exchange character the force giving rise the readjustment. 


INTRODUCTION 

The nuclear shell (one-particle) model nuclei has been used with consider- 
able success basis for the analysis isomeric y-transitions nuclei 
(Blatt and Weisskopf 1952, Chap. and Appendix Moszkowski 1953; 
Stech 1952; Weisskopf 1951). has however been noted, especially Gott- 
fried (1953), that the single-particle model tends overestimate the transition 
probability sizable factor. magnetic transitions, for which 
the whole both the data and their interpretation are most clear-cut, this 
factor 100 most cases. The discrepancy between theory and 
experiment undoubtedly significant, and requires explanation. 

has also been noted Graham and Bell (1953) that the transition 
probabilities for odd-proton transitions fall further below the predictions 
the single-particle model than those for odd-neutron transitions. study 
Table 3-1 Gottfried’s thesis shows that this also the case, the average, 
for transitions. 

explanation these two discrepancies offered, along the following lines: 
when one nuclear particle makes transition, the wave functions all 
the other particles the nucleus undergo small change, due the contribu- 
tion the potential field which they move. the final-state wave 
function the particle and its initial-state wave function, and 
the transition probability thereby decreased factor 


(1) 


where the integration over the coordinates the particle. Although indivi- 
dually the f’s may small, their large number heavy nucleus reduces the 
transition probability just such factor that noted. The greater effect 
odd-proton transitions may attributed the tendency, due exchange 
force, for protons interact most strongly the average with neutrons, and 
vice versa. The greater numbers neutrons will then accentuate the effect 


’ 


the odd-proton case.* 


received August 11, 1955. 

Contribution the Department Mathematics, McGill University, Montreal, Quebec. 

*It was once suggested the author (Bull. Am. Phys. Soc. (No. that the 
explanation lay the effect the Coulomb force. now appears that this conclusion was 
based error numerical calculation, and that the Coulomb force somewhat too weak 
for the purpose. 
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true that the readjustment the wave function all the nuclear 
particles gives additional contribution the transition probability. This, 
however, may seen small the cases which the transitions are 
ones. For suppose that the average factor which the 
readjustment single particle cuts down the matrix element 
Then the major contribution decreased factor being the 
number nuclear particles. But the case the additional contributions the 
matrix element is.in the first place factor must then 
very small. These additional contributions are then times smaller 
than that from Since see that they may neglected. 

These considerations apply the case which the transition 
single-particle. They will not apply the very strong transitions, 
nor those cases which the correct initial and final spin assignments require 
coupling the angular momenta two more particles (e.g. transitions 
odd nuclei). 


MODIFICATION THE SINGLE-PARTICLE FORMULA 


Let designate the initial- and final-state wave functions the particle 
which makes the transition and Let also designate the 
initial-state wave functions the other nucleons ya. 

Now the change the potential due the transition may written: 

where have assumed Yukawa interaction between particles. the 
strength the interaction and the Compton wavelength, 
i.e. half the classical electron radius. the Coulomb force 
contribution which shall see too small much significance. 

With regard the exchange character the force, course not well 
enough known taken into account quantitatively. should presumably 
taken have somewhat different values depending the relative sym- 
metry the wave functions and the core particle.* will seen (equa- 
tion (6) below) that the contributions from the various core particles appear 
additively exponent. shall, for first orientation, use single (average) 
Vo, which will carried through the calculations parameter. will 
found that best agreement with experiment will then obtained with 
Mev. with Mev. for the deuteron (Sachs 1953)). 
shail then consider the effect assuming the force mixture 
the ordinary and Majorana forms. 

may now apply simple perturbation theory discuss the perturbation 
the particle wave functions. have for the final-state wave function 
corresponding the initial state 


*Since the result will depend only (equation (9)) the sign therefore 
concern for our purposes. 


i 
é 
iy 
t 
| 
t 
i 
é 
— 


CANADIAN JOURNAL PHYSICS. VOL. 34, 1956 


being the matrix element between the state and any other state 
the potential. and are the energies the states and 
respectively. The factor must chosen normalize i.e. 

H. 
follows that the readjustment factor arising from the particle 
may therefore put 


( ) f (Ea—Es)° 

But the product all the factors 

where 

( ) (Ea— Es) 


The summation over occupied states, the one over all states the nuclear 
potential. 


METHOD CALCULATION 


would clearly unduly tedious for our purposes evaluate the sum- 
mation (7) for the aggregation of, say, spherical well oscillator states. 
the other hand, important take account the explicit form the 
wave functions the particle shall use, for these, harmonic oscillator 
states. For all the other particles and states shall simplify matters using 
the Fermi gas approximation. 

One must take care, however, handling the denominators 
There will, course, finite spacing between and the nearest available 
level The individual features various nuclei will fact manifested 
this point, through the details level structure. order, however, get 
reasonably accurate statistical result, shall replace the summations 
integrations, but omit from the region energy 
either side being the density energy levels the Fermi gas 
energy 

have, therefore, for equation (6) 


where the nuclear mass, and are the wave-numbers the states 
and and the nuclear volume. The k-integrations are subject the 
limitation mentioned the previous paragraph, and the r-integration runs 
only over the interior the nucleus. However, this also the range over which 
has appreciable value, and consequently serious error will 
involved extending the integration over all space. any case, the ‘‘radius 
the indefinite within the distance which falls off. 
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Making the approximation question, then, see that our result involves 
only the Fourier transform H’(r). Since convolution 
and the interparticle potential, its Fourier transform the 
product the Fourier transforms these quantities. Writing x(k) for the 


must borne mind that the integration such that 


evaluate the integral (9), introduce new variables (with range 
km, Corresponding the highest filled state), and the 
range which limited 


now easy the integration over get 


The latter two terms the square bracket the integrand are negligible 
compared with the first; for 


(14) 


where the total number particles the nucleus. comparison, the 
term k,,/2K unimportant owing the fact that, for small goes 
like that this region contributes very little the integral. The logarith- 
mic term also not important. Its singularity integrable, and any case 
occurs for 2km, which value the rest the integrand completely 
negligible. These terms have fact been included the numerical integration 
which follows, and are found not influence the result significantly. 

With respect the leading term, should noted that might have 
distinguished the contributions neutrons and protons the nucleus through- 
out. However, these contributions would additive. Different would 
appropriate for the two cases, but the sum the terms would still 

good approximation have, then, 


The angular integration may performed explicitly, and the K-integration 
may done numericaliy fur each case interest. 
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RESULTS 


shall confine our attention transitions. fact, should not 
expect any serious differences between the various cases: 

general, would necessary distinguish the factors for transitions 
involving the various states the particle For however, this 
not necessary. For fact the (magnitude)? the functions spherically 
symmetric, being given 


2 


Since only relative orientations initial- and final-state wave functions are 
important, all give the same result. take, therefore, the case 


9/2, and then find 


where (r, are polar coordinates. 
may now evaluate the Fourier transforms these functions get 


and 


being the polar angle the transform variable 

The first these wave functions corresponds the oscillator level 
the second these are the levels appropriate the known transitions 
this type. The units the wave chosen numerical 
coefficients simple possible. However, must determine their relation 
more natural units. choose the inflection point the (wave function)? 
the outer particle the boundary the nucleus, this yields nuclear radius 
our units 2.45. Now the nuclear radius conventional units 
our units. For the known transitions runs from 
115 therefore ranges from 1.67 However, this variation does not 
affect the result very much, shall consider only average value 1.9 

may now evaluate analytically the angular integrals (15), after which 
the K-integration can done numerically. obtain the result: 


climb after 0.8, reaches peak about 1.6, and becomes negligible 
after 3.2. our units and hence, remarked earlier, is, 
the whole significant range, and outside the significant 
range. Thus the neglect the latter two terms (13) was justified. 
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For the isomers this type listed Gottfried, the average 
reasonable value. 


CHARGE EFFECT 


have estimated the contribution from the Coulomb force, ignoring 
“interference” effects between Coulomb and nuclear force. This probably 
justified because the Coulomb force has its greatest effect the range which 
the nuclear force contribution insignificant, and also because the nuclear 
force must, for saturation, have one sign almost often the other. The 
Coulomb contribution then follows from (15) putting 


For the case, 50, this gives only 0.1, which too small 
serious factor. 

However, another factor giving rise difference between the odd-neutron 
and odd-proton cases the effect the exchange character the forces. Let 
assume that the nuclear force mixture ordinary and Majorana force. 
The latter will repulsive between antisymmetric pairs, attractive between 
symmetric pairs. Thus, may have two values, and V,, according the 
core nucleons form symmetric antisymmetric combinations with the outer 
particle P,. Now proton, the interaction will with one-quarter 
the core protons, with three-quarters them. For core neutrons, the 
the ratio g/A (15), now have for the odd-proton case 


while the odd-neutron case 

The average these 

this reflects the result, pointed out Blatt and Weisskopf (1952, 147), 


that there are the nucleus more antisymmetric than symmetric pairs the 
ratio 5:3. The difference the factors 
2 2 
considerably bigger than V,,* this ratio becomes 
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For the cases under consideration, this gives the average odd-proton 
transitions about five-eighths fast comparable odd-neutron ones. the 
average cases considered Gottfried, found that they are little more 
than half fast. 

may, course, argued that experimental errors make these small 
differences insignificant. content ourselves with pointing out difference 
factor almost two the data averaged over nearly transitions, 
and noting that some such effect might expected theoretically. 


COLLECTIVE MODEL 


the collective nuclear model Bohr and Mottelson (1953), reduction 
the single-particle transition matrix element due readjustment the 
nuclear surface expected, they point out (p. 107). However, obtain 
agreement with the observed results, necessary postulate that the initial 
and final surface wave functions are almost orthogonal, which seems unlikely 
even the case strong coupling. The difficulty getting large enough re- 
adjustment factors this model arises from the suppression large number 
degrees freedom, one which would contribute seriously, but which, 
taken together, would have considerable effect. 
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THE DECAY AND 


ABSTRACT 


The energies and intensities gamma rays have been measured 
studying their external conversion spectra with high resolution beta-ray 
spectromeier. analysis the beta-ray spectrum reveals groups with end 
points 2116 kev. (79%) and 1961 kev. (20 and several low energy groups 
whose total intensity the order 1%. The observed data can consistently 
accounted for terms decay scheme involving the following levels 
0.0, 155.0, 633, 1086, 1306, 1461, 1765, 1941, and 1958 kev. This scheme 
firmly pported The spins the first three 
listed are 2+, and and that the 1461 kev. level probably 


The decay accompanied gamma rays the following energies 
and quantum intensities: 122.9 kev. (0.6%), 137.2 kev. (10%), 630.8 kev. 
(0.024%), and 768.2 kev. (0.024%). These values are good agreement with 
those quoted the literature, except for the intensities the last two radiations. 


INTRODUCTION 


The reaction rhenium leads 91-hr. and 18-hr. The 
activity has been extensively studied (Beach al. 1949; Grant and 
Richmond 1949; Guss al. 1954; Koerts 1954; Langer and Price 1949; 
Metzger and Hill 1951; Porter 1955; Steffen 1951) and the measurements 
given here serve chiefly confirm existing information. The 18-hr. activity 
accompanied rich spectrum beta- and gamma-radiations which 
were only partially resolved workers using instruments moderate resolu- 
tion (Beach al. 1949; Cork al. 1948; Goodman and Pool 1947; Miller and 
Curtiss 1946; Richmond al. 1952). This paper presents the results 
series experiments carried out with high resolution beta-ray spectrometer 
(Johns al. 1953) and coincidence circuits and proposes decay scheme well 
supported coincidence measurements. Some these results have already 
appeared brief form (McMullen and Johns 1953 


APPARATUS AND METHODS 


Nine sources spec-pure rhenium powder, sealed aluminum quartz, 
were irradiated the Chalk River Brookhaven reactors for use gamma 
sources. These were shipped air express and had effective strengths 
about 300 iac. upon arrival. The gamma radiations from these sealed sources 
were studied examining the external conversion peaks which they produced 
lead, uranium, gold radiators thicknesses from 0.4 mgm./cm.? 
The annihilation radiation from the small amount copper impurity 
aluminum proved somewhat troublesome early and made necessary 
shift quartz containers. All the radiations observed decayed with either the 
18-hr. the 91-hr. half-life expected, indicating that they are associated with 
either the the nuclides. 

Sources less than mc. activity sealed quartz vials were used for 
Preparation beta sources for coincidence work. Such quartz vials were 
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crushed under nitric acid and the resulting rhenium nitrate evaporated 
dryness the desired source backing. The beta sources used were sufficiently 
thin that the conversion peaks the kev. region possessed full width 
half maximum less than momentum and showed serious degrada- 
tion the low energy side. 

The Geiger counter detector used the earlier runs (Johns 1953) was 
replaced anthracene crystal photomultiplier combination for 
later measurements. This change made possible the 100% detection elec- 
trons energy down kev. without difficulty and removed the necessity 
for making counter-window absorption corrections for conversion lines the 
kev. region. 

Details concerning the calibration the spectrometer for energy and 
intensity measurements have already been described (Johns and Nablo 1954). 
External conversion peaks produced the precisely known gamma rays 
and and internal conversion lines from the thorium active- 
deposit were used standards. The beta end points were obtained extra- 
polation from the Fermi analyses. The errors quoted all energy measure- 
ments are believed include both statistical and systematic errors. 

Two gamma-ray scintillation spectrometers using anthracene 
crystals, 6292 photomultipliers, and single-channel pulse-height analyzers 
(Bell and Graham 1953) were used coincidence circuit resolving time 
With one analyzer set record the desired portion the beta spectrum 
the photopeak one gamma ray, the gamma-ray spectrum was scanned with 
the other analyzer and both the single-channel and the coincidence spectra 
were recorded. From comparison the two, was generally possible 
draw significant conclusions concerning the decay scheme. 


EXPERIMENTAL 


The gamma rays observed are described Table (a). The energies 
the 123 and 137 kev. radiations were found both external and internal 
conversion measurements. The intensities were all measured relative the 
137 kev. gamma ray whose absolute intensity was found manner 
explained later. 

The intensities the 631 and 768 kev. gamma rays relative the 137 kev. 
radiation were measured two independent methods. The first used the 
Siegbahn spectrometer make comparison the photopeak the 137 
kev. gamma ray 1.8 mgm./cm.? uranium radiator with the photopeaks 
the high energy radiations mgm./cm.? uranium radiator. The 
calculation intensity requires the use semiempirical formula relating 
peak height radiator thickness and electron energy (Johns and Nablo 1954) 
and the assumption that the K-shell conversion efficiency five times that 
the L-shell. was impossible get useful comparison through peaks 
alone because the serious energy degradation the low energy peak. This 
method gave intensity 0.02% for the two high energy radiations. 
the second method, the photopeaks the three gamma rays under question 


TABLE 


(a) GAMMA RAYS 
Quantum intensity, Gamma-ray energy, Detection method and comment 
kev. 
122.9+0.4 1.8 mgm. and internal conversion 
137.2+0.2 1.8 mgm. and internal conversion 
0.024+0.010 630.8+0.8 mgm. coincidence with 137 
0.024+0.010 768.2+1.0 mgm. nocoincidence with 137 
(b) BETA GROUPS 
(c) INTERNAL CONVERSION DATA 
123 kev. 0.32% 0.74% 0.43+0.2 
137 kev. 4.65% 7.6% 1.4% 0.45+0.03 


were compared scintillation spectrometer whose efficiency 
function energy had previously been determined. This method gave 0.03%, 
reasonable agreement with the 0.02% value. 

Fermi analysis the beta spectrum, illustrated curves and 
Fig. showed the clear-cut break 137 kev. from the end point obtained 


Fic. Fermi analysis the beta spectra and Curves and refer 
and The curvature probably represents departure from allowed shape. 
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others (Guss 1954; Langer and Price 1949) with the values for the energies 
and branching ratios given Table 

The conversion electron spectrum shown Fig. and the infor- 
mation obtained therefrom Table I(c). quite clear from the values 
the K/L, and the ratios that both the 123 and 137 kev. radiations 
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Fic. Internal conversion peaks and The 123 and 137 kev. peaks are 
associated with the 155 and 478 kev. peaks with With expanded momentum 
scale the peaks can resolved into components. 


are character. Using the scintillation spectrometer, was found that the 
number K-shell X-rays from tungsten and osmium per 137 kev. quantum 
was Using the two scintillation spectrometers, the number 
coincidences per coincidence was found 0.43+0.03. The latter 
ratio when corrected for the fluorescence yield osmium gives value 
0.45+0.03 for the K-shell conversion coefficient the 137 kev. transition 
The difference between the measured ratios, similarly corrected, 
leads value 0.34+0.05 for the number vacancies occurring 
tungsten per 137 kev. quantum. 

From the values and given Table I(c), appears that the 
quantum intensity the 137 kev. gamma ray 10.3%. This value, combined 
with the 0.34 ratio above and corrected for the number X-rays following 
K-shell conversion the 123 kev. gamma ray, leads K-capture branching 
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ratio 3.2%, with roughly half the transitions going the 123 kev. level 
and the other half the ground state. The 10.3% figure, when added the 
number and electrons, gives intensity 24% for the 927 kev. 
beta group, excellent agreement with the values found from the Fermi 
analysis. Gamma-gamma coincidence experiments showed clearly that the 
631 kev. radiation coincidence with the 137 kev. gamma ray, while the 
768 kev. radiation not. 


EXPERIMENTAL 
(a) Gamma Ray Energies and Intensities 


The external conversion spectrum presented Fig. All the 
peaks were recorded least five times. addition, the reality the weak 
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Fic. Composite the external conversion spectrum Curve was obtained with 
mgm./cm.? lead radiator, curve with uranium radiator, and curve 
with uranium radiator. 


high-energy peaks was confirmed removing the radiator and examining 
the underlying Compton distribution. 

The energies and intensities the gamma rays shown Fig. are tabulated 
Table The errors quoted are twice the standard deviation the mean, 
and include the effects random errors the determinations the peak 
position and systematic errors due calibration. 

The relative intensities the gamma rays were measured with respect 
the 633 kev. radiation whose A-shell photoelectron peak stands cut clearly 
the spectrum. 
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TABLE 
GAMMA RAYS 


cent energy, kev. kev. 
9.0 155.0+0.2 0.0 All but 633 and 1959 kev. 
0.04 453.6+0.6 633 kev. (uncertain) 
1.00 633.0+0.1 +0.1 1133, and part 1308 kev. 
0.093 633 kev. (uncertain) 
0.33 828.0+0.2 +0.1 155 kev. 
0.012 1959 8-0 


Some difficulty was experienced comparing the intensity the 155 kev. 
radiation with that the 633 kev. gamma ray because the K-photoelectron 
peak for the former falls region where energy degradation the radiator 
serious. This comparison was made using Au, Pb, and radiators and also 
comparing the peak intensities, since these are not subject the 
effects radiator degradation. addition, the 155 kev. 633 kev. intensity 
ratio was measured means the scintillation spectrometer and 
found give value 8.9+0.5, good agreement with the mean value 9.0 
obtained from the beta-ray spectrometer determinations. Although the 
intensity estimates Table are believed self-consistent within 20%, 
the absolute values may error somewhat larger amount since these 
depend absolute intensity determination the 155 kev. radiation 
method described later. 

intensive search was made for possible cross-over transitions predicted 
from decay scheme constructed the basis the stronger lines. This search 
resulted the discovery the 453.6 kev. radiation, and set upper limits the 
intensities three other possible transitions follows: 375 kev. (<0.01%), 
1086 kev. (<0.01%), and 1460 (<0.02%). 


(b) Coincidences 

Typical coincidence spectra are shown Figs. and With channel 
set the photopeak the strong 155 kev. radiation, both the single-channel 
and the coincidence spectra were recorded. clear from Fig. that the 478 
kev. radiation coincidence with the 155 kev. gamma ray and that the 
633 kev. line not. This, together with the fact that 478.2+155.0 633.2, 
indicates that the first two excited states lie 155 and 633 kev. above 
the ground state. Since these two states are connected means the 478 
kev. transition, may expected that all gamma rays coincidence with 
the 633 kev. radiation will also appear coincidence with the 155 kev. 
gamma ray with reduced intensity (0.32 according Fig. 7). 
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single-channel spectra; curve coincidences with the 155 kev. gamma ray. 


Although the 828 and 931 kev. photopeaks are poorly resolved the single- 
channel spectrum, clear from the coincidence spectrum that the 931 kev. 
radiation full coincidence with the 155 kev. line. the other hand, the 
828 kev. radiation feeds the 155 kev. level indirectly, since its intensity the 
coincidence spectrum has been reduced roughly half that the 931 kev. 
peak. 

The weak 650 kev. peak which appears the coincidence spectrum 
partially due coincidences with the Compton pulses from high energy 
gamma rays appearing Channel but also seems include coincidences 
with the weak 673 kev. radiation, since all cases the peak seems shifted 
energy higher than 633 kev. 

Fig. shows the high energy portion the total spectrum recorded 
Counter together with the coincidence spectra associated with the 155 and 
633 kev. gamma rays Counter Counter was covered with mm. 
lead permit the use stronger sources without pile-up. For coincidences 
with the 633 kev. radiation, Counter was also covered with equal amount 
lead. With pile-up eliminated, the 828 and 931 kev. peaks are more clearly 
resolved than Fig. The 155 kev. radiation appears coincidence 
with all the high energy radiations, except the 1959 kev. gamma ray. Fig. 
shows clearly that the 828, 1133, and 1308 kev. radiations are cascade with 
the 633 kev. gamma ray and that the 931 kev. radiation not. 

was noticed that the photopeak observed for the 1308 kev. line the high 
resolution beta-ray spectrometer was consistently too broad, indicating that 
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Fic. scintillation spectrometer spectra above 600 kev. Curves and 
single-channel curves and coincidences with the 155 kev. radiation; curve 
coincidences with 633 kev. radiation. 


this line really close doublet. the coincidence curves Fig. evident 
that the ratio the 1133 the 1308 kev. photopeak changes from 0.6 for the 
155 kev. spectrum 2.0 the 633 kev. spectrum. This behavior 
expected the 1308 kev. radiation consists pair gamma rays roughly 
equal intensity and the same energy, one which feeds the 155 kev. level 
directly, while the other member the pair and the 1133 kev. radiation feed 
this level indirectly via the 633 kev. level. 

The extension curve Fig. down 300 kev. not shown the 
curve presented. This region was examined and found show weak peaks 
about 460 and 650 kev. These peaks are chiefly due coincidences between 
the 478 633 kev. radiations Channel and pulses the Compton 
distributions high energy gamma rays appearing channel From the 
apparent shifts these peaks, feel there weak evidence for 453-633 kev. 
and 673-633 kev. coincidences. 

The results these experiments are summarized the last column Table 
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(c) The Beta Spectrum 

Three separate beta sources were analyzed. This analysis was quite straight- 
forward for energies above Mev., but below this point was complicated 
the presence conversion lines and the necessity subtracting 
contribution which amounted much half the total counting rate 
some points the spectrum. Any systematic error, introduced subtracting 
this contribution, would seriously affect the shape the beta distribution 
below Mev. 

Curves and Fig. present the Fermi analysis one the sources 
studied. All showed the characteristic break point 155 kev. below the 
upper end point. Also, all showed slight tendency concave the energy 
axis. Since the amount this curvature was not strictly reproducible 
different runs, felt that might have been partly due errors sub- 
tracting the component. However, slight effect always appeared, 
above Mev., suggesting that this spectrum exhibits the type forbidden 
shape recently found (Porter a/. 1955). 

The end point the high energy group 2116+2 kev., where the error 
quoted represents the maximum deviation from the mean obtained from the 
three independent determinations. 

Instead using the Fermi analysis determine the end point the lower 
energy beta group, was assumed that this group fed the first excited state 
155 kev. and therefore must have end point 1961+2 kev. The 
analysis was then used establish the intensities the two groups, since 
was clear from the gamma measurements that all the other beta groups 
combined would not give contribution more than 1.5°% the total beta 
spectrum. the basis these runs, the intensities the 1961 and 2116 kev. 
groups are, respectively, and 79+4%, with roughly allowed for 
the lower energy groups. These intensities were determined the assumption 
that both spectra showed forbidden shape but the departures from the 
allowed shapes are small that the intensity values are not affected this 
assumption. 

Since both the beta groups represent first forbidden transitions 
yes) they may expected show small forbidden shape factors. The total 
beta spectrum may described the function 


where and are the energy-dependent shape factors for spectra 
with end points and and intensity parameter for the second 
group. calculation these factors (Banerjee and Saha 1954) for mixture 
the scalar and pseudoscalar interactions with yes, shows that 
and are the form where the parameters 
are the order 0.1 less. For our purpose was convenient re- 
arrange Equation (1) 


where and are the order 0.1 less. 
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practice was possible obtain fit the data from each source 
setting and using values and less than 0.1. Since the cal- 
culations involve great many terms which are not simply related the 
constants beta decay, and since this method does not uniquely separate the 
two spectra, seemed pointless carry this analysis further. 


TABLE III 
INTERNAL CONVERSION DATA 


Experimental values Theoretical values 


155 kev. transition 


=3.0 +0.4% 

=3.8 

=0.9 +0.15% 

K/L 0.79+0.15 M1, 7.7; 3.5; E2, 0.7 

1.00 M1, El, 4.5: 1.25: 1.0; 

478 kev. transition 


*Goldhaber and Sunyar 1951. 
al. 1952. 
al. 1951. 


(d) Internal Conversion 

The internal conversion spectrum shown Fig. almost associa- 
ted with the 155 kev. transition. The electron conversion each 
shell expressed Table III. The three components the conversion 
electrons were sufficiently well resolved allow the 
L,: ratio. 

The value might have been calculated from the data Table 
and the intensity the 1961 kev. beta group found from the analysis. 
Since this method very sensitive errors the latter quartity, the value 
was measured directly the technique used for the kev. radiation 
With the beta counter set detect only betas greater than 
all interference from was eliminated. clear from com- 
parison the experimental and theoretical values Table that the 155 kev. 
transition character. This assignment agreement with the con- 
clusions Mihelich (Scharff-Goldhaber 1953), who has also measured the 

The data presented for the 478 kev. transition, though not extensive, 
indicate that must predominantly character, though some mixture 

(e) Absolute Intensity the 155 kev. Gamma Ray 


From the measured values and Nx, the absolute intensity the 155 
kev. gamma ray Since the total number conversion electrons 
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associated with this transition 8+1%, the 1961 kev. beta group must 
represent 16+3% the disintegrations. This value reasonable agreement 
with the 20% figure found from the Fermi analyses, and leads mean value 
18+3% for this intensity. 

The best value for the absolute quantum intensity the 155 kev. gamma 
ray from these data has been taken half the 18% figure This 


has been used express all the gamma intensities Table 


absolute scale. 


(f) Angular Correlation 

Before the complexity the gamma-ray spectrum was realized, angular 
correlation studies were made, using anthracene crystals covered with mm. 
lead and fast coincidence circuit resolving time sec. Since there 
was energy selection this arrangement, except for the exclusion the 
155 kev. radiation, the apparatus was detecting the following cascades: 
828-633, 828-478, 1133-633, 1133-478, 1308-633, and 1308-478, which the 
first two are much the strongest. The strongly anisotropic angular correlation 
pattern obtained was given the equation 
unique interpretation these data does not seem possible, but clear 
that some strongly anisotropic cascades such 0-2-0, must 
involved. 


The experimental results reported for are summarized the decay 
scheme Fig. This scheme close agreement with the results Metzger 
and Hill (1951) and Steffen (1951), except few minor points. 
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Fic. Decay scheme 


The most serious discrepancy lies the intensity the two high energy 
gamma rays, our value being only one eighth the value, 0.2%, quoted 
Metzger and Hill. Their intensity would have required that the peaks 
these radiations would comparable height the 828 kev. peaks 
when our sources were first received. Actually, neither radiation was 
detectable until the activity had decayed away. 

Our value 3.2% for the number disintegrations leading A-capture 
considerably smaller than the value Steffen but reasonable 
agreement with the value given Metzger and Hill. 
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The end point the beta spectrum agrees very well with values obtained 
Langer and Price (1949) and Guss, Killion, and Porter (1954), and the 
branching ratios agree well with those found other workers (Metzger and 
Hill 1951; Porter al. 1955; Steffen 1951). found evidence for the 
departures from the allowed shape found Porter (1955) but our 
method not nearly sensitive such effects the coincidence tech- 
nique which they employed. 


The experimental results presented for can readily accounted for 
terms the decay scheme Fig. All the gamma rays assigned Table 
appear this scheme. The agreement between the observed energies and those 
predicted from the decay scheme excellent, the difference 
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1802.7 (0.027%) 
1783.0 
13079 


(0.07%) 
1132.8 (0.03%) 
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Fic. scheme Beta transitions shown solid lines were obtained from 
the Fermi analysis; the others deduced from gamma-ray intensities. 


greater than 0.5 kev. has been pointed out, the 1307.9 kev. line double, 
the unresolved lower component being the 1306 kev. radiation corre- 
sponding the transition. possible that this component also con- 
tains contribution from the transition, though direct evidence for 
this exists. 

The only level the decay scheme which unsupported both energy 
and coincidence measurements level 1305.6 kev. above the ground state. 
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Its assignment has been made purely the basis the weak 633-673 kev. 
coincidences observed. 

From the internal conversion data presented, clear that the first two 
excited states are both character. The log fot values for the 
transitions the 155, and 633 kev. levels are 8.1, 8.6, and 9.5, 
respectively, all characteristic first-forbidden transitions with 
yes. strongly suggests that the ground state spin 
impossible from the data available make unique spin assignments the 
higher lying levels The absence ground state transitions from the 
1086, 1461, and 1765 kev. levels, and the log fot values, 9.2, 8.4, and 8.3, for 
the beta transitions these levels, suggest that their spins may all zero. 
The angular correlation data provide further support this assignment 
the case the 1461 kev. level. 

The decay scheme presented here differs from that our preliminary report 
(McMullen and Johns the addition number new levels and 
the rejection the 828 kev. level. Although the energy fit the earlier 
scheme was excellent, was not supported coincidence experiments and 
impossible with the gamma energies present this spectrum design 
unique scheme energy fit alone. 

The lifetimes quoted that report are too long roughly factor 
(McMullen and Johns Our more recent measurements are consistent 
with the precise values given Sunyar (1954) and are characteristic 
transitions. 

The beta spectrum seems show departure from the allowed 
shape, was first suggested Guss, Killion, and Porter (1954). experi- 
ment such that described Porter al. (1955) for required 
make quantitative test this behavior. 

The log fot values for both and are very similar, with the low 
lying levels having value around 8.0 and the higher levels ranging from 
7.5 9.5. This suggests that the ground states both and have 
the same spin and parity, almost certainly 
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THE SPECTRUM THE VACUUM 


CARROLL? 


ABSTRACT 


The spectrum the radical has been investigated the Schumann region 
using three-meter vacuum spectrograph the fourth order. hollow cathode 
discharge through streaming helium, which trace cyanogen was added, 
was used source. Three new systems bands were found the region 
2100 These involve two hitherto unknown states: state designated 
state designated which gives rise the transition Rotational 
analyses the bands have been made and the molecular constants evaluated. 
The new data, although not decisive, support the higher value (8.2 ev.) for the 
dissociation energy CN. 


INTRODUCTION 


For years only two band systems were known the spectrum CN: 
the system’’, which was investigated Jenkins, Roots, 
and (1932) and the violet system, which was studied 
(1926) and Jenkins (1928). Quite recently, however, Douglas and 
Routly have thoroughly investigated the spectrum the region 
2000-8500 and reported three new states: and These states 
and which lie the near ultraviolet region. These authors also 
suggested that investigation the spectrum should undertaken the 
Schumann with view extending our knowledge the radical 
and the obtaining some information its dissociation energy. 
present renorts the extension observations the spectrum 
1600 and give: analyses three new systems bands the region 
2100 

EXPERIMENTAL 

The source the present experiments was essentially the same 
that described Douglas and Routly (1954, 1955). Helium was streamed 
rapidly through hollow cathode discharge tube and trace was 
added until the red appeared strongly when the discharge was viewed 
through pocket spectroscope. Under these conditions, was found that the 
new ultraviolet bands were most readily observed. The spectrum was photo- 
graphed the eter normal incidence vacuum spectrograph which has 
been described and Herzberg (1954). Ilford plates were used 
record the spectra. 

Sufficient study the fine structure the new bands could 
obtained only photographing the spectrum the higher orders. prevent 
the strong emission longer wavelengths from entering the slit the spec- 
received October 1955. 
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trograph, cylindrical lithium fluoride lens was employed the manner 
described Curry and Herzberg (1934). Although the emission most intense 
the slot the cathode, there was sufficient light from the region the elec- 
trodes render impossible achieve adequate separation wavelengths 
the use foreprism alone. However, using auxiliary plano-convex 
lens form preliminary image the source, satisfactory isolation the 
orders could achieved. This arrangement, devised Dr. Herzberg, will 
described forthcoming paper. 

The analyses described the following sections were made from plates 
taken the fourth order and measured against second-order iron arc lines. 
The wavelengths the lines were taken from the tables (Harrison 
1939) and converted vacuum wavelengths using Edlén’s tables (1953). 
The relative wave numbers unblended lines given set measurements 
are believed correct 0.06 although the absolute accuracy may 
somewhat less than this. Lines separated more than 0.45 were 
distinctly resolved, corresponding experimental resolving power 
about 120,000. 

was not possible eliminate entirely some impurity spectra from the 
source. Thus the bands nitrogen frequently ap- 
peared although they were always very weak. The bands 
were always present considerable strength. These spectra undoubtedly 
arise from the decomposition the the discharge and the subsequent 
formation molecules. Fortunately the new bands lay almost entirely 
regions free from nitrogen spectra. 


the region 1700-2050 six simple bands strongly degraded longer 
wavelengths were observed. These were readily identified transition from 
state designated the ground state The wavelengths the 
heads the bands, their origins, and visual estimates intensity are given 
Table reproduction the band shown Fig. The 2040 


TABLE 
DATA FOR THE SYSTEM 


(1-1) (0-3) (0-4) (0-5) 


*Quoted from Douglas and Routly (1954). 


band had previously been observed Douglas and Routly (1954) who con- 
sidered being most probably due unknown singlet transition 

The results the present analyses are presented Table The doublet 
fine structure not resolved any the bands although its presence 
detectable broadening the lines for higher values. The bands were 
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therefore treated consisting single and single branch and the 
rotational constants were derived the usual graphical method from the 
well-known relationship 


Table III gives the constants for the bands. the and lines the (0-2) 
band are not fully resolved, the data for this band are not accurate for 


TABLE III 
CONSTANTS FOR THE BANDS 


The state The state 


*Quoted from Douglas and Routly (1954). 


the others. Comparison the rotational constants and vibrational intervals 
for the ground state obtained the present work with the results 
previous workers (Jenkins, Roots, and Mulliken 1932 and Douglas and Routly 
1955) shows excellent agreement. The constants for the upper state are 
follows: 

58959.9 


1670.88 


only two vibrational levels are observed for state not possible 
derive vibrational formula. 

state. Fig. shows the difference between the rotational terms calculated 
from the rotational constants and the observed levels determined from both 
the and the branches. The displacement the lines accord with 


x P BRANCH 


O R BRANCH 


Ow o* 
010 oO oO 
§ 10 15 20 25 (J 


0.60 

0.50 

x 

é 


> 
- 
> 
3 
i 
x ‘ 


t 


: H i. ! 


x 
N 
| 
| 


CARROLL: SPECTRUM 


what expected for rotational perturbation and both components state 
must equally affected splitting the lines produced 
the perturbation. therefore seems that the perturbing state must 
species Clearly the origin the perturbing vibrational level lies very close 
that the perturbed level, i.e., 60631 interesting note that 
Douglas and Routly (1955) reported perturbation the state which 
they also attributed state whose energy can estimated about 
60910 Neither the perturbations can interpreted being due 
the higher levels state whose positions are known accurately 
61377.7 the high value the dissociation energy (66000 
assumed, possible that one the perturbations due high 
vibrational level the ground state 38). Unfortunately extrapolation 
the available data the ground state too long confirm disprove this 
suggestion. any event, appears that there least one other state 
which has far escaped direct observation. 


IV. THE SYSTEM 


Amongst the strongest features the spectrum were three bands the 
region 1930-2100 They were triple-headed, degraded longer wavelengths, 
and appeared fairly complex structure. The bands were approxi- 
mately equal intensity and from the vibrational intervals calculated from 
corresponding heads they could readily assigned transition from the 
new state the well-known state. This assignment was confirmed 
the fact that the wave number interval between the first and third head 
each band was close the known doublet splitting for the state. Table 
gives the wavelengths the heads the bands together with the origins. 


TABLE 


2007 2082.01 


The latter were computed from the data for the and states ob- 
tained from other band systems CN. Fig. shows reproduction the 
(1-0) band and can seen that its general appearance typical that for 
transition between state and state. The bands were not an- 
alyzed they would yield new information the states involved. 


Four bands degraded longer wavelengths were observed the regior 
single-headed, each band showed fairly complex fine 
band 1856 was assumed the band, the lower state could 
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identified Another weak feature 1670 was identified the (4-0) 
band. Table gives the head measurements the bands together with visual 
estimates intensity and the origins derived from the analyses (see below). 


TABLE 
DaTA FOR THE SYSTEM 


(0-1) (0-0) (1-0) (2-0) (3-0) (4-0) 


Analyses the four strongest bands showed that the transition was from 
new state, designated the state. The single-headed appearance 
the bands due the fact that the state conforms Hund’s case 
low values, inverted, and spin splitting almost equal that the 
state which also inverted. transition where both states 
are expected. The former gives rise three branches, and while the 
latter similarly gives P2, and branches. Furthermore, the strongest 
branches should the branches. These features are exhibited the repro- 
duction the band Fig. Strictly speaking, each line should 
doublet due the A-type doubling both the and states, but the resolu- 
tion not sufficient show the splitting which would very probably 
less than 0.10 

The analyses the (0-0), (2-0), and (3-0) bands are presented 
Table VI. determine the band origins given Table was 
plotted against which case the intercept the axis 
gives plus terms is. and which may neglected. The 
origin the band was estimated from the head. complete the data, 
the origin the band, which wis not analyzed, was computed. 

The molecular constants derived from the analyses are given Table 
Unfortunately was not possible the spin splitting constants 
for both states from the analyses. follows immediately from the 
Hill-Van Vleck formula, which gives the energy levels for doublet states with 
coupling intermediate between cases and that 
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TABLE 


WAVE NUMBERS LINES THE SYSTEM 


576.64 
551.08 
536 
522.02 
506 
472.50 
454.31 
435.41 
415.50 
395 .03 


712.11 


(0-0) Band 

55667 .90 661.18 
667 .56 658.12 
666.56 654.18 
664.75 
662.19 
654.7 
621.77 
612.74 
602.90 
630.10 592.16 
621.95 580.55 
613.01 568.18 
554.91 
592.60 
581.19 
510.21 
542.04 476.24 
527.21 457 .94 
511.59 
495.18 418.96 
397 

(1-0) Band 

56756 
56763 .09 754.00 
750. 
746.75 
741.7 
755.72 
751.99 729.12 
721.44 
712.77 
692.87 
720.92 681.71 
712.11 669.54 
702.62 656.55 
692.05 
680.70 
612.11 
595.78 
641.48 578.50 
626.70 560.15 
611.13 541.05 
521.06 
577.06 500. 13 
539.74 455 83 


763. 
762.60 
760.90 


89 
55656 
668.66 
643.12 
659.71 
630.62 654.96 
649.44 
10} 614.93 643.12 
11} 635.77 
12} 595.97 
618.67 
14} 573.77 608.99 
15} 561.41 598.33 
16} 
534.26 574.49 
18} 519.46 561.41 
19} 503.81 547.34 
487 .36 532.52 
516.83 
223 452.00 500.27 
413.31 465.12 
348.98 
n 
iS 
56754.00 
751.15 
742.67 758.65 
730.89 751.15 
702.62 723.84 746.03 
692.35 715.92 740.04 
10} 681.04 733.14 
13} 642.30 675.97 707.19 
14} 627.77 663.74 
16} 595.78 673.19 
578.50 622.32 660.07 
18} 560. 58 606 77 646.10 
194 541.69 590. 42 631.24 
521.87 573.11 615.53 
21} 501.25 554.93 598.99 
234 516.10 563 08 
523.71 
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TABLE (Concluded) 


790.14 
57758.41 
10} 747.01 
11} 762.77 
12} 
13} 740.05 
691.78 
15} 675.81 713.78 
16} 
640.91 
18} 622.18 667.61 
19} 650.40 
581.93 632.30 
21} 560.46 613.26 
538.15 
514.80 
24} 
528.08 
26} 504.47 
480 .00 
28} 

850.93 
838.22 
830.51 
58798.97 
10} 787.01 812.33 
11} 774.10 801 
12} 760.14 
13} 745.54 778.13 
14} 729.87 764.86 
15} 713.11 750.85 
16} 694.03 735 
718.49 
18} 
19} 
662.64 
21} 642.17 


(2-0) Band 


829.12 


869.17 


775.40 
764.78 


713.11 
681.79 
664.78 
646.79 
565.89 
519.96 
495.61 
470.38 
444.24 
417.12 


815.22 
792.39 
779.59 
765.79 
750. 


57823. 


where A/B and terms involving 
plotted against the intercept the resulting curve 


the ordinate gives Using the data 
quoted Herzberg (1950) for the state, the spin splitting constant was 


have been neglected. If, therefore, 


820.75 829.69 
817.36 828.81 
813.00 827.15 
807.75 824.50 
822.35 801.71 820.97 
818.30 794.58 816.41 
813.41 810.61 
804.27 
788.65 
744.98 779.40 
765.43 718.69 758.15 
754.33 704.17 746.17 
742.25 
672.40 719.44 
715.50 655.14 704.55 
617.83 672.40 
650.99 576.91 
632.51 616.97 
613.26 
530.22 
506.19 
481.28 
(3-0) Band 
864.20 
841.48 860.39 
858.73 834.06 855.45 
825.67 849.48 
832.10 794.63 
816.25 
813.18 805.45 
802.41 754.91 793.77 
739.71 
725.59 
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TABLE VII 


The state The 
state 
v= v= v= 


calculated for the upper levels. The value used for the state was 
obtained from the levels and that the value for which 
required for the present work, may somewhat different. possible, there- 
fore, that small, constant, correction may required for the 
values the state. However, the relative values the A’s given 
Table VII are believed correct 0.20 Thus, while probably 
constant for the levels and seems definitely different for 
This borne out examination the plates which show that the 
relative positions the and heads are different the and (0-0) 
bands from the relative positions these heads the other bands. Fig. 
shows the effect the heads the and (3-0) bands. 
The constants for the upper state the new system are 


64781.56 


0.180 


The four bands which were analyzed have the level the state 
common. Consequently, fairly reliable value was obtained. Combining 
this with the data Jenkins, Roots, and Mulliken (1932) one derives values 
and which are probably more accurate than those hitherto available. 
The values obtained are follows: 


There are two features interest the new system. examination the 
plates shows that the over-all intensity the band about the same 
slightly greater than that the band. But whereas the branches the 
band are well developed normal way out fairly high values, 
there rather marked fading out the branches the (3-0) band inter- 
mediate The effect can seen Fig. and particularly obvious for the 
branches. the band, the intensity the branch drops rather 
suddenly from after which further lines are observed 
while the branch fades out rather more gradually from 
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The effect may due predissociation although the presence the 
(4-0) band presents difficulty this interpretation. However, should 
noted that the band very weak, and its presence would not too 
surprising the predissociation assumed weak forbidden. 

The second point interest anomaly the fine structure the (4-0) 
band. mentioned the previous paragraph, this band weak and further- 
more somewhat overlapped nitrogen bands. Consequently there was 
little hope making full rotational analysis. However, seemed that the 
two strong branches should found. The branch was, fact, readily 
detected and traced for about members. However, the branch was not 
found even though the positions the lines could calculated quite accu- 
rately. seems unlikely that could have escaped observation despite the 
weakness and overlapping the band and must concluded that the 
branch has anomalously low intensity. The reason for this not clear. 


VI. DISCUSSION 


was hoped that investigation the spectrum the Schumann 
region might give some information the dissociation energy the molecule. 
Unfortunately decisive result has been obtained. Nevertheless the results 
present some significant indications. 

the first place, should noted that the vibrational intervals for the 
state are quite small, and therefore seemed worth while carry out 
extrapolation using the data for the origins. linear extrapola- 
tion leads dissociation limit 87500 above the ground state. The 
state can arise from the atomic combinations 
and The energies the dissociation limits corresponding 
these combinations are given Table VIII, where each the three suggested 
values for D(CN), 50000 56000 and 66000 (see Douglas 
and Routly 1955 and Gaydon 1953), has been assumed turn. seen that 
the result the extrapolation closest agreement with the assumption 
that D(CN) 66000 and that the state arises from the configuration 
This interpretation has number satisfactory features. 
would mean that the two states, and observed for the lowest 
possible configuration which can give rise states, 
which does, fact, produce two such states. extrapolation for the 
state not possible only the level observed. However, the data 
available are least consistent with the requirement that both states 
the same atomic products. Thus the state has smaller internuclear 
distance than the state and should therefore expect more stable. 
point fact, F?A lies 5098.8 below the state. 

must noted that the extrapolation fairly long and 
consequently not high accuracy. slight modification the curve 
would bring the result into agreement with the assumption that D(CN) 
56000 and that the state goes the limit 
But this would require that there least one state other than 
which lies below and which has escaped detection that there 
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violation the non-crossing rule. Quite severe curvature the 
curve necessary bring the results the extrapolation into agreement with 
the other possibilities suggested Table VIII. therefore seems that 


TABLE VIII 
DISSOCIATION LIMITS FOR THE STATE 


D(CN) 50000 56000 66000 


there are abnormalities the potential curve for the state, the simplest 
and most satisfactory interpretation the present spectroscopic results 
that D(CN) 66000 and that the states and the configuration 

was suggested Section that the fading the structure the (3-0) 
band the system might due predissociation. Unfortu- 
nately even this yields little information D(CN). The energy 
which the effect takes place (68420 above the highest value suggested 
for the dissociation energy that the predissociation does not set useful 
upper limit the value D(CN). Furthermore, not possible make 
accurate estimate the energy the dissociation limit responsible for the 
predissociation. D(CN) 66000 would mean that the predissoci- 
ation was due the ‘II state arising from the ground state products 
Consequently the predissociation would forbidden which 
least consistent with the observation that the (4-0) band the J?A 
system present though weak. 

energy level diagram for shown Fig. which all the known 
levels and transitions for the molecule are indicated. The dissociation energy 
was taken 66000 and the dissociation limits arranged accordingly. 
The constants for the eight observed states are collected Table IX. 


TABLE 
CONSTANTS FOR THE KNOWN ELECTRONIC STATES 
1107.69 13.93 1.2948 1.419 0.0207 —25.1 
53954 8.78 1.155 1.503 0.013 
2143.7 20.2 1.962 1.153 0.023 
9118.02 12.61 1.7066 1.236 0.0174 —52.2 
1.175 0.01722 


1.8903 


13.114 


Douglas and Routly (1955) have listed the electron configurations and 
resulting doublet states that are expected for CN. The two new states, 
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Fic. Energy level diagram for CN. 


and almost certainly arise from the configuration 
which produces two two and two states. 
Thus, present, only four the states listed Douglas and Routly have 
not been directly observed—two states, II, state, and state. 
very likely that the latter responsible for one the perturbations mentioned 
Section IIT. 

The assignment and symbols was suggested Douglas and Routly 
because the absence the B—A transition seemed indicate that there was 
sufficient symmetry the molecule for properties important. The 
present work however would indicate that this not the case. Thus the 
“forbidden” transition appears quite strongly whereas the 
system was not observed. 

Further work remains done the spectrum the Schumann 
region. The spectrum should investigated the wavelength range below 
1300 where Rydberg transitions might observed which would give value 
for the ionization potential the molecule. experiment even more im- 
portance but greater difficulty would the observation the absorption 
observed and these transitions should give long upper-state progressions be- 
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cause the considerable increase internuclear distance going from the 
ground state. There some hope that this way accurate extrapolations 
dissociation limits might made and D(CN) obtained. 
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THE EFFECT GROWTH CONDITIONS UPON THE 


ABSTRACT 


This investigation provides both theoretical and experimental analysis 
the factors which affect the mode solidification binary alloy. These 
factors are: (i) the concentration solute the melt (Co); (ii) the rate 
solidification (R); (iii) the temperature gradient the melt the 
interface (G). Extremely high purity lead was produced zone refining and, 
from this material, crystals containing known concentrations tin, silver, and 
gold were grown under range well-controlled growth conditions. The mode 
solidification was investigated careful study the change appearance 
the interface with change growth conditions. For crystal 
containing specific solute was cbserved that (a) the transition from 
the width the cells varied inversely and inversely (c) the transi- 
tion from cellular interface dendritic interface exhibited large orientation 
dependence, and for constant orientation breakdown occurred critical ratio 
The experimental observations confirm both the existence solute- 
rich layer liquid adjacent the interface and its quantitative 
features. From this agreement with theory the diffusion coefficients tin, 
silver, and gold liquid lead 327°C. are determined. This work serves 
illustrate the effect extremely small amounts particular solutes upon the 
development substructures during solidification. technique proposed for 
obtaining measure the purity low impurity alloys. 


INTRODUCTION 


During the course the past years considerable effort has been spent 
investigating the substructures that occur during the solidification metals. 
Chalmers (1954) has recently given very excellent account these sub- 
structures that occur during the solidification metals from the melt. 
particular, considers dendritic growth pure metals, the macromosaic 
structure called the prismatic substructure known ‘‘corruga- 
and the topography the interface. From the knowledge 
gained these investigations now possible systematize the study 
solidification from the melt that insight can gained into the eventual 
prediction all the observed phenomena associated with the 
transformation. 

All phenomena observed occur result solidification from the melt 
can broadly classified various interactions between the three important 
variables growth: 

(i) K—the crystallographic features the growing crystal, 

(ii) 7—the temperature distribution the metal, 

S—the solute distribution the metal, 
that solidification can represented symbolically growth function 
where 
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(1) G(K,T,S). 


The prominent features crystal growth cannot studied independently 
but only the interaction between two the above variables. one considers 
the solidification pure metal, only the interaction between (i) and (ii) need 
considered. 

Consider, for example, the problem dendritic growth pure metals 
from the viewpoint expressed Eq. shown Weinberg and Chalmers 
(1951), necessary condition for dendritic growth pure metals that the 
temperature must fall from the solid into the liquid, and that the direction 
dendritic growth is, course, determined the crystallography the 
growing crystal. addition, Rosenberg and Winegard (1954) have shown 
that the rate growth dendrites determined the diffusion heat 
ahead the dendrite tip. Therefore, the dendritic structure which will 
obtained any given sample the result the interplay crystallography 
and temperature distribution during freezing. 

The terraced stepped interface structures observed 
Elbaum and Chalmers (1955) may also considered similar way. These 
structures occur when solidification proceeds under conditions such that the 
temperature rises from the interface into the liquid. The detailed appearance 
the interface structure depends very markedly upon the crystallographic 
orientation the interface. Again can seen that the structure depends 
upon the factors and 

The problem growth further complicated the existence solute 
the melt metal solidified. The degree complication will depend 
the particular kinds solute the melt. Generally, the more dissimilar the 
solute from the solvent, the greater its effect upon the solidification process. 
measure its effectiveness the ‘‘distribution the solute. 
This term given the symbol and defined the ratio the concentra- 
tion solute the solid the interface the concentration 
the solute the liquid adjacent it. For the solute creates more 
distortion the solid lattice than does the liquid, and may very 
difficult distinguish, for soldification purposes, the solute atom from the 
solvent atom the lattice. then the solute raises the melting tem- 
perature the solvent and thus increases the binding between the atoms 
the solid. The atoms are therefore packed closer together. which 
generally the case, the solute weakens the binding between the atoms and their 
interatomic distance increased. The solute atom distorts the lattice locally 
atom this type the lattice causes great deal distortion and will have 
great effect upon the solidification process. 

Solutes having small distribution coefficients will create large distortions 
when they become attached the interface and will make 
difficult for solvent atoms sit their regular adjacent lattice sites. 
likely then that the degree solute-produced imperfection introduced into 
the solid the solidification process will increase decreases. this way 
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can intuitively seen that the presence solute will have large micro- 
scopic effect the process solidification. 

The presence solute will also have several large macroscopic effects 
the properties metal crystal solidified from the melt; these can now 
considered. 

one considers the solidification melt single-phase binary alloy 
then the most consistent feature that appears result the freezing 
process end-to-end segregation solute the resulting ingot casting. 
This normal segregation one dimension and results from the different 
solubilities solute the solid and the liquid. The problem the resulting 
solute distribution the solid has been treated Pfann (1952), under the 
assumption complete mixing the melt. Tiller a/. (1953) have treated 
the same problem under the assumption convective mixing the melt, 
taking into account the slow diffusion rate solute the liquid, and arrive 
different resulting distribution. The assumption convective mixing 
introduced many new features into the solidification process, one which was 
the existence layer solute-rich liquid the solid—liquid interface. 
Wagner (1954) calculated the effect natural convection due this solute- 
rich layer and showed that general the effect natural convection was 
negligible except for very slow rates growth. arrived essentially the 
same conclusions Tiller (1953), that layer solute existed the 
liquid adjacent the interface which moved with the interface and that the 
distribution this layer was governed the rate growth the crystal 
and the rate diffusion the liquid. 

general, the normal segregation obtained any particular specimen 
depends upon the solute distribution which occurs allowed occur 
the melt during freezing and upon the speed freezing. Therefore, from the 
viewpoint Eq. this phenomenon may considered resulting from 
interaction the factors and 

The existence this solute distribution and its dependence the rate 
growth sufficient explain the phenomenon observed 
Stewart al. (1951). They grew single crystals lead containing polonium 
radioactive tracer and found taking autoradiograph the crystal 
that the polonium was segregated into bands across the width the crystal. 
These alternate dark and light bands the autoradiograph showed the length 
the crystal composed bands alternately high and low polonium 
content. Tiller (1953) showed that crystal solidified intermittent 
fluctuating rates from liquid alloy, bands solute appear the solid 
result the change growth rate, and furthermore that the growth tech- 
nique Stewart al. was such that the growth was intermittent. The de- 
tailed mathematical analysis this effect has been treated Smith, Tiller, 
and Rutter (1955) who showed how may used the production 
semiconducting devices. 

The effect results from the influence change speed 
solidification the solute distribution adjacent the advancing interface 
and can, therefore, classified interaction phenomenon. 
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The existence this solute distribution and its interaction with the tem- 
perature distribution the liquid the interface led Rutter and Chalmers 
(1953) explanation the substructure which was observed 
occur impure metal crystals grown from the melt under low temperature 
gradient the liquid and rapid growth rates. They showed that the bound- 
aries are regions rich solute compared the center the cell, and that the 
cell size varied with the growth conditions. Thus, can seen that for 
certain set growth conditions the solute segregation longer one-dimen- 
sional, but occupies three-dimensional cellular network the metal. 

Tiller (1953) have predicted theoretically that the breakdown 
macroscopically smooth interface (which gives one-dimensional segregation) 
interface consisting cellular projections (three-dimensional segrega- 
tion) occurs certain critical ratio temperature gradient, the liquid 
the interface the rate growth, follows 


(2) G/R —mC,(1—k)/Dk 


where the slope the liquidus line, the initial concentration 
solute the melt, the diffusion coefficient the solute, and the 
distribution coefficient. 

Walton (1955) have investigated the onset cellular growth for lead 
solute zone-refined tin, and found conform very accurately the 
predicted results Eq. Morris (1955) extended this study show 
that for another critical ratio these two parameters growth, interface 
the cellular type breaks down into interface which exhibits dendritic 
pattern; i.e., the projections develop side branches. 

Considering the foregoing investigations (Morris al. 1955; Rutter and 
Chalmers 1953; Walton a/. 1955), was felt that systematic and carefully 
conducted study should made provide extensive quantitative information 
concerning the stability and the significant features cellular and dendritic 
growth from the viewpoint expressed Eq. This work 
designed definitely establish the existence solute-rich layer liquid 
the interface and the quantitative features this solute distribution. 
provides theoretical and experimental analysis the three-variable problem 
exemplified the solidification binary alloy. 

The present work also attempt awaken full appreciation the 
the effect which very small amounts particular solutes have 
upon the development substructures that result from solidification from the 
melt. 


EXPERIMENTAL TECHNIQUES 
This experimental work was designed study quantitatively the effect 
varying amounts particular solute the mode solidification metal 
crystal. order carry out such experiments the following requirements 
must met: 
(1) Production extremely pure metal and subsequent alloying with 
particular solute. 
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(2) Solidification these alloys under carefully controlled conditions. 

(3) Measurement the growth conditions. 

(4) Observation the effect solidification. 
Purification and Alloying 

The solvent metal used this research was 99.965% pure Fisher reagent 
lead which was zone-refined (Pfann 1952) very high order purity. 
The zone refining was carried out in. diameter vycor tube which had 
been coated with carbon film prevent the lead from sticking it. This 
carbon film was produced burning acetylene the tube. The tube contain- 
ing the bar was then placed movable chassis shown Fig. This allowed 


Fic. Zone refining apparatus. 


the tube pass through induction heating coil constant rate order 
produce and move molten zone. small flow argon gas was passed 
through the tube order prevent oxidation the sample. The argon used 
was first purified passing over titanium sponge 850°C. 

Since lead melts 327°C., the natural cooling not sufficient maintain 
short molten zone; therefore, cooling coils were placed either side the 
induction coil extract enough heat produce molten zone reasonable 
length. order insure that there was time for complete mixing the 
liquid zone, the vycor tube was passed slowly through the induction coil. 
rate 0.3 cm./min. was found satisfactory. Each sample was given 
passes insure that enough purification had occurred. With respect 
solutes having 0.5, 0.01, 0.001 had estimated impurity content 
less than atomic respectively. 

The zone-refined lead was tested for impurity content anticipating the 
validity Eq. crystal the purified lead was grown under the most 
rapid growth rate 1.5 cm./min.) and the lowest temperature gradient 
2°C./cm.) that were attainable with the apparatus used these experi- 
ments. The decanted interface the crystal showed evidence cellular 
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structure even under these severe conditions. This indicates that the purity 
was better than that specified above. 

this starting material various amounts spectrographic tin were added 
lytically refined silver were also produced from weight. 
single alloy gold weight was also produced. 


Control Conditions Solidification 


large number crystals various solute contents were grown with 
range growth rates and under range temperature gradients the 
liquid. The crystals were not grown using the Chalmers’ technique (Chalmers 
1953) horizontally travelling furnace because this method has inherent 
disadvantages when one wishes vary independently both the rate growth 
and the temperature gradient the liquid. new technique was therefore 
developed which will referred the valve’’ technique. This 
will discussed detail below. These crystals were grown under neutral 
atmosphere eliminate oxidation the metal that clean outer surface 
would obtained for examination. critical stage during the growth 
the crystal the liquid was rapidly separated from the solid the decanting 
technique described Elbaum and Chalmers (1955). This provided means 
investigating the appearance the interface. The interface was 
studied for evidence cellular dendritic structure provide information 
concerning the stability different types interfaces under range 
growth conditions. 


Principles Thermal Valve Technique 


Consider cylinder metal, with efficient cooling coil soldered one 
end, positioned between two independently powered heating coils shown 
Fig. These coils may raised sufficiently high temperature melt 
some the metal and produce temperature distribution similar that 
shown Fig. the solid and the liquid. can readily seen that 
coil the ‘‘thermal was made decrease temperature, then the 
temperature gradient the solid behind the interface would increased. 
This increased rate heat extraction through the solid would produce the 
freezing some the liquid order satisfy the equation continuity 
heat flow given 


where thermal conductivity the solid, 
thermal conductivity the liquid, 
temperature distribution the solid, 
temperature distribution the liquid, 
latent heat fusion, 
distance measured from the interface origin. 


> 
> 
| 
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Cooling Coil Heating Coil Heating Coil 


Interface 


o, — Temperature x 


Distance 
(b) 


Fic. (a) Schematic illustration thermal valve technique. 
(6) Representative temperature distribution produced thermal valve. 


However, the interface moved forward and the temperature coil 
remains unchanged, then the temperature gradient the liquid will increase 
and growth will stop. therefore necessary lower the temperature the 
two coils appropriate rates that the crystal will solidify constant 
rate and constant temperature gradient will maintained the liquid 
all stages growth. This condition has been treated theoretically 
Appendix and can realized the temperatures the two coils and 
Fig. follow the cooling rates 


dT,/dt —GR 


and 


The furnace shown Fig. was built containing three heating coils wound 
ceramic tubes which form cage about the heating zone. These coils were 
powered three variacs capable delivering kva. each. Coils and 
were controlled one temperature, while Coil was controlled anothe 
temperature. 

grow crystal, the power distribution supplied the coils was first 
adjusted melt portion the sample and establish the desired temperature 
gradient the liquid. grow crystal any rate the temperatures 
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the coils were then lowered the rates predicted Eq. This was accom- 
plished feeding e.m.f. series with each the control thermocouples 
rate that, when translated into the temperature scale the thermo- 
couples, was the same that given Eq. This tended cause the tem- 
perature indicators the controllers scale this rate; however, 
the controllers decreased the power the coils order keep the tempera- 
ture indicating pointer the control point. Since the added e.m.f. plus the 
e.m.f. from the thermocouple keeps the controllers the control points, then 
the temperature the coils must fall rate compensate for that which 
e.m.f. added; i.e., the rate predicted Eq. 

The e.m.f. was supplied two gear-driven potentiometers which were 
connected shaft rotating constant speed. The rate delivery 
e.m.f. the respective thermocouple circuit was controlled setting the 
current the potentiometer windings suitable value. circuit diagram 
for the apparatus shown Fig. 

Single-crystal seeds were not used these experiments since they were 
found very difficult produce the proper orientation and the proper 
solute concentration the standard techniques. zone-refined lead, grain 
boundaries have extremely high mobilities, that stray nucleated 
may move rapidly into the rest the crystal. Seed crystals containing five 
ten crystallites were used this investigation and were found adequate. 


Measurement Growth Conditions 


The temperature gradient was measured with the aid six 
thermocouples placed cm. intervals the upper level the graphite 
boat shown Fig. The thermocouple leads were passed out the gas-tight 
system through fusite terminal cold junction Dewar flask filled with 
ice and water, and connected through rotating switch 
recorder. The rotating switch served connect each thermocouple turn 


Fic. thermal valve furnace with side bricks removed show windings. Graphite 
boat assembly used for growing specimens shown foreground. 
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Fic. Electrical circuit crystal growing apparatus. 


the recorder for period five seconds once every sec. The temperature 
gradient could calculated analyzing the plot the 
thermocouple readings (Morris 1955). The crystal was grown until the 
interface reached one the thermocouple positions and then the liquid was 
decanted. This allowed the temperature gradient that position the 
interface accurately measured from the slope this thermocouple plot. 

The rate growth the crystals was measured putting ripple marks 
the upper surface the crystal every minute half minute 
vibration the liquid. The distance between these time markers could 
measured accurately and, therefore, the rate growth accurately determined. 

The experimental features this work have been described more detail 
Tiller (1955). 

EXPERIMENTAL OBSERVATIONS 


large number samples were solidified under wide range growth 
conditions. was found that the manner which the solidification takes 
place characterized the appearance the solid—liquid interface during 
solidification; the interface may smooth (Fig. 9), may exhibit cellular 
structure (Fig. 17), may exhibit dendritic structure (Fig. 21). Therefore, 
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the experimental observations discussed refer largely the features 
the interface preserved for observation the decanting technique 
previously mentioned. All the photomicrographs appearing this section 
were taken under bright field illumination. 


Critical Growth Conditions Governing the Transition from Smooth Cellular 
Interface 

series lead crystals containing tin concentrations ranging from 
weight were grown under conditions the neighborhood 
the critical values predicted Eq. Crystals were produced over small 
range conditions for each concentration order determine, precisely 
possible, the conditions for which smooth interface first became unstable. 
series lead crystals containing silver concentrations ranging from 
weight and single alloy gold lead were similarly 
investigated. 

The results this series experiments have been tabulated Table 
for alloys, and Table for alloys and the lead-gold alloy. 


TABLE 
Co, G/R, Appearance 
wt. cm./sec. °C./cm. interface 

0.01 0.0167 5.4 Irregular cells 

0.01 0.0106 9.0 Nocells 

0.01 0.0167 Irregular cells 

0.01 0.0092 5.7 5.7X10? 

0.01 0.0106 4.35 4.1X10? Pox 

0.03 0.00735 8.2 Irregular cells 
and pox 

0.03 0.0110 8.7 Cells 

0.03 0.011 8.4 7.6X10? Cells 

0.03 0.00767 8.8 

0.03 0.00750 9.0 Nocells 

0.05 0.00416 4.0 Cells 

0.05 0.00334 5.4 Pox 

0.05 0.00867 13.8 Pox irregular 
cells 

0.05 0.00918 15.0 Pox 

0.06 0.0100 7.5 Cells 

0.06 0.00750 12.2 irregular 
cells 

0.06 0.00717 12.0 and 
irregular cells 

0.06 0.0102 8.2 Cells 


0.10 0.00700 13.5 Irregular cells 
0.10 15.0 and pox 
0.10 0.00567 12.0 Irregular cells 
and pox 
0.10 0.00835 7.5 Regular cells 


0.00567 13.0 Cells 

0.00284 13.0 Pox 

0.00333 13.0 Irregular cells 
0.00333 16.2 Pox and cells 


> 
q 
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TABLE 
Co, G/R, Appearance 

wt. cm./sec. °C./cm. interface 

0.0183 Elongated cells 

0.0167 4.7 Elongated cells 

0.0167 7.0 Pox 

0.020 4.2 Elongated and 
regular cells 

0.020 6.3 Regular cells 

0.015 6.8 Elongated cells 

0.007 8.5 Beginnings pox 

0.0067 10.8 Irregular cells 

0.008 9.0 Elongated cells 

0.0057 12.7 Smooth interface 

0.0055 10.8 Dense pox 

0.0050 10.5 Pox 

0.0037 10.7 Dense pox 

0.0046 12.5 Pox, elongated 
cells, irregular 
cells 

Co, G/R, Appearance 

wt. cm./sec. °C./cm. °C./cm.*/sec. interface 

0.020 5.3 cells 


Fic. Critical growth conditions governing the transition from smooth cellular 
interface for alloys. 
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The breakdown conditions for each concentration have been plotted Fig. 
for the alloys, and Fig. for the alloys. The criterion 
used establish the breakdown conditions for smooth interface was the 
general occurrence the pox-like structure shown Fig. 11. These points 
lie straight line predicted Eq. Eq. states that breakdown 


Fic. Critical growth conditions governing the transition from smooth cellular 
interface for lead-silver alloys. 


depends the ratio G/R and not these parameters independently. 
verify this, crystals were grown with tin solute content 
and breakdown was observed occur for the following conditions: 0.2 
15.0°C./cm. The rate growth and the temperature gradient 
have been changed factor three and yet the breakdown still occurs 
constant ratio These conditions cover approximately the total 
range gradients and speeds used these experiments, thus verifying that 
the important factor the ratio and not these parameters independently. 


Topography Cellular Interfaces 

was found possible follow the change from smooth interface 
cellular interface, and subsequently study the variation the appearance 
the cells with changes growth conditions. For convenience, conditions 
solidification will specified with reference Fig. stating the portion 
the diagram which represents them. 


o 
Cells 
= — 
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Crystals grown under conditions the right the transition boundary 
have smooth interface which exhibits the general microstructure investigated 
Elbaum and Chalmers (1955). There evidence instability with this 
type interface shown Fig. (there constitutional supercooling). 
the boundary approached, the interface develops the pox-like structure 
shown Fig. 10. With the aid dark field illumination was possible 
determine that the pox are projections the interface, and that they vary 
size from cm. width. They first appear widely scattered but 
the boundary approached their density and size increases shown 
Fig. 11. 

Just beyond the boundary the interface was observed develop mixture 
irregular cells and pox. conditions change away from the boundary into 
the cell region the pox tend disappear while the irregular cells increase 
both number and regularity. this region was observed that the irregular 
cells subdivide into elongated cells. Fig. illustrates interface exhibiting 
irregular cells and elongated cells. Note the density the pox the irregular 
cells compared their density the elongated cells. The elongated cells 
appear form abrupt faulting the interface produce boundary 
depression with its long axis certain crystallographic direction each 
crystal shown Fig. 13. Deforming the crystal produce multiple slip 
‘gave all the possible (111) traces the interface. Fig. typical example 
showing that the long axis the cell approximately parallel (111) 
trace. general, the long axis within +5° (111) trace the interface. 
However when the interface was (111) plane the long axis would slowly 
curve take range directions which were not necessarily (111) traces. 
Fig. shows interface elongated cells which are about 0.5 cm. long and 
0.0075 cm. wide. 

conditions change move further the left the transition 
boundary, the elongated cells subdivide into array nearly regular cells. 
This illustrated Fig. 16. Those crystals having their (100) axis perpendicu- 
lar the interface are the first develop array close-packed cells. 
The crystallographic variation this subdivision fairly small. 


Variation Cell Size and Shape with Conditions Solidification 

number specimens having concentration 0.15% tin lead 
were grown with various and order investigate the change cell 
size with these variables. was desired study the variation cell size 
for interface constant crystallographic orientation. was observed 
during preliminary experiments that only crystals with their axis the 
direction exhibited perfectly regular cells. Therefore, this experiment, 
only interfaces were considered which exhibited the regular hexagonal struc- 
ture. this way was possible certain that all cell size measurements 
were made interfaces the same orientation. 

decanted interface was obtained each sample and the average cell size 
about i000 cells was measured from photomicrograph such that 
represented Fig. 17. The results these experiments have been tabulated 
Table and plotted Fig. cell size vs. 1/RG. 
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Fic. Effect growth conditions upon the cell width (a). 
TABLE III 

wt. %Sn cm./sec. cm. 
0.15 4.65 15.1 
0.15 4.85 12.4 


0.15 7.0 6.2 


was found extremely difficult determine quantitative depen- 
dence cell size because crystallographic effects that arose during 
the investigation. was observed, however, that the cell size tended 
increase the concentration increased. 

very pronounced feature that was observed when crystals were grown 
under conditions lying the range which produced cellular interfaces was the 
effect the orientation the interface the shape the cell. mentioned 
previously, the regular cells shown Fig. occur only when the crystal 
axis close the dendrite axis. the orientation the crystal changed 
the appearance the cells also changed. Some crystals exhibit four- and 
tive-sided cells and the cell size may differ factor three four owing 
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the difference orientation the interface presented the liquid. This 
shown Fig. 18. 

When the growth conditions are made more severe, i.e., increased, 
increased, and decreased, the cells become less regular and they project 
further into the liquid developing definite cap-like appearance shown 
Fig. 19. These cells appear some cases project from the interface 
angle. likely that the axis the cell tends point the dendrite direc- 
tion and this direction makes angle with the interface normal, then the 
cells would not appear symmetrical. 


Critical Growth Conditions Governing the Transition from Cellular 
Dendritic Interface 


series lead crystals containing tin concentrations ranging from 0.25% 
2.0% weight were grown under conditions the region instability 
cellular interface. The appearance the interface was studied for evidence 
cellular dendritic growth. During growth the top surface the liquid 
was scraped clean, and after decanting the top surface was also studied for 
evidence which implied cellular dendritic type growth. 

Crystals were produced over small range conditions for each concentra- 
tion order determine, precisely possible, the conditions for which 
‘cellular interface first becomes unstable. These conditions will termed the 
“breakdown for the particular concentration question. The 
criterion used establish the onset dendritic growth was the general 
beginning side-branching the cells shown Fig. 20. should noted 
that this type dendritic growth somewhat different from the classical 
dendritic growth where the dendrites grow well out into supercooled liquid 
and develop pronounced side branches. The type dendritic growth discussed 
here found limited dendritic growth where the dendrite heads 
project about cm. from the interface. 

experimenter would faced with almost insurmountable task 
determine the breakdown conditions for every crystallographic orientation. 
However, was observed that the first crystal become unstable exhibited 
high index interface, i.e., its orientation falls the center the stereographic 
triangle standard (100) projection, and the last crystal break down 
exhibited the dendrite orientation. was therefore possible find the break- 
down conditions for these two limiting cases. This could done observing 
the conditions for which the cells any crystallite first became unstable, 
and the conditions for which the last cellular crystallite became unstable. 
was assumed that one the crystallites would have orientation which 
falls near the center the stereographic triangle, and another would have 
orientation near the dendrite orientation. general, the shape the cells, 
the orientation the interface relative the dendrite orientation could 
closely determined. 

The results this series experiments have been tabulated Table IV; 
the lower and upper breakdown limits have been plotted Fig. straight 
line can drawn through the points which represent the instability high 
index interface and second line through the points which represent the 
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Fic. growth conditions upon the transition from cellular dendritic 
interface. The upper and lower curves show the limits the orientation dependence. 
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TABLE (Concluded) 


Interface 


Co, 
wt. cm./sec. appearance 
0.75 7.5 Cells and dendrites 
0.75 9.6 Cells and dendrites 
0.75 12.0 139 *All cells (very 
close 
breakdown) 
0.75 5.6 43.1 **Dendrites, last 
region just 
breaking down 
1.0 12.0 208 *Cells, one region 
just breaking 
down 
1.0 12.0 176 Cells, few 
dendrites 
1.0 10.8 127 Cells, dendrites 
1.0 8.5 83.5 Dendrites, one 
region cells 
1.0 9.1 67.5 **Dendrites, one 
region 
unstable cells 
1.0 49.8 All dendrites 
103 Cells, unstable 
cells, dendrites 
region cells 
1.25 9.4 72.6 **Dendrites, one 
region 
unstable cells 
7.4 Dendrites, one 
region 
unstable cells 
1.50 10.4 84.4 Dendrites 
1.50 97.2 Dendrites, one 
region 
unstable cells 
1.50 10.9 91.0 Dendrites, one 
region 
unstable cells 
1.50 10.0 83.8 **Breakdown 
2.0 11.0 99.1 **Dendrites, last 
region cells 
breaking down 
2.0 10.2 83.4 All dendrites 
2.0 11.0 102 Dendrites and 
cells 
2.0 


*High index breakdown limit. 
index breakdown limit. 


instability interface the dendrite orientation. Between these two 
boundary lines the conditions are such that interfaces with orientations 
between these two extremes become unstable. The parameter was used 
for plotting purposes since was determined that the breakdown conditions 
for certain depend the ratio and not independently. 
Using Cy) 1.25%, the upper breakdown was observed for the following 
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PLATE 


Representative smooth interface. 

Interface exhibiting well-developed pox and irregular cells. 
Interface exhibiting irregular cells subdividing into 


ites 


PLATE 


Fic. exhibiting elongated cells two adjacent grains. (40X) 

Fic. Interface exhibiting elongated cells with traces duplex slip. The vertical slip 
traces are approximately parallel the direction the elongated cell boundary. 

Fic. exhibiting extremely long cells, 0.5 cm. X0.0075 cm. (40X) 

Fic. Interface exhibiting subdivision elongated cells nearly regular hexagonal 
(100X) 


III 


Fic. exhibiting regular cells. Representative photomicrograph used for 
measurement cell size. 

Fic. 18. Interface exhibiting effect crystallographic orientation upon the cell size and cell 
shape. (50X) 

Fic. 19. Interface well-developed cellular projections. 


Fic. exhibiting initial instability cellular interface resulting 
branching. 

Fic. exhibiting further developmert side-branching. 

exhibiting well-developed dendritic features. (100X) 

Fic. 23. Dendritic interface exhibiting pyramidal projections inclined the specimen 
axis. 


Fic. 24. Interface exhibiting the effect orientation upon the formation dendrites. 
Fic. 25. Approximately (100) type interface exhibiting dendritic 
Fic. 26. Approximately (100) type interface exhibiting truncated pyramidal projections. 
(100X) 


— 4 


Fic. Upper surface exhibiting initial formation dendrites evidenced the 
tion the corrugation boundary. 

Fic. Upper surface exhibiting stages dendrite formation adjacent grains. 


429 
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cm.; 0.4 cm./min., 5.8°C./cm. These conditions cover approxi- 
mately the total range gradients and speeds used the experiments dealing 
with the transition from cellular dendritic interface, thus serving 
indicate that the important factor the ratio and not these parameters 
independently. was also found that plotting versus gave linear 
plot. 


The Topography Dendritic Interfaces 

was found possible follow the changes from cellular interface 
dendritic interface, and subsequently study the variation appearance 
dendritic interfaces with changes the growth conditions. mentioned 
previously, crystals grown under conditions the right the lower break- 
down curve exhibit interface cellular nature shown Fig. 19. 
the boundary approached the square type cells, characteristic the 
orientation the center the stereographic triangle, become unstable and 
develop curved boundaries. The next stage the development that small 
regions boundary form the sides the cells shown Fig. 20. This 
the first evidence side branching the cells. the lower transition bound- 
ary crossed, the cell boundaries become greatly curved and develop further 
side branching (Fig. 21). conditions move into the region between the 
two curves, the dendrites, they should now called, develop very 
characteristic appearance (Fig. 22), with curved boundaries, side branches, 
and central projection extending considerable distance into the liquid. 
This illustrated best Fig. 23, which shows the projection extending into 
the liquid angle the specimen axis. 

shown Fig. 24, the breakdown has large orientation dependence. 
can seen that crystal one side the grain boundary exhibits 
well-developed dendrites, whereas the crystal the other side exhibits 
stable cells. 

the upper breakdown curve approached observed that only those 
crystallites having the dendrite orientation exhibit stable cells. Interfaces 
this type have been observed exhibit definitely characteristic micro- 
structure such that studied Elbaum and Chalmers (1955). the growth 
conditions approach the upper boundary the cell walls begin curve and 
breakdown occurs what appears the overlapping one the cells 
into another. This illustrated Fig. 25. the upper transition boundary 
reached, the cells develop truncated pyramid the center the cell 
which extends into the liquid. This illustrated Fig. 26. 


Topography the Free Surface the Transition Zone 
The upper surface the crystal was investigated function the growth 
conditions the transition region supplement the interface study. was 
found that the upper surface the crystal supplied additional information 
which could incorporated present three-dimensional picture the 
transition from cellular substructure dendritic branched substructure. 
Figure illustrates the appearance the corrugations the upper sur- 
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face. can seen that one crystal the corrugations intersect the upper 
surface angle, that the interface would present the unsymmetrical, 
large cells the liquid. another crystal the corrugations are running parallel 
the upper surface the crystal and, therefore, present the symmetrical 
cells the liquid. The corrugations these crystals appear perfectly 
constant size and shape, that the interfaces the crystals represented 
these corrugations would stable cellular form. 

However, when the interface tends become unstable, this observed 
the top surface oscillation the corrugation boundaries shown 
Fig. 28. Fig. illustrates the crystallographic effects dendritic breakdown 
the upper surface. can seen that the crystal with corrugations inter- 
secting the upper surface large angle exhibits well-developed dendrites, 
whereas the crystal with the corrugations running parallel the upper 
surface has only just broken down. 


SUMMARY OBSERVATIONS 


(1) The breakdown macroscopically smooth interface cellular 
interface exhibits the linear variation A(G/R) predicted Eq. 
where the constant proportionality. was determined that the ratio 
G/R was the important parameter this relationship, and not inde- 
pendently. 

(2) The above relationship exhibits small crystallographic dependence. 
The breakdown conditions are slightly different for interfaces different 
orientation, both crystals containing the same 

(3) the transition region, the interface develops pox-like structure 
stage development prior cell formation. 

(4) The cells first appear irregular and elongated just inside the 
transition zone. the cellular region approached, the elongated cells 
subdivide form smaller, almost regular, cells. 

(5) The direction the elongated cell boundaries provides information 
concerning the crystallographic orientation the crystal. The direction 
the elongated cell boundary coincides with (111) trace the interface. 

(6) The variation cell size for the dendrite orientation, function 
shown have the form B(1/RG), where the cell size and 
the constant proportionality. The cell size observed increase with 
increase 

(7) The magnitude the crystallographic effect increases the cellular 
region. Interfaces exhibit four-sided and five-sided cells well hexagonal 
cells. The cells for one orientation may three four times larger than the 
cells for another orientation. The high index interface exhibits the enlarged 
cells. 

(8) The breakdown from cells dendrites depends upon the crystallo- 
graphic orientation the interface. Crystals having orientations which fall 
the center the stereographic triangle standard (100) projection 
break down first. Breakdown for these high index interfaces can described 
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the linear relationship where the constant pro- 
portionality. Crystals having the dendrite orientation break down last. 
Breakdown for this low index interface was found exhibit the linear varia- 
found that was the important parameter these relationships, and not 
independently. 

(9) this transition region the cells project further into the liquid than 
previously. The cell boundaries curve and develop side-branches they 
break down into dendrites. the cells break down into dendrites they develop 
pyramidal projections extending into the liquid. 

(10) Observations were made the top surface the crystal which supple- 
mented and corroborated the interface observations. 


DISCUSSION 


The experiments carried out the smooth interface cellular interface 
transition serve endorse the validity the theoretical considerations 
concerning the stability smooth interface. Thus, the authors have definitely 
established that solute-rich layer liquid adjacent the interface exists, 
and that the distribution solute this layer exhibits the quantitative 
features predicted Tiller a/. (1953), i.e., 


where represents the solute distribution the liquid, and refers 
distance measured from the interface into the liquid. 

This work has also served verify that there possibility producing 
constitutionally supercooled zone liquid adjacent the interface, and 
that this unstable zone exists the interface will alter shape attempt 
eliminate the instability. From the viewpoint expressed Eq. can 
seen that the development cellular interface with its subsequent three- 
dimensional segregation mainly due the interaction two the impor- 
tant variables growth and The effect very slight for the break- 
down smooth interface cellular interface. This small dependence 
probably due the slightly different surface area presented the liquid 
crystals different orientation (different surface topography), shown 
Fig. 

The experiments show that although the rate solidification changes 
order magnitude the breakdown smooth cellular interface can 
described constant distribution coefficient Therefore, the rate 
deposition atoms from the liquid onto the solid for the rates considered 
does not appear affect the fundamental processes involved equilibrium 
freezing exemplified the equilibrium diagram. This reasonabie 
conclusion since even for the slow rates growth used this investigation 
the net rate deposition atoms line perpendicular the interface 


Figures and show that the breakdown smooth interface cellular 
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interface described linear relationship between and the ratio G/R 
predicted Eq. From the slope the transition boundary the diffusion 
coefficient for the solute liquid lead the melting temperature lead can 
calculated. The results this calculation have been tabulated Table 


TABLE 
Solute Co, G/R, 
°C./wt.% wt.% 


The accuracy these calculations depends very largely upon the accuracy 
the parameters and The experimental determination the breakdown 
conditions can made accurate within 5%, whereas and must taken 
from equilibrium diagram data for low concentrations and assumed 
unchanged for dilute solutions. 

measurement the diffusion coefficients tin silver lead are 
available for comparison; however, they fall the expected range. Jost (1952) 
value can extrapolated the melting temperature lead using the 
activation energy for viscous flow the activation energy for diffusion. 
Using 2.5 kcal./mole the activation energy, the value diffusion coeffi- 
cient becomes cm.?/sec., which quite good agreement with 
this measured value. 

often found very difficult estimate accurately the slope the 
liquidus line for dilute solutions; however can eliminated from Eq. 
utilizing the relation which relates the differences the slopes 
the liquidus and solidus lines the equilibrium diagram the latent 
heat fusion and the absolute melting temperature the solvent 
that 


Eq. then becomes transformed 


where the universal gas constant, and measured mole fractions. 

Equation may tested using the values and tabulated 
much 50% for these values and However, Eq. though not 
accurate Eq. very useful for purposes. 

has been shown Walton al. (1955), this technique may used 
give measure the purity any metal. Consider the following example. 
Assume that for lead crystal the growth conditions were 
1°C./cm., and that the decanted imterface had smooth appearance. 


(7) G/R 
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Let approximately cm.?/sec. for all possible solutes. Then 
inserting these values Eq. can stated that the lead purer than 


atomic fork 0.1, 
atomic fork 0.01, 
atomic fork 0.001. 


The lead will purer than indicated above since breakdown had not 
occurred under these conditions. This provides analysis the impurity 
content terms the distribution coefficient; chemical analysis can 
arrived grouping the elements into their respective ranges Since 
the ratio constant for any crystallographic system (entropy fusion), 
can seen that the higher the melting point the solvent the greater 
the order the purity that can analyzed for. 

hoped that these results have served illustrate the importance 
very minute quantities particular solutes gold lead) 
producing substructures result the solidification process. 


Interface Topography 

The presence the pox-like structure the region the transition zone 
between smooth and cellular interfaces appears the first manifestation 
the attempt, the interface, redistribute solute order eliminate 
small amount constitutional supercooling. The first pox appear the 
interface are probably chance fluctuations conditions which are able 
produce local regions constitutional supercooling. the general conditions 
approach the instability, the density the pox increases. likely that when 
the amount solute redistributed becomes too great the pox are re- 
placed the cellular structure. The density the pox the irregular cells 
observed much greater than the elongated cells. When the amount 
cell boundary increased there need for the pox and they disappear. 

all considerations dealing with the size and the shape the cells 
necessary make two assumptions: (i) The cell will try eliminate the 
constitutional supercooling providing solute sink into which the excess 
solute can diffuse. (ii) The crystal all times will try minimize its free 
energy; this end will try minimize its interfacial energy minimizing 
its solid—liquid surface area. The cell width will therefore try large 
the growth conditions will permit; this reasonable since the smaller the cell 
width for constant cell depth the greater the interfacial energy the 
system. 

When the elongated cells form they appear faulting the 
interface. They not appear development from the pox but seem 
quite separate attempt relieve the constitutional supercooling 
forming solute sink. Contrary the information Blaha (1955) these 
boundaries have orientation characteristics. The formation these 
elongated cells increases the interfacial area; however, the increase inter- 
facial energy minimized since the cell boundaries are traces (111) planes 
and the cell walls are thus low energy (111) surfaces. 
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When the amount constitutional supercooling eliminated increased, 
the elongated cells subdivide provide more cell boundary into which solute 
can diffuse. Fig. shows the gradual distortion the elongated cell 
produce more regular cells, and then the abrupt faulting produce sharp 
boundary these points. 

Figure represents the experimental results which illustrate the variation 
cell size, function the growth parameters and These results 
show that the cell width varies inversely the temperature gradient and 
inversely the rate growth agreement with the work Rutter and 
Chalmers (1953). The cell width inversely proportional since the rate 
movement the interface governs the time allowed for the diffusion 
solute which must redistributed produce stable cellular interface. 
The cell width inversely proportional since the temperature gradient 
governs the extent which the cell projects into the liquid. The further 
does the cell project into the liquid the greater the lateral diffusion gradient 
solute, that the solute can diffuse further the time allowed for redistri- 
bution. The longer the projection, the wider can and still redistribute 
the solute. Therefore, since the length the projection inversely propor- 
tional toG, then also the width. 

was observed that the size and shape the cells bear marked depen- 
dence upon the orientation the interface. This anomalous effect probably 
dependent upon the interfacial energy the particular orientation considered. 

Rutter and Chalmers (1953) have shown that crystals having their 
dendrite axis angle the growth direction are grown under conditions 
which fall the middle the cellular range, then the corrugations align 
themselves the growth direction. However, the growth conditions are 
made more severe, the corrugations try align themselves along their 
dendrite direction. This same observation was made the present investiga- 
tion. the conditions approach the upper limit the cellular range, the 
cells become asymmetrical the interface and the projections exhibit 
axis which appears the dendrite axis the crystal making angle 
the interface normal. The traces these inclined projections the upper 
surface the crystal (Fig. 27) further substantiate the observation that the 
cellular shape for crystals having nondendritic orientations asymmetrical 
with respect the interface normal. probably this asymmetrical shape 
the cell that makes unstable and break down dendrites before the 
symmetrical cell. very likely that the intercell region, because its 
asymmetry, becomes unstable for the same reasons but under less severe 
conditions than for the symmetrical intercell region. 


Breakdown Cellular Interface toa Dendritic Interface 


Figure shows that the lower and upper breakdown curves are linear when 


plotted the form 


where and are the proportionality constants. Therefore, 
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will exhibit cells and dendrites; 


(iii) the interface completely dendritic. 


has been shown that both increasing the solute concentration the melt 
and decreasing the temperature gradient tend lengthen the projections 
the cellular interface. has also been shown that increasing the rate 
growth tends decrease the cell width. Therefore, decreasing the gradient 
and increasing both the rate growth and the solute concentration the depth 
width ratio the cells can increased. immediately apparent that 
the further the projections extend ahead the rest the interface, the more 
solute will there the base, and therefore the greater will the chance 
that will become constitutionally supercooled. This effect enhanced 
the fact that the longer the cells become with respect the width the more 
will the isotherms become distorted this intercell region. must remem- 
bered that latent heat being evolved the surface the cell solidification 
proceeds. Some this heat will diffuse the intercell region and tend 
flatten the temperature gradient. Both these effects will tend produce 
constitutional supercooling this intercell region. The cells cannot decrease 
width lower the concentration the base the cells, and thus eliminate 
the instability, because this would only produce greater instability. When 
constitutional supercooling produced this region the cellular interface 
will unstable. This instability will reach critical value for certain ratio 
depth width the cell, thus causing the onset dendritic growth was 
observed. increased the projection will become longer well wider. 
It. was observed that the increase projection length much greater than 
the increase cell width; therefore, for the same and the ratio depth/ 
width the cell increases with This indicates that increased the 
ratio does not need great produce breakdown was observed 
experimentally. 

These results can compared with the work Morris al. (1955) the 
same system. They studied the breakdown from cells dendrites observing 
the conditions under which the corrugations the upper surface began 
develop noticeable side-branches. They observed only the upper surface, and 
considered only corrugations which ran parallel the axis the crystal, 
that they were effectively observing the upper breakdown limit. Their results 
were plotted the ratio G/R with the emphasis being placed this ratio 
rather than the parameter was observed this investigation. However, 
was observed during the course this investigation, the observations 
the upper surface provide much coarser criterion for interface instability 
than observations the interface itself. For this reason the half power 
dependence the rate growth could not probably have been deduced. 
the results Morris al. are plotted versus they produce 
wide scatter about the upper breakdown limit. can said, therefore, that 
there real disagreement between the results these workers and the 
present investigation. The present investigation provides much more 
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complete, and therefore more accurate, study the instability cellular 


interface. 

can concluded from this work that the phenomena associated with the 
solidification binary alloy are essentially caused interactions between 
the three important variables growth. The crystallographic variable 
seen increase importance from insignificant effect for macroscopically 
smooth interface very pronounced effect for the onset dendritic growth, 
whereas the temperature distribution and the solute distribution have 


pronounced effect throughout. 
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APPENDIX 


constant temperature gradient the liquid lowering the temperatures 
the coils and shown Fig. 2(a) specified rates. The gradients the 
solid and the liquid the interface must maintained constant during 
growth for this occur. Consider the following schematic diagram: 


where the temperature coil 
the temperature the interface, 
the interface from the beginning the specimen, 
the total length the charge, solid and liquid. 


From the diagram can seen that 
can seen that for constant rate Rt. Therefore, 
The continuity heat flow must satisfy Eq. the 
LpR the interface. 
Therefore 
and 


find the rate which and must decrease order have con- 
with respect time. Therefore, order solidify crystal constant 
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rate under constant temperature gradient the temperature coils 
and respectively, must decrease the following rates: 
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OBSERVATIONS STARS AND HEAVY PRIMARIES RECORDED 
EMULSIONS FLOWN VIKING ROCKET NO. 


HERMAN YAGODA 


ABSTRACT 


Methoas for the preparation fresh nuclear emulsions for use the study 
cosmic radiation rocket exposures are described. Shattering and mechanical 
damage minimized substituting transparent plastics for the glass support. 
casting the emulsion temporary microporous medium such plaster 
Paris, emulsions mm. thick can dried about six days without appreciable 
loss sensitivity. Strong, thin-walled housings for the emulsions are described. 
The star production rate emulsions flown the mesosphere peak altitude 
217 km. 1180+280 and 1.45+0.17 times greater than similar 
measurements the stratosphere (30 km.). The augmentation nuclear 
disintegrations can accounted for the opening the effective solid 
angle for star-producing radiation during the rocket flight. Energetic stars 
with three more shower particles record with frequency 278+30 
about three times greater than balloon flight the same geomagnetic 
latitude. For large showers the N/P ratio about 0.17. The observa- 
tional ratio nuclei based spectrum 238 heavy primary tracks 
When corrected for fragmentation surrounding matter this ratio 
reduced 2.5+0.5. The omnidirectional flux heavy primaries 
the top the atmosphere 7.62 analysis the 
zenith angle distribution the heavy primaries shows that the unidirectional 
flux fair agreement with geomagnetic theory for zenith angles <60°. The 
directional flux for 60° about two times greater than predicted theory. 


INTRODUCTION 


High altitude rockets have been employed extensively the study cosmic 
radiation above the sensible atmosphere, employing telemetering systems for 
transmitting the data during flight ground level receiving stations (Newell 
1953). The counter telescopes, usually sensitive over small solid angle, 
provide poor statistics during the brief period high altitude flight, and the 
complex apparatus must considered expendable nature. While nuclear 
emulsions would appear ideal for these studies because their small weight, 
large sensitive volume, and their ability record tracks incident from all 
directions, comparatively little work has been done with them, largely be- 
cause the necessity physical recovery the emulsions after flight termina- 
tion. 

Following the first recovery plates from V-2 rocket (Yagoda al. 
1950) series successive failures were encountered which terminated with 
the firing Aerobee NRL-10 May 1952. Study the emulsion recovered 
from this flight 2050 thick, cm. disk made from gel) showed 
star intensity prongs) and the passage nine 
heavy primary nuclei. 

This experiment suggested that the recovery larger areas emulsion 
important information could secured the charge spectrum the heavy 
primaries reaching the top the atmosphere, and that possibly definitive 
evidence could secured for the abundance lithium, beryllium, and boron 


received August 1955. 
Contribution from the National Institutes Health, Bethesda, Maryland, U.S.A. 
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nuclei the primary beam. was also difficult reconcile the high intensity 
stars recorded the mesosphere, about six times greater than that observed 
stratosphere balloon flights conducted the same geomagnetic latitude, 
with theoretical expectations. 

shortcoming the nuclear emulsion brief rocket exposures its 
continuous sensitivity between manufacture and development. The rocket 
flight time about seven minutes’ duration insignificant fraction the 
total sensitive age which usually exceeds three weeks the plates are manu- 
factured abroad. For certain phenomena, such heavy primary tracks and 
large shower stars, the accumulated low altitude background does not inter- 
fere with the measurements. The low energy star background greatly reduced 
preparing the emulsions the launching grounds shortly before flight. 
This could done conjunction with the launching Viking rocket No. 
from White Sands, New Mexico, December 15, 1952. 


EXPERIMENTAL 


The emulsions were prepared from gel flown the launching area 
paper-insulated containers which kept the temperature below 20°C. the 
elevation White Sands (1230 meters, 41° N.) the star-producing 
radiation about three four times greater than sea level. order 
reduce the background during the drying the plates, the darkroom was 
located the basement massive concrete building. While awaiting flight, 
measurements were made star production the darkroom and free air 
outdoors. The respective rates were and 2.64+0.34 stars 
the two localities? showing that the heavy building served attenuate the 
star-producing radiation factor about 2.5. 

When the plates were dry—two six days, depending thickness—the 
loaded cameras were stored the control room adjoining the launching 
platform. This very massive structure concrete and steel for which 
the average absorption length estimated 1800 gm. Under 2000 
gm. clay George and Evans (1951) observed the formation 0.0066 
stars 55° This figure, normalized the ratio star 
rates above the absorber, suggests that the star production rate the control 
room was the order 0.015 

Preparations intended for use ground level were coated subbed-glass 
secured from Ltd. The preweighed glass was mounted circular 
plastic frame lubricated with thin film Dow-Corning Silicone stopcock 
grease. The molds were leveled and filled with molten gel 50+3°C. taking 
care remove air bubbles pouring through several folds moist cheese- 
cloth. When set the molds were covered with 1.5 mm. thick rigid sheet 
plaster Paris and dried circulating air 50% R.H. When nearly dry 
the molds were taken apart and the plates stored above potassium carbonate 


free air exposure can compared with independent measurement (Pomeroy 1951) 
made Albuquerque, New Mexico, elevation 1580 meters. Pomeroy found the 
intensity stars prongs) between 1.31 and 1.21 after correction for 
accumulated background. the same basis our White Sands measurements correspond 
1.25 +0.23 
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equilibrium with its saturated solution until they reached constant weight, 
+0.1 gm. This procedure not only reduces background eliminating the high- 
altitude airplane flight from England, but also provides precise measure the 
quantity emulsion deposited. For computing the volume and thickness 
the coatings density 3.80 gm. was employed, based experimental 
measurements emulsion dried under the above conditions. 

Several the glass-coated preparations were flown glass-to-glass sand- 
wiches together with thick film Duco cement. This procedure increases the 
rigidity the unit and the flexible bonding possibly reduced shattering. 
One these emulsions was recovered perfect condition, and 
second, while the glass shattered, the emulsion coatings were intact. 
general, best avoid the use glass backings and coat the emulsion 
transparent plastics. Disks Lucite Plexiglass are cut fit the camera 
(Fig. 1), are degreased with carbon tetrachloride and washed warm 


Fic. Cross-sectional view two-plate rocket camera. 

Bottom cup machined from solid disk Dural alloy 0.5 in. thick, provided with steel 
centering pin Subbed lucite carrying emulsion layers. Bakelite spacer ring. Layer 
sponge rubber. Aluminum cover secured cup means machine screws. 


detergent until the surface loses its water repellency. The disks are rinsed 
distilled water and dipped freshly prepared solution gelatin 
water 40°C. When dry the attached the mold with 
the aid Dow-Corning Silicone stopcock grease and filled with molten gel. 
This subbing procedure permits the emulsion dry uniformly without 
appreciable lateral creep for preparations whose final thickness less than 
600 The adherence the coating nonpermanent. The emulsion coating 
usually separates result impact the instrument section. still found 
adhering, the emulsion disk readily separated flexing the lucite, 
drying the emulsion above desiccant. 

For the preparation emulsion disks whose thickness 1000 employ 
microporous sheets plaster Paris the supporting base. This novel 
medium has the advantage that the gel can dry rapidly through both the upper 
and lower interface, thereby preserving the sensitivity the gel which ordin- 
arily tends deteriorate during long drying periods 25° 28°C. Also. 


B 
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after high speed impact the instrument section, the plaster appears take 
the shock landing, and the pulverized base, unlike shattered glass, does not 
cut into damage the emulsion. Residual adhering plaster, being slowly 
water soluble, readily removed during the early stages photographic 
processing. 

The plaster supports are prepared mold composed two sheets 
plate glass separated aluminum frame mm. thick which hole 
102 mm. diameter machined with slightly tapering walls. Forty grams 
plaster Paris added cc. 25% ethyl alcohol contained 100-cc. 
beaker which drop aerosol wetting agent has been added. This mixed 
thoroughly and air bubbles which reach the surface are destroyed directing 
stream ether vapor from inverted wash bottle. The alcohol increases 
the setting time, permitting opportunity for the ascension and destruction 
occluded air. When the mix reaches the consistency sour cream 
transferred the mold aiming for cone-like deposition. The second plate 
glass placed the apex the cone and gently depressed until reaches 
the aluminum frame. The plaster allowed set (two three hours) with 
the upper glass weighted down with kgm. mass. When completely dry the 
plaster disk subbed immersion solution containing gelatin and 
0.1% glycerin weight 40°C. dried constant weight relative 
humidity about 50%. The plaster disk attached the lucite-ring mold 
means molten paraffin wax. 

For purposes high altitude exposure the emulsions must placed 
light-tight housing sufficiently strong withstand stresses accompanying the 
impact the instrument Part the damage secondary character 
and caused direct hits instruments torn loose from their mountings. 
The shock landing per does not appear augment fog weaken the 
latent image. recovery the instrument section delayed, fading the 
latent image occurs, accelerated the high temperature prevailing the 
desert. order minimize fading hydrogen peroxide vapors the aluminum 
alloy cameras are coated with film paraffin wax, and all direct metal-to- 
emulsion contacts are avoided. 

Cameras with the greatest life-expectancy must have one member designed 
withstand large stresses, such the incorporation thick I-beam 
ring-structure. The main features the latter are shown Fig. The basic 
design can also modified house rectangular emulsion sheets suitable 
machining thick slab aluminum alloy (Dural ST). These housings 
introduce about gm. cm.~? absorber. The absorber can reduced 0.2 
gm. cm.~?, with increased risk mechanical damage, taping the emul- 
sion the rocket skin described Fig. The isothermal processing pro- 
cedure for the development mm. thick emulsions (Yagoda 1955) was 
employed all the emulsions recovered. 


impact velocity reduced about 300 mi. per hr. severing the instrument section 
from the main body means timed explosive charge when the descending vehicle 
about miles altitude (Rosen al. 1954). 
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Fic. for mounting emulsions Rocket skin for cosmic-ray exposures with 
minimum absorber. 

Rocket skin. Corrugated paper cardboard. Emulsion wrapped black paper and 
water-impervious Plastic backing plate, 1.5 mm. thick lucite. Layer 
white rubberized adhesive tape surrounding Layer rubberized 
for securing unit inside rocket skin. 


III. FLIGHT CHARACTERISTICS 


While the flight high altitude rocket brief compared with strato- 
sphere balloon exposure, the characteristics the essentially parabolic 
trajectory can established with high degree precision. For purposes 
comparing intensities cosmic radiation the mesosphere with those available 
from measurements the stratosphere, convenient consider the portion 
the rocket trajectory AKB Fig. when the vehicle above km. The 
time interval then constitutes the duration the mesosphere cosmic-ray 
exposure, which for this flight 407.5+ sec. 

The intensity events, such cosmic-ray stars, produced during the meso- 
sphere exposure secured the difference population equal volumes 
flown and control emulsions identical age and history. This differential has 
corrected for stars accumulated the flown material while ascending 
and falling from the km. level. These corrections are computed 
graphical integration and with the aid function star production 
given elevation, with atmospheric depth, For the White Sands 
region, our measurements indicate that emulsions star 
production adequately expressed, for correction purposes, the equation 


The correction for the ascending limb straightforward path 
usually known precisely from telescopic and radar measurements. For Viking 
this correction 0.15 stars (>3 prongs) With the aid drag curve 
secured the fall V-2 No. estimated that 0.29 stars accumu- 
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ALTITUDE 


Fic. trajectory Viking No. 


lated during the descent the instrument section over path BP. The total 
ascent and descent correction employed 0.45+0.1 stars 

The control plates remain the launching site and the flown material 
resides the impact point until recovery effected. the two points differ 
elevation third correction necessary. This can computed from the 
difference atmospheric depth between the impact and launching points and 
the elapsed time. and differ +300 meters, then the corrections 
evaluated from Eq. (1) are —0.56 and +0.24 respectively. the 
present exposure the impact point was only meters below the elevation 
the launching area and since recovery was effected less than hr., this 
correction was found too small warrant application. 


IV. STAR PRODUCTION THE MESOSPHERE 


The emulsions were scanned with 22X Leitz oil immersion objective and 
wide field oculars allowing 25% field overlap between successive 
traversals. Since the population stars and tracks was small, believed 
that detection efficiency close 100% was achieved. The flown and the 
control emulsions were scanned the same observer, using the identical 
criteria for small star acceptance. 

comparative histogram the star-size distribution the flown and con- 
trol group five-day-old emulsions shown Fig. Stars given number 
black and gray tracks are further differentiated into: 


(a) shower tracks (<1.5 minimum) absent, 
(b) stars having least one minimum ionization track, 


(c) stars initiated charged primary and exhibiting shower 
minimum ionization tracks. 


600 
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CONTROLS ROCKET FLOWN 


TRACKS ABSENT 
MINIMUM TRACK 


CHARGED PRIMARY 
Ns 33 


TOTAL No. STARS MULTIPLICITY 


Fic. Comparative histograms star size distribution flown and control emulsions. 
The control histogram composed background stars recorded 3.88 cc. the a-group 
emulsions, 7.65 cc. the and 1.26 cc. emulsion exposed the darkroom White 
Sands. The rocket-exposed histogram based observations 13.24 cc. freshly made 


The rocket-flown emulsions are characterized large fraction and 
(c) stars, and pronounced tail high energy stars 15. 

The relative position the emulsions shown Fig. and Table 
summarizes the star populations observed the several preparations 
similar age and position. This study shows essentially uniform star produc- 


TABLE 
VARIATION STAR DIFFERENTIAL WITH POSITION ROCKET 


Differential between 
flown and controls 


Position Star Emulsion Stars, 
count vol., cc. observed corrected 
a-Controls 3.88 8.75+1.5* 
Nos. 57, 70, 71, and 

Nos. 63, 64, 65, 

Average all 

a-emulsions 282 19.09 5.5641.4 


*The errors are statistical based the square root the star counts. 
probable error the estimation the differential star population was estimated 


as 


where the average the volumes rocket and control emulsions scanned and and 
represent the star populations per unit volume the exposed and control preparations, 


respectively. 
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Fic. disposition emulsions upper instrument section Viking No. 
All cameras were mounted vertically close proximity the rocket skin, which measured 
0.040 in. above line aa’, and 0.050 in. the section below it. The total mass above line 
bb’ estimated 304 Ib. with its center mass located 29.3 in. above These emulsions 
can considered embedded spherical space equivalent the volume the cone 
containing average density 0.21 gm. absorber and average absorption path 
11.5 gm. No. and were mounted near the motor section about 450 in. 
below the rocket tip the vicinity about 1500 Ib. metal. 


tion rate all levels the instrument section. Preparation No. located 
near the motor section close proximity large masses metal shows 
somewhat higher star intensity. Within the limits the statistical uncertainty 
all four estimates Table are consistent with average value AP, 
corrected for rise and fall the km. level, 5.56+1.30 the 
basis flight time 407 sec. above km., this differential corresponds 
average star production rate the mesosphere 
Two independent measurements secured from preparations which had aged 
while awaiting flight the vehicle (exposures and Table are fair 
agreement with this estimate. 

The several rocket and balloon flights made geomagnetic latitude 41° 
are compared Table II. While the rocket exposures yield consistently 
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TABLE 
COMPARISON ROCKET AND BALLOON FLIGHT STAR INTENSITIES 


Altitude km. Star 


Vehicle Date Exposure, intensity, 
Max. Av. sec. 


Rocket exposures White Sands, New Mexico 


No. 9-29-49 150.7 80.4 324 
Aerobee 128.3 65.4 287 4800 1000 
Aerobee 2-12-53 137.0 71.3 300 4460+ 1200 
Viking 12-15-52 217 124 407 1180+280 
1070+300 
732 3-27-52 33.1 22,800 
738 30.2 26,600 610+42 
740 5-22-52 31.0 19,800 930+73 
879 8-29-52 30.6 31,800 760+60 
Average balloon flight 815+104 


*The original value 5000+800 stars cited Yagoda al. (1950) has been 
revised conform more accurate trajectory for the flight. 

tLaunched from Goodfellow Air Force Base, San Angelo, Texas. The geomagnetic latitude 


higher. star intensity than the balloon observations, the three rocket flights 
are not good agreement with each other. The deviations among the rocket 
estimates may originate from the following factors: 

(1) Errors inherent securing the differential when the absolute star 
population the control plate high and the volume scanned relatively 
small. two-fold deviation the magnitude appears frequently 
individual plates small voiume. This sampling uncertainty may respons- 
ible for the high values indicated the V-2 No. and Aerobee No. 
flights. The emulsions exposed Aerobee No. were manufactured 
England and accumulated large background 74.7+4.1 stars the 
time development which made the determination especially difficult. 

(2) The instrument section the Aerobee rocket small 
loaded with equipment. Hence emulsions exposed close proximity the 
high density matter the Aerobees are more apt receive secondary 
star-producing radiation than those exposed the spacious Viking. 

(3) The sampling time rocket flights brief and individual readings may 
possibly reflect fluctuations the flux star-producing radiation reaching 
the top the atmosphere. 

Comparing the Viking star intensities with balloon observations about 
km. indicates that ratio 1.45+0.24 between the mesosphere 
stratosphere exposures observed. This increment largely originates from the 
following factors: 

(1) from the opening the solid angle within which star-producing radiation 
can reach the emulsion, 

(2) from secondary radiation impinging the detector from matter within 


the instrument section, 
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(3) from albedo particles originating the heavy motor section and reflected 
upwards into the emulsions some ft. above. 


The star increment secured from the Viking No. observations can 
employed estimate the flux the star-producing radiation J,° assume 
that the range the star-producing radiation emulsion the same 
magnitude observed condensed matter such ice (Harding 1949) 
aluminum (Bernardini a/. 1949) 200 gm. cm.~?, and ignore transition 
effects: 


(2) 


and 6.2 steradians. This flux can interpreted average 
unidirectional flux, and such can compared with the average unidirec- 
above the atmosphere rockets Van Allen (1953) and Gangnes al. (1949). 
For zenith angles <63° Perlow (1952) reports lower value for the aver- 
age flux 0.079+.005 The close agreement the flux 
magnitude, evaluated from the present emulsion study, with the rocket 
represents the correct normal that the higher estimates secured 
the early V-2 and Aerobee rocket exposures may erroneous represent 
abnormal fluctuations. 


STAR AUGMENTATION MECHANISM 


The effective solid angle incident radiation mean free path reaching 
the detector located atmospheric depth can computed from the 
Gross transformation: 

where the exponential integral xfe evaluated Rossi 
(1952). Since the vehicle constantly altering the altitude the detector 
average solid angle for the entire flight above computed 
numerical integration the trajectory terms 

The variation with atmospheric depth (in units shown 
Fig. The circles indicate independent estimates made from the experi- 
mental rocket measurements air masses function zenith angle 
reported Pressly (1953). These were computed numerical integration 


w/2 
0 


where the air mass from the vertical fixed altitude. The two methods 
are good agreement, but the direct air mass measurements yield solid 
angle about three per cent greater than estimated the Gross transformation. 
observationally large but finite value. 
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The variation during the flight Viking shown Fig. computed 
for collision mean free path high-energy nucleons air 110 gm. 
The value for the rocket flight above the km. level 6.22 steradians, 
whereas the effective solid angle subtended emulsions balloon flight 
about km. 4.83 steradians. can therefore anticipate star augmen- 
tation 6.22/4.83 1.29 result the greater number cosmic-ray 
primaries reaching the rocket plates. The measured ratio star intensities for 
the two exposures 1.45+0.17. The close agreement 
between these ratios, the basis the above simple model, indicates that 
the greater part the measured value originates from primary radiation 
and that only about per cent attributable secondary star- 
producing radiation created the rocket body and albedo effects. Measure- 


100 200 v0 400 Sec. 


Fic. Variation effective solid angle for the star-producing radiation with time after 
launching Viking rocket No. 
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ments the albedo (Winckler and Anderson 1954) 40° indicate 
that not more than five per cent the incident flux for particles 0.7 


any upward direction. 


VI. SHOWER STARS 


shown Fig. large shower stars are virtually absent the 
control plates and constitute 1.31 the fresh rocket exposed emulsions. 
therefore reasonable assume that this category stars recorded during 
the high altitude flight the vehicle, and that their frequency, will largely 
independent the age and low-altitude history the plates. total 
shower stars were observed the scanning 66.3 cc. flown emulsion. 
Their distribution according shower size and the charge incident particle 
shown Fig. Collectively, these energetic stars are produced the 


Fic. Histogram shower stars function number shower particles, The 
stars are classified accordance with the nature the particle initiating the disintegration. 
The black blocks represent events initiated heavy primary nuclei from which three 
more singly charged minimum ionization particles emerge. 


mesosphere with frequency About per cent these 
stars are jets associated with four less gray and black tracks 4), 
typical example which illustrated Fig. 13. 

For comparative purposes, 2.90 cc. volume emulsion flown the 
stratosphere (Flight 337 Table II) was carefully examined for shower 
stars using the same acceptance criteria the examination the rocket 
plates. These occurred with frequency indicating that 
stars with record about three times more frequently the rocket 
than the balloon exposure. noteworthy that this increment cannot 
explained its entirety the increase the effective solid angle 
rocket elevations. The large effect may originate from the small sample 
comparative shower stars available from the Texas balloon flight. The discre- 
pancy diminished comparing with the Bristol data (Fig. which are 
greater statistical significance (Camerini 1951). 


a 
= 
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MESOSPHERE, 


STRATOSPHERE, 


SHOWER STARS 


Fic. Integral shower size distribution for shower stars with function the 
number shower particles, 


VII. THE NEUTRON-PROTON RATIO 


Owing their instability, not likely that appreciable flux neutrons 
enters the top the atmosphere. However, high ratio neutron proton 
induced stars has been observed even the highest altitudes reached 
balloons. Since known that the incident primary radiation composed 
predominantly charged particles, this implies that the observed neutron 
component must created locally exceedingly rapid rate. therefore 
interest compare the N/P ratio (we define star initiated 
charged primary associated with three more shower particles) 
observed Viking exposure with available data balloon elevations. 

The data summarized Table III indicate that the N/P ratios for the 
rocket and balloon exposures 41° are indistinguishable within the 


TABLE III 
VARIATION THE N/P STAR RATIO WITH ELEVATION 


Integral value N/P ratio 


Exposure Maximum 
and Haber- 


hain (1954) ~30 0.14 


\ O 
O 
100 
-O 
O 
O 
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limits statistical error. The lowest N/P ratio reported that Engier and 
Haber-Schain (1954) who limited the observations primaries high energy 
(65 with The large number neutron induced shower stars 
observed the rocket flight can attributed the large mass associated 
instrumentation (~300 Ib.) proximity the emulsions. The present data 
(Fig. indicate that the ratio diminishes with increasing N,, and that 
when the ratio 0.17 for the rocket exposure. further discussion 
this ratio best postponed until the data Viking become available. 
this flight the emulsions were exposed the very tip the rocket with 
only the rocket skin and the camera mass directly above them, and the 
probability for charge exchange thus greatly reduced. 


VIII. OBSERVATIONS SLOW MESONS 


examining the emulsions for nuclear evaporations tally was also made 
events and o-stars that came rest the sensitive layers. correcting 
for starless negative mesons (1.35 and for pions that recorded low 
altitude, rough estimates the frequency pion termination and the 
ratio have been made, and compared with similar measurements balloon 
elevations. Table shows that the ratio mesons terminating 


TABLE 


COMPARISON OF SLOW MESON INTENSITIES AND PION 
CHARGE RATIOS ON ROCKET AND BALLOON FLIGHTS 


Pions terminating 


Exposure 


Balloons 55+8 


essentially the same the rocket and balloon emulsions and that negative 
predominate. The data suggest that the charged pion production 
rate times greater the rocket than the balloon exposures. This not 
unreasonable view the large mass the rocket and the known copious 
production mesons the interaction cosmic-ray primaries with heavy 
target nuclei. 

order reach the emulsions, the mesons observed terminate traversed 
the average ~12 gm. matter (see Fig. enclosed the instrument 
section. This suggests that the pions tallied were ejected from nuclear inter- 
actions with kinetic energy ~46 Mev. The local production the mesons 
further indicated the observation eight and two events 
which originated from stars within the emulsion. Collectively this corresponds 
the production 21+7 pion ejects during the rocket flight. 
set emulsions comparable thickness (Table II, flight 732) exposed 
internal slow meson production rate was 

The ratio slow pions Mev.) stars the rocket flight 0.35, 
which compared with 0.068 balloon flights the same geomagnetic 
latitude. Comparing these ratios with similar data evaluated Camerini 


q 


136 CANADIAN JOURNAL PHYSICS. VOL. 34, 1956 


al. (1949) shows that the rocket while flight serves efficient 
generator: 


Flight Slow meson-star ratio 
0.062 with lead 
0.040 without lead 


~3.4km., Jungfraujoch 0.11 


The scanning the emulsions did not reveal any tracks attributable 
small number the latter category may have 
been incorporated the tally, since mass measurements were 
made meson tracks terminating with the production single gray near 
minimum prong. 


IX. HEAVY PRIMARY SPECTRUM 


Measurements heavy primary tracks were confined chiefly thick 
emulsions (Table order provide statistically significant number 


TABLE 
HEAVY PRIMARIES RECORDED DURING FLIGHT VIKING No. 
Emulsion No. tracks Tracks 
Emulsion Station thickness, per flight 

450 2300 0.32 


*Emulsions and were mounted the sun follower and had fixed azimuthal orientation 
for the first 200 sec. flight. Beyond this point the spin the rocket proved too large for the 
operation the sun follower. 


along all measurable tracks. The linear density was established 
for using Koristka long working distance objective 
and eyepieces. primary charge calibration the density along 
minimum relativistic alpha particle tracks was measured using narrow 
angle beams such particles which originated from the collision heavy 
primaries. Since the rocket exposure was conducted 41° the earth’s 
magnetic field deflects nonrelativistic heavy primaries with minimum energies 
arrival less than 3.2 Bev. per nucleon. Hence, was assumed that all 
heavy primary tracks recorded were moving relativistic velocity and that 
their charge was given 


where and represent delta rays per 100 track (corrected for dip) 
along the heavy primary tracks and the alpha particle calibration tracks. 
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Owing the brevity the exposure and the freshness the emulsions 
correction for track-coincident background electron tracks was found 
too small warrant application. 

Because the small number heavy primary tracks recorded was not 
practical base the spectrum exclusively tracks which recorded small 
angles with the emulsion plane. Comparison histograms and Fig. 


NUCLEAR CHARGE, 


Fic. showing distribution measured heavy primary charges. 


shows that essentially the same charge distribution secured using only the 
long the steep tracks emulsions thickness quantitative comparison 


Ratio M/H 
Long tracks, 3.10 
Steep tracks, 2.45 


The principal advantage the long flat tracks resides the fact that the 
peaks adjacent light elements, such the and nuclei, appear 
resolved. all subsequent discussion the flux the data all tracks have 
been combined. 

the scanning the emulsions for tracks number tracks 
were observed which the basis their delta-ray densities might Li, 
quantitatively, and the members this group shown Fig. probably 
constitute small fraction the total primary and secondary flux. That 
some were picked under low power survey suggests that very few tracks 
escaped detection. Originally more careful independent survey for 
the was planned, but the course this work large block 
emuision was recovered from Viking No. which was flown the very tip 
the nose cone, virtually free from surrounding matter. the present 
experiment the emulsions were surrounded about 300 Ib. irregularly 
distributed instrumentation, and since this circumstance would make very 


STEEP TRACKS 
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difficult recognize primary flux among the secondary fragmentation 
products, was decided expend the effort the more favorable exposure. 
Further discussion the relative abundance Z-nuclei the primary 
cosmic-ray beam will postponed until the results from the new flight become 
available. 

well known, certain velocity range the tracks slow He, Li, 
Be, and nuclei ejected fragments from nuclear disruptions may 
confused with the tracks relativistic particles charge Owing 
the large thickness the emulsions very few spurious tracks are believed 
included the primary charge spectrum shown Fig. 11. most 


VIKING CHARGE SPECTRUM 


No. PARTICLES CHARGE 


NUCLEAR CHARGE, 


instances locally produced heavy fragments could observed either 
terminate within the emulsion layer, traced back star, and they 
crossed both emulsion surfaces, the range was usually adequate note the 
absence high energy delta rays. When relativistic heavy primary made 
collision the emulsion, any resulting heavy fast fragments were excluded 
from the spectrum. Likewise, the eight tracks recorded emulsion were 
omitted because the exposure favored the recording secondary particles. 
The following conclusions can drawn from the analysis the residual 238 
heavy primary tracks: 

The ratio medium heavy nuclei directly recorded during the 
flight Viking No. The ratio corrected for collisions 
residual air and surrounding instrumentation estimated 2.5+0.5 
the top the atmosphere. This value may compared with ratio 
deduced Noon and Kaplon (1955) from heavy primary collisions observed 
emulsions flown the stratosphere above White Sands. 


0 
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The medium group nuclei are composed chiefly carbon, nitrogen, 
and oxygen which appear diminish that order. This also agreement 
with the Rochester balloon observations (Kaplon al. 1954). 

Fluorine and neon constitute 12.6% all nuclei doubtful 
whether the delta-ray count sufficiently reliable distinguish between 
the two elements, and possible that some the tracks listed fluorine 
may represent the tail ends the oxygen and neon peaks. the fluorine peak 
given full significance its cosmic-ray abundance particularly striking 
view the very small astrophysical abundance this element. The rocket 
observations suggest that the ratio fluorine silicon tracks about three, 
whereas from Brown’s data (1949) astrophysical abundances the ratio 
0.009. small fraction the tracks identified fluorine may have 
originated locally the instrument section result the fragmentation 
heavier nuclei, but this would not account for the 300-fold discrepancy between 
the cosmic and astrophysical abundances. 

33/23 1.43. Eisenberg (1954) has computed value 2.1 for this ratio. 
general, the rocket exposure indicates larger abundance the very 
heaviest nuclei than usually observed balloon observations. This 
expected view the small collision mean free path air for heavy nuclei. 
Particles with charge constitute about seven per cent the 
flux. Proper identification this region very difficult owing the exceed- 
ingly high density delta rays, and the charge assigned individual tracks 
may error charge units. 


HEAVY PRIMARY FLUX THE TOP THE ATMOSPHERE 


Table shows that the number heavy primary tracks recorded per unit 
area emulsion varies with respect the position the detector the 
rocket. particular, casting No. located near the massive (1500 motor 
section exhibits flux about per cent smaller than the plates the instru- 
ment section. Plates and mounted the very tip the nose cone probably 
represent the closest approximation the true incident flux, the local 
absorber there minimum estimated gm. cm.~? (chiefly aluminum). 
The average uncorrected flux indicated these plates 
6.39+0.68 the assumption that the M/H ratio 
2.66, the corresponding observational fluxes for the medium and heavy nuclei 
are, the same units, 4.64and 1.75. 

Correcting for collision loss the aluminum instrumentation (Ay 37, 
gm. cm.~?) and the residual integrated air mass (Qy 6.12, 
6.11 steradians) yields the primary flux the top the atmosphere: 


2.16+0.5, 


Z>6 


This value compared Table with other estimates the heavy 
primary flux secured the same geomagnetic latitude ionization chamber 
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TABLE 
HEAVY PRIMARY FLUX TOP ATMOSPHERE 41° 


Particles 


Exposure Maximum 
altitude, km. M-nuclei H-nuclei 

Van Allan (1951) 

V-2 No. 101 
Peters, 

eynolds, and 

Kaplon, Noon, and 


flown the nose cone V-2 rocket No. Van Allen (1951) and with two 
stratosphere balloon flights White Sands carrying emulsions. The highest 
values the flux are reported Van Allen, but these are stated upper 
limits, and the pulse data may not have been corrected adequately for bursts 
induced nuclear disintegrations. The flux estimated from the Viking No. 
observations general somewhat lower than indicated the two balloon 
estimates. The largest uncertainty originates the correction for collision 
loss surrounding matter. the omnidirectional flux based all 238 
tracks and the average absorption path taken 11.5 gm. (see Fig. 
for basis this estimate), then the flux heavy primaries in- 
creased 8.1 which agrees more closely with the balloon 
observations. 
XI. ZENITH ANGLE DISTRIBUTION HEAVY PRIMARIES 

Zenith angles the 238 heavy primary tracks with respect the vertical 
rocket axis have been made measuring the angle between the track and 
the apparent vertical, and the dip angle plates exposed with the emulsion 
plane perpendicular the horizon, the zenith angle given 


With this geometric arrangement can shown that the directional flux 
related the number tracks the angular interval passing 
through the vertical plate area by: 


(6) N= 
where 


The results this analysis (Table VII) show that the unidirectional flux 
fairly uniform over the entire upper hemisphere. Only the 15° angular 
interval there large deviation from the mean, and this attributable 
statistical fluctuations only four tracks were observed within the interval. 

interest compare the unidirectional flux the heavy primary 
nuclei with predictions geomagnetic theory, particularly because, the 
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TABLE VII 
UNIDIRECTIONAL FI.UX HEAVY PRIMARIES 


Flux top 
atmosphere 


0.0236 16.9 +8.5 
15-30 0.1576 8.3 
30-45 0.3896 8.2 +1.5 
45-60 0.6576 6.4 +1.0 
60-75 0.8896 8.5 +1.0 
75-90 1.0236 6.9 +0.8 
0-30° 0.181 9.4 +2.3 
30-60 1.047 7.1 +0.8 
60-90 1.913 147 
0.5708 
45-90 2.5708 189 


rocket exposure, experiment and theory are not complicated the presence 
appreciable atmosphere. The large spin the vehicle precluded the 
study azimuthal asymmetry. However, integration over 360° 
azimuth, was possible compare the zenith angle distribution the 238 
heavy primary tracks with deductions from the 

point observation with the geomagnetic latitude the minimum 
energy arrival charged particles units entering zenith 
angles and azimuth given 


This Stérmer cone valid only for 40° 41° N.; larger zenith 
angles correction has introduced for the earth’s shadow. For 40° 
values combining the cone and the earth’s shadow were available 
from computations Alpher (1950) 40°, and these were corrected 
41° the launching locality interpolation. The Stérmer values were 
converted into cutoff energies for protons the relationship: 


which evaluates 


using 8.06 gauss for the earth’s dipole moment and the DuMond 
and Cohen (1953) values for the atomic constants. pointed out Vallarta 
(1948), the cutoff energy subject two small corrections, one taking 
account the distance between the earth’s magnetic center and the point 
observation; and angular correction for the angle between the magnetic 
vertical and the gravitational vertical. For the White Sands locality (a) 

available aspect data the rocket flight indicate that the longitudinal axis the 
rocket was oriented within +5° with respect the vertical. maximum deviation 23° 


occurred 150 sec. after launching result fuel burnout. This lasted for about sec. 
and the zenith angles small number tracks recorded during this period are erroneous. 
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(10) would tend increase the unidirectional flux about 0.2 per cent, 
and view its large uncertainty has not been applied. Likewise, the 
angular correction <3°, quantity well within the error our zenith 
angle determinations. 

equation (10) represents the distance from the dipole the point 
observation. White Sands the distance between the dipole and sea 
6.371 cm. which has increased the relatively small, but signifi- 
cant, distance reached the rocket. the peak elevation, the average 
elevation the rocket above the earth’s surface, assuming perfect parabolic 
trajectory, For the Viking No. exposure the final value 
used equation (10) estimated cm. For heavy primaries 


The unidirectional flux was evaluated combining equation (11) with 
the integral energy spectrum for heavy primary nuclei deduced Kaplon, 
Peters, Reynolds, and Ritson (1952): 


ANGLE 


Fic. Comparison observational zenith angle distribution (diamonds) unidirectional 
flux with predictions from geomagnetic theory (line). 


value was computed U.S. Army Map Service, from measurements 
the distance between the center and White Sands, New Mexico, viz. 6,372,238 
meters. 


4 
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fixed zenith angle, each converted into corresponding flux value, and 
these averaged over 360° azimuth, 

This function, based the Stérmer and shadow cones, plotted Fig. 
and compared with the experimental estimates the unidirectional flux. 
For zenith angles <40° the magnitude the computed and observed flux 
surprisingly good agreement. The largest deviation from geomagnetic theory 
occurs zenith angles >60°, region which our measurements have the 
smallest statistical uncertainty. 

combine our data into two large angular groups 45° and 
then would appear that the corresponding fluxes and are very rough 
agreement with the main features the computed geomagnetic function, 
since about per cent less than general, Fig. suggests ap- 
preciably larger flux 60° than can anticipated from geomagnetic 
theory. The correlation between theory and experiment would improved 
the following factors: 

the exponent the power law equation (12) was error ~10 per 
cent and lower limit was employed the computation. 

The application precise longitudinal correction factor for the White 
Sands region would reduce the magnitude the energy cutoff, and hence 
increase the unidirectional flux. 

has been suggested Simpson (1954) that the pertinent geomagnetic 
latitude that the point observation and not the one usually taken 
the surface. This shift the magnetic latitude coordinate might 
the order the very high elevations reached the rocket. 
assume that 44° points high above the atmosphere then the plateau 
flux increased 9.8, and large zenith angles the flux raised 5.0 
particles 


Directional studies the total charged cosmic-ray intensity above White 
Sands have been made other investigators employing rocket telemetering 
techniques. Experiments this nature have the advantage that the direction 
the counter telescope can correlated during flight with the aspect the 
vehicle means photographs the earth taken known time intervals. 
Measurements made Singer (1950), Gangnes (1949), and Perlow al. 
(1952) are consistent exhibiting high values the unidirectional intensity 
large zenith angles Singer attributes this albedo component 
about per cent the total charged particle intensity. Perlow and co-workers 
ascribe the increased intensity large zenith angles decreased shadowing 
the bottom their telescope the earth, permitting radiation enter 
the counter from the lower hemisphere. 

This effect not noticeable from our measurements the heavy primary 
component. interacting with air nuclei the relativistic secondaries from the 
fragmentation proceed downward and hence not contribute the splash 
albedo. Likewise, owing the large collision cross-section the heavy 
primaries with air nuclei, particles approaching from the lower hemisphere 
are destroyed collision travcising the large air mass gm. 
before reaching the emulsions even the peak the rocket trajectory. 
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XII. HEAVY PRIMARY COLLISIONS 
the examination the group thick y-emulsions flown Viking No. 
interactions heavy primary nuclei have been observed which the 


particle destroyed catastrophically, emerges from the collision with 
reduction charge >1. These events are described Table VIII, and while 


TABLE VIII 
CHARACTERISTICS HEAVY PRIMARY COLLISIONS 
Charge No. fragmentation products 


Event Emulsion incident 
particle 


their number statistically small, they provide interesting information 
the fragmentation tendency the first generation heavy nuclei entering 
the atmosphere. 

(1) three the events the primary gave rise heavy relativistic 
secondary proceeding the same direction, but reduced charge. particu- 
lar interest are events No. and which boron and fluorine nucleus 
are produced locally. This indicates that the charge spectrum compiled 
Fig. may adulterated with secondary particles originating from collision 
processes the rocket skin and from massive matter within the instrument 
section. the examination one track associated closely parallel beam 
relativistic alpha particles was observed (Fig. 14) clearly indicating that this 
carbon nucleus had originated distant collision process, and was excluded 
from the data. This secondary phenomenon, however, can only recognized 
when the singly and doubly charged satellite particles are fairly close the 
member and limits the elimination process fairly close-by collisions. 

(2) eight the events there evidence for fast meson production 
the number singly charged minimum particles observed must assigned 
fast protons order conserve the estimated charge the primary. 

(3) The ratio relativistic protons alpha tracks (excluding fast mesons) 
produced these collisions about the same for incident and nuclei 
and averages 48/13 3.7. 

(4) the basis the length heavy primary track surveyed rough 
estimate can made the collision mean free paths. The results are com- 
pared Table with the more statistically significant results reported 
Noon and Kaplon (1955) from balloon observations White Sands. The 


of 


4. 
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TABLE 


COLLISION CROSS SECTIONS OF HEAVY NUCLEI IN EMULSION 


Mean charge 
Exposure incident flux No. collisions gm. 


Stratosphere* 59.4+7.8 


Stratosphere 


*Based observations Noon and Kaplon (1955) balloon flight White Sands. 


available rocket data suggest somewhat larger collision cross section than 
that observed the balloon exposures, particularly the case the 
nuclei. 
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LETTERS THE EDITOR 


Under this heading brief reports important discoveries physics may published. These 
reports should not exceed 600 words and, for any issue, should submitted not later than six weeks 
previous the first day the month issue. proof will sent the authors. 


Recent work (Douglas al. 1955) the reaction indicates Q-value for this 
reaction —0.96 Mev. rather than +0.12 Mev. previously reported (Ajzenberg and Laurit- 
sen 1955). This value may combined with the Q-value for 
Mev. give the mass difference 9.732 Mev. The position the 3/2 state 
corresponding the ground state then calculated 11.79+0.10 Mev. when 
corrected for coulomb energy the manner suggested Peaslee (1954). the authors’ 
paper (Bartholomew al. 1955) the reaction broad state spin 1/2+ 
excitation 11.61 Mev. was reported which had the characteristics 3/2 state. The 
newly reported Q-value for the reaction now makes seem probable that this 
3/2 state corresponds the ground state rather than excited state 
was suggested the above paper (Bartholomew al. 1955). this interpretation correct, 
implies that with one neutron and two holes the proton p-shell, has ground state 
spin 1/2+ similar which has three particles. 
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The Role the 1.30 Mev. Gamma-ray the Decay 


Since 1949, large number papers have been published dealing with the A-capture 
component the decay 50-day particular, workers this laboratory (Johns, 
McMullen, Donnelly, and Nablo 1954; Johns, Waterman, and Cox 1953) reported 
gamma-rays the following energies and quantum intensities: 556+1 kev. (3.6 
57643 kev. 722+1 kev. kev. and kev. 
coincidence lead absorption experiment indicated that the 556 kev. radiation 
was cascade with both the 722 and 1300 kev. gamma-rays and the authors postulated 
level 1856 kev. account for these results. view the work Lu, Kelly, and 
Wiedenbeck (1954), who failed find any evidence for 556-1300 kev. coincidences, and the 
disturbing fact that our decay scheme could not reconciled with the results obtained 
studies capture gamma-rays (Kinsey and Bartholomew 1953; Motz 1953), was 
decided repeat both the energy and coincidence measurements with the more refined tech- 
niques developed during the last two years. 

The gamma-rays were studied external conversion large magnetic spectrometer 
(Johns and Nablo 1954; Johns, Waterman, Mac/ Askill, and Cox 1953), using mgm./cm.? 
uranium radiator place the mgm. radiator employed the earlier work and 
with much improved counting statistics. These experiments give value 1299+2 kev. 
for the highest energy radiation, excellent agreement with the earlier 
determination. However, the two weak gamma-rays energy 576 and 1271 kev. did not appear, 
upper limits for their intensities being set, respectively, and quanta per 
disintegration. 

The results the coincidence measurements are shown Fig. Two Nal scintillation 
gamma-rav spectrometers with single pulse-height were used coincidence 
circuit resolving time. Each counter was covered with 7.5 gm./cm.? lead elimi- 
nate pile-up created the strong 192 kev. radiation. Curve shows the pulse-height 
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Fic. Gamma-gamma and beta-gamma coincidence spectra 


distribution each the single spectrometers. The 1300 kev. peak much augmented 
this spectrum solid-angle pile-up. Curve shows the true height this peak obtained 
reducing the solid angle the value used for Curve estimate the intensity 
the 1300 kev. radiation based Curve good agreement with the value obtained 
external conversion. Curve represents the bremsstrahlung spectrum produced source 
pure beta-ray emitter with approximately the same end-point clear that 
the sharply-falling background Curve above 1.1 Mev. due bremsstrahlung caused 
stoppage Mev. beta-rays the lead absorber. 

The gamma-gamma coincidence spectrum obtained when the fixed channel was set the 
722 kev. photopeak and the window the other spectrometer was moved across the spectrum 
shown Curve The results obtained the same way with the fixed channel set the 
1300 kev. photopeak are presented Curve The data for the latter case have been corrected 
for small chance contribution which was point the spectrum greater than 0.02 counts/ 
min. These curves show clearly that the 556 kev. gamma-ray coincidence with the 722 
kev. radiation and not coincidence with the 1300 kev. radiation. Curve represents the 
spectrum 1300-556 kev. coincidences expected the basis the measured intensity the 
1300 kev. radiation and the decay scheme reported our previous work. clear that less 
than 1.5% the 1300 kev. gamma-rays can cascade with the 556 kev. radiation. 

The spectrum beta-rays coincidence with the 1300 kev. radiation was studied the 
same way the gamma-gamma coincidence spectrum described above. For this purpose, 
one the crystals was replaced anthracene crystal covered with mgm./cm.? 
aluminum prevent detection the abundant conversion electrons from the 192 kev. 
transition. The beta spectrum was scanned with kev. window while the fixed channel 


| 
j 


148 
556 
A 
722 
1300 
x100 
0.6 ot 
| 
5+ 2.18 
1+ 1.98 
556 2+ |p+ 
O 246+ 0+ 

| 


JOHNS DECAY 149 


remained focused the photopeak the 1300 kev. gamma-ray. The coincidence rates 
obtained were very low (~0.1 per minute 500 kev.) but were from times the chance 
rate and times the cosmic coincidence rate. eliminate any possibility that the observed 
coincidences were due scattering associated with the bremsstrahlung spectrum, comparison 
experiment was carried out with source about the same strength. The chance coinci- 
dence rate was nearly the same for but there were true coincidences the region 
above 150 kev. seems clear that the coincidences observed indium are due beta-rays 
feeding level 1300 kev. calculation the number expected the assump- 
tion that this group had the same intensity the 1300 kev. gamma-ray gave value within 
30% the number actually recorded, which represents good agreement view uncertainties 
the solid angles and efficiencies detection the counters. 

Fermi plot the coincidence spectrum, corrected for chance, presented Curve 
The spectrum appears well-behaved and has end-point 710+50 kev. Both the beta 
and gamma spectrometers were calibrated frequently during the month required build 
reasonable statistics, using sources the indium source position. Drifts the electronics 
over this period did not change the calibration more than kev. 

decay scheme proposed the basis the above experiments presented Fig. 
The spins the three levels have recently been uniquely determined through 
the polarization angular correlation experiments Brazos and Steffen (1955). Consequently 
the 1.278 Mev. level must fed from the state With the rejection the 1.856 
Mev. level originally proposed (Johns, McMullen, Donnelly, and Nablo 1954), this decay 
scheme now concordant with the results capture gamma-ray studies. K-capture transitions 
the 556 kev. level expected but seem weak. 

was originally pointed out Nussbaum and Van Lieshout (1953) that there should 
beta-transition excited state energy ~1.2 Mev. and they set upper limit 
its intensity particles per disintegration. From our experiments, clear that 
this group has intensity end-point 700 kev. Such transition has 
log value 5.5, characteristic allowed spectrum. Consequently, the 1.30 Mev. level 
represents the first excited state that nuclide and almost certainly has spin and even 
parity. 

recent communication from Dr. Motz the Brookhaven National Laboratory indicates 
that has arrived essentially the same conclusions through independent and somewhat 
different group experiments. 

The financial assistance the National Research Council Canada and the Research 
Council Ontario gratefully acknowledged. 
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Line Source with Progressive Phase Shift 


evaluating the energy radiation from acoustic line source with progressive phase 
shift, Thiessen (1955) obtained integral representation (his equation [2]). Instead the 
somewhat awkward development equation [4] the integral can readily expressed 
simple closed form containing tabulated functions, and furthermore asymptotic approximations 
are readily made for the cases where large. 

have quite generally 


which may evaluated for all values the anda 

When and both integers the first two sine terms vanish and the simple expression 
Thiessen this latter expression may used for equal zero any multiple 0.2, 0.1, 


| 

i 

| 

4 

i 

| 

| 

| 

| 

| 

| 


150 CANADIAN JOURNAL PHYSICS. VOL. 34, 1956 


and 0.05 respectively, and this provides sufficiently closely-spaced values derive the curves 
the type given Fig. 


For 
For 


Some discrepancies Thiessen’s curves near may noted comparison with this 
formula. These forms show the transition from the almost constant value radiated energy 
tow tothe rapidly decreasing region when exceeds unity. 
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Transitions Dilute Gadolinium Ethyl Sulphate 


The seven transitions dilute gadolinium ethyl sulphate have been 
extensively studied Bleaney, Scovil, and Trenam (1954) using single field modulation 
video paramagnetic resonance However, they did not report the detection 
any higher order transitions. The new double field modulation spectrometer developed 
Buckmaster and Scovil (unpublished) possesses sufficient sensitivity wide band operation 
(90°K.) with signal-to-noise ratios about 1000:1. The values the spectroscopic splitting 
factor (g) and splitting parameters for the spin-Hamiltonian developed Elliott and Stevens 
(1953) are given this letter for the transitions both the parallel and perpendicular directions 
the crystal symmetry axis. 

The transitions corresponding arise from the admixing the wave functions 
describing the energy levels due the largest off-diagonal term the spin-Hamiltonian 
(Elliott and Stevens 1953). Theoretically, this term not present when the direction the 
magnetic field parallel the hexagonal symmetry axis the crystal. This has been experi- 
mentally verified. The intensity these transitions this direction less than that detectable 
with the maximum sensitivity our spectrometer gm. 290°K.). 

and transitions are given. These measurements were performed the same crystal 
gadolinium ethyl sulphate magnetically diluted with the isomorphous, diamagnetic lanthanum 
salt the ratio about one part two hundred. The magnetic field values used for the 
calculations for the parallel direction from transitions were obtained extrapola- 
tion from measurements small angles each side this direction. 


TABLE I 
“*g'’ VALUES FOR DILUTE GADOLINIUM 


ETHYL SULPHATE 


(°K.) 4AM = +1 AM = +2 AM = +1 4M = +2 
290 1 9908 +0.001 1.9905 +0.001 


90 1. 9880+0 001 1.9902 +0. 002 1.9901 +0.001 1.9901+0.001 


The splitting parameters, which are defined the same manner Bleaney, Scovil, and 
Trenam (1954) and calculated from the same data used for Table are given Table II. 
The values these parameters for the perpendicular direction have been corrected for the 
effect the largest off-diagonal term 
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TABLE 
SPLITTING PARAMETERS FOR DILUTE GADOLINIUM ETHYL SULPHATE 
M = +2 +206 .2+2 —3.98+0.3 +0.4140.3 
+204.7+2 


Fig. shows oscillogram the six transitions when the magnetic field 
oriented with respect the crystal symmetry axis that the total separation these transi- 
tions about five hundred gauss. The low transition appears the right Fig. 


more complete study these transitions presently being undertaken. Measurements, 
intervals, the behavior these transitions function the angle between the 
direction the magnetic field and the symmetry axis the crystal are being made. 
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